杂硼原子簇XB6+ (X=C,Si, Ge,Sn, Pb)的稳定性和芳香性

利用从头算MP2方法和密度泛函理论B3LYP和B3PW91方法, 研究了杂硼原子簇XB₆⁺ (X=C, Si, Ge, Sn, Pb)的结构、稳定性及化学键合情况. 对C, Si, Ge, B使用6-311+G(d)基组, 对Sn和Pb使用LANL2DZ赝势基组. 研究结果表明, 具有C_s对称性的假平面XB₆⁺ (X=C, Si, Ge, Sn, Pb)结构是势能面上的全域极小点, 其稳定性要高于C_6v对称性的锥形结构和C₂对称性的假锥形结构. 在B3LYP水平上, 对这些异构体的势能面的极小点进行了自然键轨道(NBO)的分析; 对最稳定构型的最高占据分子轨道(HOMO)和最低空轨道(LUMO)能级差、分子轨道(MO)和核独立化学位移(NICS)进行了计算和讨论. 分析了杂原子和硼原子间、相邻硼原子间的键合情况, 讨论了最稳定构型的芳香性质.     

Density functional theory study of 11-atom germanium clusters: effect of electron count on cluster geometry

Density functional theory (DFT) at the hybrid B3LYP level has been applied to the germanium clusters Ge(11)(z) (z = -6, -4, -2, 0, +2, +4, +6) starting from eight different initial configurations. The global minimum within the Ge(11)(2-) set is an elongated pentacapped trigonal prism distorted from D(3)(h) to C(2v) symmetry. However, the much more spherical edge-coalesced icosahedron, also of C(2v) symmetry, expected by the Wade-Mingos rules for a 2n + 2 skeletal electron system and found experimentally in B(11)H(11)(2-) and isoelectronic carboranes, is of only slightly higher energy (+5.2 kcal/mol). Even more elongated D(3)(h) pentacapped trigonal prisms are the global minima for the electron-rich structures Ge(11)(4-) and Ge(11)(6-). For Ge(11)(4-) the C(5v) 5-capped pentagonal antiprism analogous to the dicarbollide ligand C(2)B(9)H(11)(2-) is of significantly higher energy (approximately 28 kcal/mol) than the D(3h) global minimum. The C(2v) edge-coalesced icosahedron is also the global minimum for the electron-poor Ge(11) similar to its occurrence in experimentally known 11-vertex "isocloso" metallaboranes of the type (eta(6)-arene)RuB(10)H(10). The lowest energy polyhedral structures computed for the more hypoelectronic Ge(11)(4+) and Ge(11)(6+) clusters are very similar to those found experimentally for the isoelectronic ions E(11)(7-) (E = Ga, In, Tl) and Tl(9)Au(2)(9-) in intermetallics in the case of Ge(11)(4+) and Ge(11)(6+), respectively. These DFT studies predict an interesting D(5h) centered pentagonal prismatic structure for Ge(11)(2+) and isoelectronic metal clusters.     

Theoretical investigation on the electronic and geometric structure of GaN2+ and GaN4+

The electronic and geometric structures of gallium dinitride cation, GaN2+ and gallium tetranitride cation, GaN4+ were systematically studied by employing density functional theory (DFT-B3LYP) and perturbation theory (MP2, MP4) in conjunction with large basis sets, (aug-)cc-pVxZ, x = T, Q. A total of 7 structures for GaN2+ and 24 for GaN4+ were identified, corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces, and bonding mechanisms for some low-lying electronic states. The calculated dissociation energy (De) of the ground state of GaN2+, X1Sigma+, is 5.6 kcal/mol with respect to Ga+(1S) + N2(X1Sigmag+) and that of the excited state, ?3Pi, is 24.8 kcal/mol with respect to Ga+(3P) + N2(X1Sigmag+). The ground state and the first excited minimum of GaN4+ are of 1A1(C2v) and 3B1(C2v) symmetry with corresponding De of 11.0 and 43.7 kcal/mol with respect to Ga+(1S) + 2N2(X1Sigmag+) for X1A1 and Ga+(3P) + 2N2(X1Sigmag+) for 3B1.     

