金属硼化物结构与稳定性的理论研究

用ＨＦ／３－２１Ｇａｂｉｎｉｔｉｏ法对金属硼化物ＭＢ２／ＭＢ２（Ｍ＝Ｌｉ，Ｎａ，Ｂｅ，Ｍｇ，Ａｌ）的７５个电子态结构进行能量梯度法优化，再用大基组二次且上互作用ＱＣＩＳＤ（Ｔ）／６－３１１Ｇ进行单点计算，得到了结构参数总能量，为了考察各原子簇的稳定性，还对２４个碎片的７０多个电子态，求得相应的ＱＣＩＳＤ（Ｔ）能量，在此基础上计算了原子化能、电离能、离解通道和碎片化能，得到了原子簇的稳定性规律。     

磷的四原子簇合物的从头算研究

在ＳＴＯ－３Ｇ、３－２１Ｇ、３－２１Ｇ、４－３１Ｇ和４－３１Ｇ基组水平上，对磷的四原子簇合物Ｐ４＾２＋，Ｐ４，Ｐ４＾２－，Ｐ４＾４－做了ＳＣＦ－ＭＯ从头算研究。几何优化表明，对于已知结构的族合物分子计算结果与实验测定了几何形一致，同时确定了实验上结构未知的簇合物分子的稳定２几何。根据计算结果对四原子磷簇合物系列分子的相对稳定性和几何形了讨论。     

CnAl-(n=1~11)结构的量子化学从头计算

在RHF／ROHF6－311G*水平上，优化了直链构型CnAl－（n＝1～11）的几何参数，计算了体系的能量以及它们失去C，C2，C3，Al，CA1，C2Al等六种解离方式的解离能，并与CnP－的计算结果加以对比．计算结果证实了我们以往的实验观察，即CnAl的结构稳定性具有奇偶交替的变化规律，其中n为偶数的族离子相对比较稳定．还根据对计算结果的分析，探讨了CnAl－的其它结构特性．     

直线型CnP^—（n=1～11）结构的理论研究

以激光溅射方法产生了一系列含一个磷原子的碳原子簇负离子．针对其在实验中呈现的奇偶变化规律进行了量子化学从头算研究．在HF／6－311G*水平上（对单重态为RHF，多重态为ROHF）优化了直线型CnP-（n=1～11）键长和能量，计算了相邻簇离子的能量差与成簇碳原子的平均结合能及从CnP-分别解离C、C2、C3、P、CP、C2P等6种通道所需的能量．计算发现，n为奇数的CnP-单重态（1Σ）最稳定，而，n为偶数的簇离子则以三重态（3Σ）的能量较低．所计算的CnP-各结构参数均表现出奇偶交替的变化规律，n为奇数的CnP-相对稳定．     

AlCB原子簇的从头算研究

用ab initio能量解析梯度法, 有UHF(RHF)/6-31G^*水平上优化得到AlCB的11个电子态, 并从CISD能量、振动分析、原子平均结合能以及原子簇的碎片化和碎片化能等四方面研究了AlCB的稳定性。     

铍硼四原子簇合物的量子化学从头算

铍和硼簇合物的实验和理论研究表明,它们有特殊的成键方式,BBe2,B2Be2,B6Be和B4Be等簇合物结构已由晶体粉末法测得^[1]。这类化合物可用于带电粒子的活化分析和低能电子的遏止实验,有一定应用前景^[2],对BnBem簇电子结构和成键性质进行研究,不仅对理解大簇和凝聚相的形成机制及微观结构等有重要意义,而且有一定的应用价值^[3-6],B2Be/B2Be^ 的结构^[7]已有论证。本文选用高精度大基组二次组态相互作用QCISD（T）／6－311G^**方法对BBe3,B2Be2和B3Be的电子结构进行了计算,并在HF／6－31G^*水平上作了频率计算。     

杂硼原子簇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)进行了计算和讨论. 分析了杂原子和硼原子间、相邻硼原子间的键合情况, 讨论了最稳定构型的芳香性质.     

类金属P对原子簇NinB2(n=1～6) 性质影响的理论研究

依据实验事实设计组成可调的系列原子簇 NinBP、NinB2 (n = 1～6)模型, 利用密度泛函(DFT)方法在 B3LYP/Lanl2dz 水平上对所设计构型进行优化,得到各原子簇的最稳定构型,进而研究类金属元素 P 对非晶态合金 Ni-B 性质的影响。结果表明:三元非晶态合金 Ni–B–P 的稳定性强于二元非晶态合金 Ni-B;类金属 P 起到调节金属 Ni 原子 3d 轨道布居数的作用。这些均与实验结果一致,说明模型能反映合金的性质。     

