不饱和类硅烯H2C=SiNaCl的密度泛函理论研究

采用密度泛函理论方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类硅烯H2C=SiNaCl的结构.结果表明,不饱和类硅烯H2C=SiNaCl共有4种平衡构型,其中非平面的p-配合物型构型能量最低,是不饱和类硅烯H2C=SiNaCl存在的主要构型.对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒.计算预言了最稳定构型的振动频率和红外强度.     

不饱和类硅烯H_2CSiNaCl的密度泛函理论研究

采用密度泛函理论方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类硅烯H2CSiNaCl的结构.结果表明,不饱和类硅烯H2CSiNaCl共有4种平衡构型,其中非平面的p-配合物型构型能量最低,是不饱和类硅烯H2CSiNaCl存在的主要构型.对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒.计算预言了最稳定构型的振动频率和红外强度.     

不饱和类硅烯H2C=SiNaF的DFT研究

用密度泛函理论方法, 在B3LYP/6-31+G(d, p)水平上研究了不饱和类硅烯H2C=SiNaF的结构. 结果表明, 不饱和类硅烯H2C=SiNaF共有四种平衡构型, 其中非平面的p-配合物型构型能量最低, 是不饱和类硅烯H2C=SiNaF存在的主要构型. 对平衡构型间异构化反应的过渡态进行了计算, 求得了转化势垒. 计算预言了最稳定构型的振动频率和红外强度.     

类卡宾H2C=CLiF的构型及其异构化

本文用解析梯度和RHF/STO-3G方法研究了锂卤类卡宾H2C=CLiF势能面的主要特征, 得到了四种平衡构型及其异构化的过渡态构型. 平衡构型的能量按三元环构型、四元环构型、σ-配合物构型和p-配合物构型的顺序递增. 文中分析了各构型的特点, 给出了偶极矩、电荷分布和前线分子轨道, 并对该类卡宾的化学反应性质进行了讨论.     

不饱和类硅烯H2C＝SiMBr (M＝Li, Na)的DFT研究

采用密度泛函理论方法, 在B3LYP/6-311G (d,p)水平上研究了不饱和类硅烯H₂C＝SiMBr (M＝Li, Na)的结构. 结果表明, 不饱和类硅烯H₂C＝SiLiBr与H₂C＝SiNaBr各有三种平衡构型, 其中非平面的p-配合物型构型能量最低, 是这两种不饱和类硅烯存在的主要构型. 对平衡构型间异构化反应的过渡态进行了计算, 求得了转化势垒. 计算预言了最稳定构型的振动频率和红外强度.     

溴代类卡宾H2CMBr(M=Na,K)的结构与稳定性

本文用RHF/3-21G^*解析梯度方法研究了溴代类卡宾H2CMBr(M=Na,K)势能面的主要特征, 得到了它们各自的三种平衡构型及两种异构化的过渡态构型, 并用二级微扰方法进行了能量的相关能校正, 计算了振动频率, 以确证平衡构型和过渡态。计算结果表明, 溴代类卡宾H2CMBr(M=N,K))的三元环构型是最稳定的几何构型, 其它的两种构型σ-配合物和p-配合物, 也是势能面上的极值点, 但能量相对较高, 不稳定。本文分析了各构型的特点, 对其化学反应特性进行了评估。     

类硅烯H2SiClBeCl的结构及异构化反应

用DFT和QCISD方法研究了类硅烯H2SiClBeCl的构型与异构化反应. 结果表明, 类硅烯H2SiClBeCl有三种平衡构型, 其中四面体构型能量最低, 是其存在的主要构型. 对平衡构型间异构化反应的过渡态进行了计算, 求得了转化势垒. 计算模拟了最稳定构型的红外光谱.     

溶液中类锗烯H2GeClMgCl的结构与异构化反应

用密度泛函理论(DFT)和QCISD(quadratic configuration interaction with single and double excitations)方法研究了类锗烯H2GeClMgCl在气相和五种溶剂中的构型与异构化反应. 结果表明, 类锗烯H2GeClMgCl有三种平衡构型. 其中p-配合物型构型能量最低, 是其存在的主要构型. 讨论了溶剂效应对结构、能量与异构化反应的影响. 计算模拟了最稳定构型的红外光谱.     