Theoretical study of gallium nitride molecules, GaN2 and GaN4

The electronic and geometric structures of gallium dinitride GaN 2, and gallium tetranitride molecules, GaN 4, were systematically studied by employing density functional theory and perturbation theory (MP2, MP4) in conjunction with the aug-cc-pVTZ basis set. In addition, for the ground-state of GaN 4( (2)B 1) a density functional theory study was carried out combining different functionals with different basis sets. A total of 7 minima have been identified for GaN 2, while 37 structures were identified for GaN 4 corresponding to minima, transition states, and saddle points. We report geometries and dissociation energies for all the above structures as well as potential energy profiles, potential energy surfaces and bonding mechanisms for some low-lying electronic states of GaN 4. The dissociation energy of the ground-state GaN 2 ( X (2)Pi) is 1.1 kcal/mol with respect to Ga( (2)P) + N 2( X (1)Sigma g (+)). The ground-state and the first two excited minima of GaN 4 are of (2)B 1( C 2 v ), (2)A 1( C 2 v , five member ring), and (4)Sigma g (-)( D infinityh ) symmetry, respectively. The dissociation energy ( D e) of the ground-state of GaN 4, X (2)B 1, with respect to Ga( (2)P) + 2 N 2( X (1)Sigma g (+)), is 2.4 kcal/mol, whereas the D e of (4)Sigma g (-) with respect to Ga( (4)P) + 2 N 2( X (1)Sigma g (+)) is 17.6 kcal/mol.     

(RXNR)_4(X=B，Al，Ga)簇合物的结构与化学键性质

用自洽场理论（HF）和密度泛函理论（DFT）的B3LYP方法,在6-311G~*水平上 ,首次从理论上研究了(RXNR)_4 (R = B,Al,Ga;R = H,CH_3,NH_2,OH)簇合 物及其先驱化合物（RXNR）_2的几何构型、电子结构和化学键性质,并与其异构体 及其等电子化合物进行了比较。结果表明,(RBNR)_4为环状骨架结构,(RAlNR)_4 和(RGaNR)_4为建立骨架结构。     

Mechanisms of [1,3]-Sigmatropic Migrations of the Nitroso Group in the ON-X-CH=X Systems (X = O,S, Se,NH, CH<Subscript>2</Subscript>)

[1,3]-Sigmatropic migrations of the nitroso group in the systems ON-X-CH=X (X = O, S, Se, NH, CH₂) were studied by MP2(fc)/6-311+G** and B3LYP/6-311+G** quantum-chemical calculations. The energy barrier in the process was estimated at 2.4 (2.5), 20.0 (25.0), and 22.3 (23.4) kcal/mol for X = O, NH, and CH₂, respectively. The energy minima for X = S and X = Se correspond to cyclic structures with two-coordinate NO group, which are more stable than acyclic structures by 9.3 (4.3) (X = S) or 13.1 (5.7) kcal/mol (X = Se).     

Theoretical study on structures and stabilities of N4X (X = O,S, Se,Te) series

The stable and transition structures of N₄X (X = O, S, Se, Te) series with singlet state are optimized with the ab initio (MP2) and density functional theory (B3LYP) methods using the 6‐311+G(d) basis set. The ring isomers are found to be the global minima for N₄O, N₄S, N₄Se, and the chain isomer is the minimum for N₄Te. The stabilities are studied by evaluating the dissociation barriers with respect to dissociation. The reactants and products connected by transition structures are determined by applying the intrinsic reaction coordinate (IRC) calculations. The C_2v, C_3v and ring isomers decompose into linear NNX and N₂ molecules, however, the chain isomers decompose into cyclic N₂X and N₂ firstly. A new possible isomerization mechanism between the cyclic and linear structures of N₂X series is studied. The cyclic structures of N₂X convert into linear structures easily with the very low barriers. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