超价化合物F_nLi_(n+1)(n=1～3)一阶超极化率的理论研究

采用UMP2/6-311+G(3df)方法得到了FnLin+1(n=1～3)体系的8个全部实频的异构体构型,并研究了这些体系的非线性光学性质.结果表明,这类超价化合物具有较大的一阶超极化率值(β0=2299.2～12322.3a.u.).超价化合物中弥散的电子云使体系具有较低的跃迁能,从而决定了FnLin+1具有较大的一阶超极化率.提示非中心对称的超价化合物将是一类潜在的非线性光学分子.     

超价化合物M_(n+1)F_n(n=1,2;M=Na,K)体系的非线性光学性质的理论研究

采用UMP2/6-311+G(3df)方法研究了一系列化合物M_(n+1)F_n(n=1,2;M=Na,K)体系的几何构型及非线性光学性质(NLO).结果表明:这些体系均为超价化合物,均拥有较大的第一超极化率(β_0).尤其是Na_3F_2体系中结构2c的β_0值为29.16×10~(-49)C~3·m~3·J~(-2),是已知超价化合物Li_3F_2的2c结构的β_0值2.22×10~(-49)C~3·m~3·J~(-2)的13.1倍.     

封闭型硼烷B15H15^2—（D3h）的结构及稳定性的从头算分析

利用量子化学从头算方法在ＨＦ／６－３１Ｇ水平上，优化计算了带２个负电荷的封闭型硼烷Ｂ１５Ｈ１５＾２－（Ｄ３ｈ）的结构、能量和正则振动频率。该硼烷的简正振动频率计算表明，它是ＨＦ／６－３１Ｇ势能超曲面上的真实稳定点。理论预测它在化学上和动力学上都是稳定的。     

Theoretical Studies on Structures and Properties of Endohedral Fullerenes Complexes： XHn@C32（X=F,O, N,C; n=l-4）

Theoretical studies on structures and properties of endohedral fullerene complexes formed by encapsulating small molecules of HF, H20, NH3, and CH4 in a C32 fullerene cage, were carried out by ab initio method. Current calculations reveal that these processes to encase them in fullerene are energetically unfavorable because of the small cavity size of C32. The red shift in the F-H stretching frequency indicates the potential existence of hydrogen bonding between the HF molecule and the carbon cage.     

Theoretical Study on the Stable Structures and Aromaticities for Pnictogen Dianionic Sb4^2-, Bi4^2-, and （SbBi）2^2- Clusters

Through the theoretical calculation of structural optimization, vibrational frequencies and atomization energies with one method of density functional theory （B3LYP） and two post- Hartree-Fock approaches （MP2, CCSD（T））, several stable isomers for new three pnictogen dianionic Sb4^2-, Bi4^2-, and （SbBi）2^2- species were determined. For two homoatomic Sb4^2- and Bi4^2- species, there are three stable isomers： square （D4h）, roof-shaped （C2v-1）, and C2v-2 structure with the square isomer being the ground state. For the heteroatomic dian- ionic （SbBi）2^2- species, there are also three stable isomers： rhombus （D2h）, roof-shaped （C1）, and C2v structures with the rhombic isomer being the ground state. The calculated NICS values show that nucleus-independent chemical shifts （NICS） values of roof-shaped isomers for Sb4^2-, Bi4^2-, and （SbBi）2^2- species are all negative, consequently indicating that these roof-shaped isomers possess aromaticities. NICS values for the planar ring isomers are all positive, suggesting that these three planar ring isomers have antiaromatic characters. The aromaticity for the two stable roof-shaped and square isomers are preliminarily explained and discussed with MO analysis.     

Sn_nSm(n=1～9)团簇结构稳定性与磁学性质的理论研究

采用Saunders全局优化随机踢球模型与密度泛函理论相结合的方法,在B3LYP/SDD理论水平下研究了锡基原子团簇Sn_n(n=2~10)及锡基稀土原子钐掺杂团簇Sn_nSm(n=1~9)的几何结构、稳定性、电子性质和磁性.结果表明,团簇Sn_nSm(n=1~9)中掺杂的钐原子通常位于主团簇的表面,掺杂团簇的基态结构更倾向于具有较高对称性的三维结构;二元锡基混合团簇的平均结合能变大,稳定性增强,这主要归因于Sn—Sm键比Sn—Sn键的键能大,具有更强的相互作用;掺杂团簇具有较高的磁性,其总磁矩主要源自于钐原子4f电子的贡献;随着团簇尺寸的增加,二元团簇的总磁矩呈现出趋于饱和的现象.     