溶液中类锗烯H2GeClMgCl的结构与异构化反应

用密度泛函理论(DFF)和 QCISD (quadratic configuration interaction with single and double excitations)方法研究了类锗烯 H2CeCMgCl在气相和五种溶剂中的构型与异构化反应.结果表明,类锗烯 H2GeClMgCl有三种平衡构型.其中p-配合物型构型能量最低,是其存在的主要构型.讨论了溶剂效应对结构、能量与异构化反应的影响.计算模拟了最稳定构型的红外光谱.     

不饱和类锗烯H_2CGeNaBr的密度泛函理论研究

应用密度泛函理论DFT方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类锗烯H2CGeNaBr的结构及异构化反应.结果表明,不饱和类锗烯H2CGeNaBr有3种平衡构型,其中非平面的p-配合物型构型能量最低,是其存在的主要构型.并对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒,预言了最稳定构型的振动频率和红外强度.     

Hypoelectronic dirhenaboranes having eight to twelve vertices: internal versus surface rhenium-rhenium bonding

Fehlner, Ghosh, and their co-workers have synthesized a series of dirhenaboranes Cp(2)Re(2)B(n-2)H(n-2) (n = 8, 9, 10, 11, 12) exhibiting unprecedented oblate (flattened) deltahedral structures. These structures have degree 6 and/or 7 rhenium vertices at the flattest regions on opposite sides of an axially compressed deltahedron thereby leading to Re═Re distances in the range 2.69 to 2.94 ? suggesting internal formal double bonds. These experimental oblate (flattened) deltahedral structures are shown by density functional theory to be the lowest energy structures for these dirhenaboranes. In some cases the energy differences between such oblate deltahedral structures and the next higher energy structures are quite considerable, that is, up to 25 kcal/mol for the nine-vertex Cp(2)Re(2)B(7)H(7) structures. The higher energy Cp(2)Re(2)B(n-2)H(n-2) structures are of the two types: (1) Most spherical (closo) deltahedra having unusually short 2.28 to 2.39 ? Re-Re edges with unusually high Wiberg bond indices suggesting formal multiple bonds on the deltahedral surface; (2) Deltahedra having one or two degree 3 vertices and 2.6 to 2.9 ? Re-Re edges. The latter deltahedra are derived from smaller deltahedra by capping Re(2)B faces with the degree 3 vertices.     

Structure, energy, and IR spectra of I2*-.nH2O clusters (n=1-8): a theoretical study

The authors report theoretical results on structure, bonding, energy, and infrared spectra of iodine dimer radical anion hydrated clusters, I(2) (-).nH(2)O (n=1-8), based on a systematic study following density functional theory. Several initial guess structures are considered for each size cluster to locate minimum energy conformers with a Gaussian 6-311++G(d,p) split valence basis function (triple split valence 6-311 basis set is applied for iodine). It is observed that three different types of hydrogen bonded structures, namely, symmetrical double hydrogen bonding, single hydrogen bonding, and interwater hydrogen bonding structures, are possible in these hydrated clusters. But conformers having interwater hydrogen bonding arrangements are more stable compared to those of double or single hydrogen bonded structures. It is also noticed that up to four solvent H(2)O units can reside around the solute in interwater hydrogen bonding network. At the maximum six H(2)O units are independently linked to the dimer anion having four double hydrogen bonding and two single hydrogen bonding, suggesting the hydration number of I(2) (-) to be 6. However, conformers having H(2)O units independently linked to the iodine dimer anion are not the most stable structures. In all these hydrated clusters, the odd electron is found to be localized over two I atoms and the two atoms are bound by a three-electron hemi bond. The solvation, interaction, and vertical detachment energies are calculated for all I(2) (-).nH(2)O clusters. Energy of interaction and vertical detachment energy profiles show stepwise saturation, indicating geometrical shell closing in the hydrated clusters, but solvation energy profile fails to show such behavior. A linear correlation is observed between the calculated energy of interaction and vertical detachment energy. It is observed that formation of I(2) (-)-water cluster induces significant shifts from the normal O-H stretching modes of isolated H(2)O. However, bending mode of H(2)O remains insensitive to the successive addition of solvent H(2)O units. Weighted average energy profiles and IR spectra are reported for all the hydrated clusters based on the statistical population of individual conformers at room temperature.     