内含式复合物X@(HAlNH)12 (X=Be,Mg, Ca,Zn, Al+, Ga+)的结构和稳定性

在研究闭式多面体(HAlNH)12簇合物几何构型及稳定性的基础上, 用DFT的B3LYP方法在6-31G(d)的水平上, 对其内含式复合物X@(HAlNH)12 (X= Be, Mg, Ca, Zn, Al+, Ga+)进行了构型优化和能量计算, 并讨论了稳定结构的几何构型、自然键轨道(NBO)、振动频率、能量参数及NMR数据与结构的关系. 用Gaussian 03的QST3方法确定了客体X通过笼面6-元环的迁移过渡态(TS)结构, 并用IRC方法对所得TS结构进行了验证. 最后得到内含式复合物X@(HAlNH)12结构在热力学和动力学上的稳定性信息, 其中复合物Ga+@(HAlNH)12的结构相对最稳定.     

Microsolvation of hydrogen sulfide: exploration of H2S.(H2O)n and SH-.H3O+.(H2O)n-1 (n = 5-7) cluster structures on ab initio potential energy surfaces by the scaled hypersphere search method

The scaled hypersphere search method was applied to ab initio potential energy surfaces of the H2S.(H2O)n/SH-.H3O+.(H2O)n-1 system with n = 5-7. Local minima databases including 121, 326, and 553 structures for n = 5-7, respectively, were obtained based on calculations at the MP2/6-311++G(3df,2p)//B3LYP/6-31+G** level. In these small cluster sizes, the SH-.H3O+.(H2O)n-1 type is still unstable relative to the H2S.(H2O)n type, and the global minima for H2S.(H2O)n are very similar to those of pure water clusters of (H2O)n+1. Thermodynamic simulations based on the present databases showed a structure transition from the well-mixed (H2O)n+1-like global minimum at low temperatures to unmixed complexes between H2S and (H2O)n at high temperatures.     

Structure and bonding in cyclic isomers of B2AlH(n)(m) (n = 3-6, m = -2 to +1): a comparative study with B3H(n)(m), BAl2H(n)(m) and Al3H(n)(m)

The structure, bonding and energetics of B(2)AlH(n)(m) (n = 3-6, m = -2 to +1) are compared with corresponding homocyclic boron, aluminum analogues and BAl(2)H(n)(m) using density functional theory (DFT). Divalent to hexacoordinated boron and aluminum atoms are found in these species. The geometrical and bonding pattern in B(2)AlH(4)(-) is similar to that for B(2)SiH(4). Species with lone pairs on the divalent boron and aluminum atoms are found to be minima on the potential energy surface of B(2)AlH(3)(2-). A dramatic structural diversity is observed in going from B(3)H(n)(m) to B(2)AlH(n)(m), BAl(2)H(n)(m) and Al(3)H(n)(m) and this is attributable to the preference of lower coordination on aluminum, higher coordination on boron and the higher multicenter bonding capability of boron. The most stable structures of B(3)H(6)(+), B(2)AlH(5) and BAl(2)H(4)(-) and the trihydrogen bridged structure of Al(3)H(3)(2-) show an isostructural relationship, indicating the isolobal analogy between trivalent boron and divalent aluminum anion.     

Theoretical study of aromaticity in small hydrogen and metal cation clusters X (X=H,Li, Na,K, and Cu)