六氢吡啶和氨形成的氢键团簇C5H10NH(NH3)n(n=1～3)结构的从头算研究

在RHF/6-31G(d)水平下,对C5H10NH(NH3)n(n=1～3)氢键团簇的平衡构型进行了从头算研究,优化得到各种可能的平衡构型.C5H10NH(NH3)为线型氢键结构,而C5H10NH(NH3)2为三元环结构,C5H10NH(NH3)3为四元环结构.在MP2/6-31G(d)//B3LYP/6-31G(d)水平下,对最稳定构型C5H10NH(NH3)n(Ⅰ)(n=1～3)的分子轨道进行布居分析,并且对相应的占据轨道进行指认.C5H10NH(NH3)n(Ⅰ)(n=1～3)垂直电离势的计算结果表明,形成氢键团簇后,分子的垂直电离势降低.     

Ab Initio Theoretical Investigation on the Geometrical and Electronic Structures of AlAun-/0 (n = 2-4) Clusters

A systematic density functional theory and wave function theory investigation on the geometrical and electronic properties of AlAu n 0/-(n=2-4) clusters has been performed in this work. AlAu n-anions prove to possess ground states of the V-shaped C2v AlAu2 - , umbrella-shaped C3v AlAu3- , and perfect tetrahedral T d AlAu4- , while their neutrals favor the V-shaped C2v AlAu2 , perfect planar triangular D3h AlAu3 , and severely distorted C s AlAu4 , respectively. Aluminum aurides appear to be analogous to the corresponding aluminum hydrides, expect C s AlAu4 . Molecular orbitals (MOs) analyses also support this conclusion. Detailed orbital analyses indicate that Au 6s makes 94-96% and Au 5d makes 6-4% contribution to the Au-based orbitals in Al-Au bonds, which is smaller than the BAu n0/- series, partially reflecting the relativistic effect of gold. The one-electron detachment energies of the anions and characteristic stretching vibrational frequencies of Al-Au bonds between 100-400 cm -1 have been calculated to facilitate future experimental characterization of these clusters.     

Noble‐Gas Complexes: Theoretical Investigation of Multicenter Polynuclear Species

Ab initio calculations at the MP2 level of theory disclose the conceivable existence of neutral complexes containing four or five distinct noble gases (Ng) each bound to a distinct Be‐atom. These multicenter polynuclear Ng molecules are formally obtained by replacing the H‐atoms of CH₄ and but‐2‐yne with ? NBeNg moieties, which behave as independent monovalent ‘functional groups’. Our investigated complexes include the five homotetranuclear [C(NBeNg)₄] complexes 1 – 5 (Ng=He? Xe), the five heterotetranuclear complexes [CN₄Be₄(He)(Ne)(Ar)(Kr)] ( 6 ), [CN₄Be₄(He)(Ne)(Ar)(Xe)] ( 7 ), [CN₄Be₄(He)(Ne)(Kr)(Xe)] ( 8 ), [CN₄Be₄(He)(Ar)(Kr)(Xe)] ( 9 ), and [CN₄Be₄(Ne)(Ar)(Kr)(Xe)] ( 10 ), and the heteropentanuclear complex [HC₄N₅Be₅(He)(Ne)(Ar)(Kr)(Xe)] ( 11 ). We also investigated the five model complexes [H₃CNBeNg] (Ng=He? Xe) containing a single ? NBeNg moiety. The geometries and vibrational frequencies of all these species, invariably characterized as minimum‐energy structures, were computed at the MP2(full)/6‐31G(d,p)/SDD level of theory, and their stability with respect to the loss of the various Ng‐atoms was evaluated by single‐point calculations at the MP2(full)/6‐311G(d)/SDD level of theory. The beryllium‐Ng binding energies range from ca. 17 (Ng=He) to ca. 63 (Ng=Xe) kJ/mol, and the results of natural‐bond‐orbital (NBO) and atoms‐in‐molecules (AIM) analysis reveal that the Be? Ng interaction is essentially electrostatic for helium, neon, argon, and krypton, and has probably a small covalent contribution for xenon.