Limited occurrence of isocloso deltahedra with 9 to 12 vertices in low-energy hypoelectronic diferradicarbaborane structures

Theoretical studies show that the 10-vertex system Cp(2)Fe(2)C(2)B(6)H(8) is the only one of the 2n skeletal electron Cp(2)Fe(2)C(2)B(n-4)H(n-2) systems (n = 9, 10, 11, 12) for which a true isocloso deltahedron having a single degree 6 vertex is highly favored over alternative structures. This is demonstrated by the occurrence of only the 10-vertex isocloso deltahedron as the central Fe(2)C(2)B(6) polyhedron in all nine of the Cp(2)Fe(2)C(2)B(6)H(8) structures within 8 kcal/mol of the global minimum. Low energy isocloso structures are also observed for the 11-vertex Cp(2)Fe(2)C(2)B(7)H(9). However, interspersed with these isocloso structures are Cp(2)Fe(2)C(2)B(7)H(9) structures based on deltahedra having two or more degree 6 vertices. For the 12-vertex Cp(2)Fe(2)C(2)B(8)H(10), the six lowest energy structures all have central Fe(2)C(2)B(8) deltahedra with two degree 6 vertices, one for each iron atom. The Cp(2)Fe(2)C(2)B(8)H(10) structures having a central Fe(2)C(2)B(8) icosahedron with all degree 5 vertices lie at significantly higher energies, starting at 17.8 kcal/mol above the global minimum. The 9-vertex Cp(2)Fe(2)C(2)B(5)H(7) system appears to be too small for isocloso structures to be favorable, although three such structures are found at energies between 5.5 and 8.0 kcal/mol above the global minimum. Five Cp(2)Fe(2)C(2)B(5)H(7) structures based on the tricapped trigonal prism lie in an energy below the lowest energy isocloso structure. The lowest energy Cp(2)Fe(2)C(2)B(5)H(7) structure and two higher energy structures within 8.0 kcal/mol of the global minimum have central Fe(2)C(2)B(5) deltahedra with a degree 6 vertex for each iron atom.     

Monte Carlo temperature basin paving with effective fragment potential: an efficient and fast method for finding low-energy structures of water clusters (H2O)20 and (H2O)25

Determining low-energy structures of large water clusters is a challenge for any optimization algorithm. In this work, we have developed a new Monte Carlo (MC)-based method, temperature basin paving (TBP), which is related to the well-known basin hopping method. In the TBP method, the Boltzmann weight factor used in MC methods is dynamically modified based on the history of the simulation. The states that are visited more are given a lower probability by increasing their temperatures and vice versa. This allows faster escapes from the states frequently visited in the simulation. We have used the TBP method to find a large number of low-energy minima of water clusters of size 20 and 25. We have found structures energetically same to the global minimum structures known for these two clusters. We have compared the efficiency of this method to the basin-hopping method and found that it can locate the minima faster. Statistical efficiency of the new method has been investigated by running a large number of trajectories. The new method can locate low-energy structures of both the clusters faster than some of the reported algorithms for water clusters and can switch between high energy and low-energy structures multiple times in a simulation illustrating its efficiency. The large number of minima obtained from the simulations is used to get both general and specific features of the minima. The distribution of minima for these two clusters based on the similarity of their oxygen frames shows that the (H(2)O)(20) can have different variety of structures, but for (H(2)O)(25), low-energy structures are mostly cagelike. Several (H(2)O)(25) structures are found with similar energy but with different cage architectures. Noncage structures of (H(2)O)(25) are also found but they are 6-7 kcal/mol higher in energy from the global minimum. The TBP method is likely to play an important role for exploring the complex energy landscape of large molecules.     