Five cation clusters X (X=H, Li, Na, K, and Cu) with two possible isomers, i.e., regular trigonal structure (D_3h) and linear structure (D_∞h), have been investigated using four methods: B3LYP, B3PW91, MP2, CCSD(T) and basis set 6‐311+G(3df). The calculations show that only the regular trigonal structure (D_3h) is stable. The related neutral clusters X₃Cl (X=H, Li, Na, K, and Cu) are also investigated using two methods: B3LYP, MP2, and basis set 6‐311+G(3df). For H₃Cl species, there is no a stable structure to be found. For other four X₃Cl (X=Li, Na, K, and Cu) species, there are two stable isomers, for which the bidentate structures (C_2v‐1) [see Fig. 1 (d)] are global minima. According to the general criteria for aromaticity including resonance energy (RE) and nucleus‐independent chemical shift (NICS), the five trigonal isomers exhibit a higher degree of aromaticity. Molecular orbital analysis reveals that the five trigonal X(X=H, Li, Na, K, and Cu) isomers possess only σ‐aromaticity originating from s orbitals. For the Cu ring the d orbitals do not play a significant role in the electron delocalization effects. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

五原子平面四配位碳分子[CAl_3X]和[CAl_3X]~-(X=Sn,Pb)的理论研究

在B3LYP计算水平下,使用混合基组6-311+G(d)LANl2dz首次对含Sn,Pb的五原子平面四配位碳分子[CAl3X]和[CAl3X]-(X=Sn,Pb)两个体系(分别具有17和18价电子)的异构体结构及能量进行了量子化学计算研究,并与文献已经报道的[CAl3X]和[CAl3X]-(X=Si,Ge)作了比较.结果表明,具有平面四配位碳结构的异构体int1是能量的全局最小点.相对于17e的[CAl3X]体系来说,在得到一个电子以后形成18e的[CAl3X]-int1异构体在热力学上变得更稳定.本项研究将进一步丰富五原子平面四配位碳化学.     

Structure and bonding in cyclic isomers of BAl2Hnm (n=3-6, m=-2 to +1): preference for planar tetracoordination, pyramidal tricoordination, and divalency

The structure and energetics of cyclic BAl2Hnm (n=3-6, m=-2 to +1), calculated at the B3LYP/6-311+G** and QCISD(T)/6-311++G** levels, are compared with their corresponding homocyclic boron and aluminium analogues. Structures in which the boron and aluminium atoms have coordination numbers of up to six are found to be minima. There is a parallel between structure and bonding in isomers of BAl2H(3)2- and BSi2H3. The number of structures that contain hydrogens out of the BAl2 ring plane is found to increase from BAl2H3(2-) to BAl2H6+. Double bridging at one bond is common in BAl2H5 and BAl2H6+. Similarly, species with lone pairs on the divalent boron and aluminium atoms are found to be minima on the potential energy surface of BAl2H(3)2-. BAl2H4- (2 b) is the first example of a structure with planar tetracoordinate boron and aluminium atoms in the same structure. Bridging hydrogen atoms on the B--Al bond prefer not to be in the BAl2 plane so that the pi MO is stabilised by pi-sigma mixing. This stabilisation increases with increasing number of bridging hydrogen atoms. The order of stability of the individual structures is decided by optimising the preference for lower coordination at aluminium, a higher coordination at boron and more bridging hydrogen atoms between B--Al bonds. The relative stabilisation energy (RSE) for the minimum energy structures of BAl2Hnm that contain pi-delocalisation are compared with the corresponding homocyclic aluminium and boron analogues.     

Structures and stabilities of group 13 adducts [(NHC)(EX3)] and [(NHC)2(E2X(n))] (E=B to In; X=H, Cl; n=4, 2, 0; NHC=N-heterocyclic carbene) and the search for hydrogen storage systems: a theoretical study

Quantum chemical calculations using density functional theory at the BP86/TZVPP level and ab initio calculations at the SCS-MP2/TZVPP level have been carried out for the group 13 complexes [(NHC)(EX(3))] and [(NHC)(2)(E(2)X(n))] (E=B to In; X=H, Cl; n=4, 2, 0; NHC=N-heterocyclic carbene). The monodentate Lewis acids EX(3) and the bidentate Lewis acids E(2) X(n) bind N-heterocyclic carbenes rather strongly in donor-acceptor complexes [(NHC)(EX(3))] and [(NHC)(2)(E(2)X(n))]. The equilibrium structures of the bidentate complexes depend on the electronic reference state of E(2)X(n), which may vary for different atoms E and X. All complexes [(NHC)(2)(E(2)X(4))] possess C(s) symmetry in which the NHC ligands bind in a trans conformation to the group 13 atoms E. The complexes [(NHC)(2)(E(2)H(2))] with E=B, Al, Ga have also C(s) symmetry with a trans arrangement of the NHC ligands and a planar CE(H)E(H)C moiety that has a E=E π bond. In contrast, the indium complex [(NHC)(2)(In(2) H(2))] has C(i) symmetry with pyramidal-coordinated In atoms in which the hydrogen atoms are twisted above and below the CInInC plane. The latter C(i) form is calculated for all chloride systems [(NHC)(2)(E(2)Cl(2))], but the boron complex [(NHC)(2)(B(2)Cl(2))] deviates only slightly from C(s) symmetry. The B(2) fragment in the linear coordinated complex [(NHC)(2)(B(2))] has a highly excited (3)(1)Σ(g)(-) reference state, which gives an effective B≡B triple bond with a very short interatomic distance. The heavier homologues [(NHC)(2)(E(2))] (E=Al to In) exhibit a anti-periplanar arrangement of the NHC ligands in which the E(2) fragments have a (1)(1) Δ(g) reference state and an E=E double bond. The calculated energies suggest that the dihydrogen release from the complexes [(NHC)(EH(3))] and [(NHC)(2)(E(2)H(n))] becomes energetically more favourable when atom E becomes heavier. The indium complexes should therefore be the best candidates of the investigated series for hydrogen-storage systems that could potentially deliver dihydrogen at close to ambient temperature. The hydrogenation reaction of the dimeric magnesium(I) compound [LMgMgL] (L=β-diketiminate) with [(NHC)(EH(3))] becomes increasingly exothermic with the trend B相似文献   

Redox energetics of hypercloso boron hydrides B(n)H(n) (n = 6-13) and B12X12 (X = F, Cl, OH, and CH3)

The reduction potentials (E°(Red) versus SHE) of hypercloso boron hydrides B(n)H(n) (n = 6-13) and B(12)X(12) (X = F, Cl, OH, and CH(3)) in water have been computed using the Conductor-like Polarizable Continuum Model (CPCM) and the Solvation Model Density (SMD) method for solvation modeling. The B3LYP/aug-cc-pvtz and M06-2X/aug-cc-pvtz as well as G4 level of theory were applied to determine the free energies of the first and second electron attachment (ΔG(E.A.)) to boron clusters. The solvation free energies (ΔG(solv)) greatly depend on the choice of the cavity set (UAKS, Pauling, or SMD) while the dependence on the choice of exchange/correlation functional is modest. The SMD cavity set gives the largest ΔΔG(solv) for B(n)H(n)(0/-) and B(n)H(n)(-/2-) while the UAKS cavity set gives the smallest ΔΔG(solv) value. The E°(Red) of B(n)H(n)(-/2-) (n = 6-12) with the G4/M06-2X(Pauling) (energy/solvation(cavity)) combination agrees within 0.2 V of experimental values. The experimental oxidative stability (E(1/2)) of B(n)X(n)(2-) (X = F, Cl, OH, and CH(3)) is usually located between the values predicted using the B3LYP and M06-2X functionals. The disproportionation free energies (ΔG(dpro)) of 2B(n)H(n)(-) → B(n)H(n) + B(n)H(n)(2-) reveal that the stabilities of B(n)H(n)(-) (n = 6-13) to disproportionation decrease in the order B(8)H(8)(-) > B(9)H(9)(-) > B(11)H(11)(-) > B(10)H(10)(-). The spin densities in B(12)X(12)(-) (X = F, Cl, OH, and CH(3)) tend to delocalize on the boron atoms rather than on the exterior functional groups. The partitioning of ΔG(solv)(B(n)H(n)(2-)) over spheres allows a rationalization of the nonlinear correlation between ΔG(E.A.) and E°(Red) for B(6)H(6)(-/2-), B(11)H(11)(-/2-), and B(13)H(13)(-/2-).     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