High-stability hydrogenated silicon-carbon clusters: a full study of Si2C2H2 in comparison to Si2C2, C2B2H4, and other similar species

The structural and electronic characteristics of the Si2C2H2 and Si2C2 clusters are studied by ab initio calculations based on coupled cluster and density functional theory using the hybrid B3LYP functional. In addition, similar species, such as SiC2H2 and Si3C2H2, are also studied for comparison. It is illustrated that the lowest energy structures of all three hydrogenated clusters, which have the general form Si(n)(CH)2, n = 1, 2, 3, are fully analogous to the structures of the corresponding organometallic isovalent carboranes. The most stable structure of Si2C2H2 is obtained by attaching two hydrogens onto the carbon atoms of a higher energy (+1.5 eV) planar trapezoidal structure of Si2C2, followed by geometry optimization which leads to puckering of the planar structure. Furthermore, it is demonstrated that Si2C2H2 and the other two "similar" hydrogenated clusters are much more stable than the corresponding bare nonhydrogenated clusters. Comparison of Si2C2H2 and C2B2H4 shows that their structural and bonding similarity includes also nuclear rearrangement similarity. The two species are isomerizable with an energy difference between their lowest energy puckered 1,2- and 1,3-isomers of about +/-0.3 eV. It is suggested that SiC2H2, Si2C2H2, and Si3C2H2 are special cases of a larger class of stable clusters. It is speculated on the basis of the calculated infrared spectrum that Si2C2H2 and perhaps other members of this class of clusters could be found in appreciable abundance in interstellar space.     

Density functional study of 8- and 11-vertex polyhedral borane structures: comparison with bare germanium clusters

Density functional theory (DFT) at the hybrid B3LYP level has been applied to the polyhedral boranes B(n)H(n)(z) (n = 8 and 11, z = -2, -4, and -6) for comparison with isoelectronic germanium clusters Ge(n)(z). The energy differences between the global minima and other higher energy borane structures are much larger relative to the case of the corresponding bare germanium clusters. Furthermore, for both B(8)H(8)(2-) and B(11)H(11)(2-), the lowest energy computed structures are the corresponding experimentally observed most spherical deltahedra predicted by the Wade-Mingos rules, namely the D(2)(d) bisdisphenoid and the C(2)(v) edge-coalesced icosahedron, respectively. Only in the case of B(8)H(8)(2-) is there a second structure close (+2.6 kcal/mol) to the D(2)(d) bisdisphenoid global minimum, namely the C(2)(v) bicapped trigonal prism corresponding to the "square" intermediate in a single diamond-square-diamond process that can lead to the experimentally observed room temperature fluxionality of B(8)H(8)(2-). Stable borane structures with 3-fold symmetry (e.g., D(3)(h), C(3)(v), etc.) are not found for boranes with 8- and 11-vertices, in contrast to the corresponding germanium clusters where stable structures derived from the D(3)(d) bicapped octahedron and D(3)(h) pentacapped trigonal prism are found for the 8- and 11-vertex systems, respectively. The lowest energy structures found for the electron-rich boranes B(8)H(8)(4-) and B(11)H(11)(4-) are nido polyhedra derived from a closo deltahedron by removal of a relatively high degree vertex, as predicted by the Wade-Mingos rules. They relate to isoelectronic species found experimentally, e.g., B(8)H(12) and R(4)C(4)B(4)H(4) for B(8)H(8)(4-) and C(2)B(9)H(11)(2-) for B(11)H(11)(4-). Three structures were found for B(11)H(11)(6-) with arachno type geometry having two open faces in accord with the Wade-Mingos rules.     

Ab initio surface reaction energetics of SiH4 and Si2H6 on Si(001)-(2 x 2)