Theoretical characterizations of HAsXH (X = N, P, As, Sb, and Bi) isomers in the singlet and triplet states

The lowest singlet and triplet potential energy surfaces for all group 15 HAsXH (X = N, P, As, Sb, and Bi) systems have been explored through ab initio calculations. The geometries of the various isomers were determined at the QCISD/LANL2DZdp level and confirmed to be minima by vibrational analysis. In the case of nitrogen, the global minimum is found to be a triplet H(2)NAs structure. For the phosphorus case, singlet trans-HAs==PH is found to be global minima surrounded by large activation barriers, so that it should be observable. For arsenic, theoretical investigations demonstrate that the stability of HAsAsH isomers decreases in the order singlet trans-HAs==AsH > triplet H(2)AsAs > singlet cis-HAs==AsH > triplet HAsAsH > singlet H(2)AsAs. For antimony and bismuth, the theoretical findings suggest that the stability of HAsXH (X = Sb and Bi) systems decreases in the order triplet H(2)AsX approximately singlet trans-HAs==XH > singlet cis-HAs==XH > triplet HAsXH > triplet H(2)XAs > singlet H(2)AsX > singlet H(2)XAs. Our model calculations indicate that the relativistic effect on heavier group 15 elements should play an important role in determining the geometries as well as the stability of HAsXH molecules. The results obtained are in good agreement with the available experimental data and allow a number of predictions to be made.     

Density functional theory study of eight-atom germanium clusters: effect of electron count on cluster geometry

Density functional theory (DFT) at the hybrid B3LYP level has been applied to the germanium clusters Ge8z(z=-6, -4, -2, 0, +2, +4) using nine initial geometries. For Ge8(2-) the D2d bisdisphenoid structure predicted by the Wade-Mingos rules is not computed to be the global minimum but instead lies 3.9 kcal mol-1 above the Td tetracapped tetrahedron global minimum predicted to exhibit spherical aromaticity. The hyperelectronic clusters Ge(8)4- and Ge8(6-) have nido B8H12 and square antiprism structures, respectively, as global minima in accord with the Wade-Mingos rules and experimental data on E(8)2+(E=Sb, Bi) cations. Hypoelectronic eight-vertex clusters isoelectronic and isolobal with Ge8, Ge8(2+) and Ge(8)4+ are not known experimentally. Their computed structures include smaller polyhedra having one or more capped triangular faces as well as more open non-polyhedral structures.     

Quantum chemistry study of H+(H2O)8: a global search for its isomers by the scaled hypersphere search method, and its thermal behavior

The structures of the protonated water cluster H+(H2O)8 have been globally explored by the scaled hypersphere search method. On the Hartree-Fock potential energy surface 174 isomers were found, among which 168 were computed to be minima at the B3LYP/6-31+G** level, and their energies were further refined at the level of MP2/6-311++G(3df,2p). The global minimum on the potential energy surface computed at the B3LYP/6-31+G** level shows a cagelike structure with the "Eigen" motif, while the lowest-free-energy isomer has a five-membered-ring structure at 170 K and a chain form at 273 K. The present results are well in line with previous experimental findings. In addition, the ADMP (atom-centered density matrix propagation) simulation indicates that the extra proton in the lowest-free-energy isomer (170 K), which has a five-membered ring and the "Zundel" feature, is often in an asymmetrical hydrogen bond.