First-principles pseudopotential calculations, within a simple dynamically constrained scheme, have been performed to investigate the reaction of 0.25 ML coverage of SiH4 and Si2H6 with the Si(001)-(2 x 2) surface. The silane molecule (SiH4) is adsorbed on to the surface at a number of different sites (on dimer, interrow, or intrarow) with varying barrier heights. Two distinct structures, which are similar in energy, arise from the initial dissociative reaction SiH4-->SiH3(silyl) + H, where the dissociated species are adsorbed either on the same dimer components or on adjacent dimer components. Several further decays of silyl from SiH4 are presented in two separate regimes of high and low ambient hydrogen coverages. The decomposition of silyl can form two different bridging structures: an on top or an intrarow bridging structure in both of the two hydrogen coverage regimes. The disilane molecule (Si2H6) is also adsorbed upon this surface with varying energy barriers, resulting in a dissociation reaction where two SiH3 species are adsorbed on one dimer or in an adjacent dimer configuration. Plausible energy reaction paths for the above models are presented. The stability of the SiH2 species is also discussed.     

Gold as hydrogen. An experimental and theoretical study of the structures and bonding in disilicon gold clusters Si2Au(n)- and Si2Au(n) (n = 2 and 4) and comparisons to Si2H2 and Si2H4

In a previous communication, we showed that a single Au atom behaves like H in its bonding to Si in a series of Si-Au clusters, SiAu(n) (n = 2-4) (Kiran et al. Angew. Chem., Int. Ed. 2004, 43, 2125). In this article, we show that the H analogy of Au is more general. We find that the chemical bonding and potential energy surfaces of two disilicon Au clusters, Si(2)Au(2) and Si(2)Au(4), are analogous to Si(2)H(2) and Si(2)H(4), respectively. Photoelectron spectroscopy and ab initio calculations are used to investigate the geometrical and electronic structures of Si(2)Au(2)(-), Si(2)Au(4)(-), and their neutral species. The most stable structures for both Si(2)Au(2) and Si(2)Au(2)(-) are found to be C(2)(v), in which each Au bridges the two Si atoms. For Si(2)Au(4)(-), two nearly degenerate dibridged structures in a cis (C(2)(h)) and a trans (C(2)(v)) configuration are found to be the most stable isomers. However, in the neural potential energy surface of Si(2)Au(4), a monobridged isomer is the global minimum. The ground-state structures of Si(2)Au(2)(-) and Si(2)Au(4)(-) are confirmed by comparing the computed vertical detachment energies with the experimental data. The various stable isomers found for Si(2)Au(2) and Si(2)Au(4) are similar to those known for Si(2)H(2) and Si(2)H(4), respectively. Geometrical and electronic structure comparisons with the corresponding silicon hydrides are made to further establish the isolobal analogy between a gold atom and a hydrogen atom.     

Structures, energetics, vibrational spectra of NH4+ (H2O)(n=4,6) clusters: Ab initio calculations and first principles molecular dynamics simulations

Important structural isomers of NH(4) (+)(H(2)O)(n=4,6) have been studied by using density functional theory, Moller-Plesset second order perturbation theory, and coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)]. The zero-point energy (ZPE) correction to the complete basis set limit of the CCSD(T) binding energies and free energies is necessary to identify the low energy structures for NH(4) (+)(H(2)O)(n=4,6) because otherwise wrong structures could be assigned for the most probable structures. For NH(4) (+)(H(2)O)(6), the cage-type structure, which is more stable than the previously reported open structure before the ZPE correction, turns out to be less stable after the ZPE correction. In first principles Car-Parrinello molecular dynamics simulations around 100 K, the combined power spectrum of three lowest energy isomers of NH(4) (+)(H(2)O)(4) and two lowest energy isomers of NH(4) (+)(H(2)O)(6) explains each experimental IR spectrum.     

Gitonic and distonic alkanonium dications (diprotonated alkane dications C(n)H(2n+4(2+)), n = 1-4)(1)

The structures and stabilities of gitonic and distonic alkanonium dications, i.e., diprotonated alkane dications C(n)H(2n+4)(2+) (n = 1-4), were investigated at the MP4(SDTQ)/6-311G**//MP2/6-31G** level. The global minimum energy structures (2, 4, 7, and 10) of the C(n)H(2n+4)(2+) dications are double C--H protonated alkanes to give structures with two two electron three-center (2e-3c) bonds. Two different dissociation pathways for the dications, viz deprotonation and demethylation, were also computed. Demethylation was found to be the favorable mode of dissociation.