化学学报 ›› 2016, Vol. 74 ›› Issue (9): 758-763.DOI: 10.6023/A16050230 上一篇    下一篇

研究论文

[BxAl13-x]-(x=0~13)二元团簇的密度泛函理论研究

黄敏, 徐畅, 程龙玖   

  1. 安徽大学化学化工学院 合肥 260301
  • 投稿日期:2016-05-11 发布日期:2016-08-10
  • 通讯作者: 程龙玖 E-mail:clj@ustc.edu
  • 基金资助:

    项目受国家自然科学基金(Nos.21273008,21573001)资助.

Density Functional Theory Studies of the Binary Systems[BxAl13-x]- (x=0~13)

Huang Min, Xu Chang, Cheng Longjiu   

  1. School of Chemistry and Chemical Engineering, Anhui University, Hefei 260301
  • Received:2016-05-11 Published:2016-08-10
  • Supported by:

    Project supported by the National Natural Science Foundation of China (Nos. 21273008, 21573001).

采用遗传算法(Genetic algorithm,GA)对B-Al二元团簇[BxAl13-x]-x=0~13)进行了全局结构搜索,并在密度泛函理论(Density Functional Theory,DFT)下对其进行了优化计算. 结果表明在二元团簇中,随着硼原子的数量逐渐增加,团簇的结构由二十面体(3D)逐渐转化为准平面(2D)且团簇结构由B与Al原子的比例所决定的. 当x=0~7时,Al原子占多数,团簇保持3D结构;反之团簇呈现2D构型;在x=7~8时,团簇的结构发生由3D向2D的转变. 能量决定着团簇的稳定性(例如相对能量). 能量越小,其结构越稳定. 在计算相对能量Erel,当x=1时,其相对能量最低,团簇结构最稳定. 为了进一步了解团簇的稳定性,计算了[BxAl13-x]- 团簇的HOMO-LUMO能级差(EH-L)和垂直电离能(Vertical Detachment Energies,VDE),并且发现它们的值在整体上是随着B原子数量的增加而减小,表明其团簇的稳定性逐渐减弱. 在所有的团簇中,BAl12-的能级差最大,结构最稳定. 因此文章中对其进行了分子轨道分析,发现当一个B原子替代了Al13-团簇中的中心Al原子时,所得到的BAl12-的电子壳层结构的1s2和1p6几乎和Al13-团簇保持一致. 此外,对于x=13时的准平面全硼团簇B13-,文章中用适应性自然密度划分(Adaptive Natural Density Partitioning,AdNDP)对其化学成键进行了分析,结果显示B13-有8个π电子,具有π反芳香性.

关键词: 密度泛函理论, [BxAl13-x]- (x=0~13)团簇, 优化, 分子轨道, 适应性自然密度划分

In this paper, the global minimum search and structural optimization for the B-Al binary clusters [BxAl13-x]- (x=0~13) are performed using the genetic algorithm (GA) method coupled with density functional theory (DFT). The effects of composition on the atomic structures, electronic properties including the energy gaps and vertical detachment energies of B-Al binary clusters are discussed. The results distinctly reveal a three dimensional (3D) to (quasi-)planar (2D) structural transition as a function of x upon increasing the number of boron atoms. When x is in the range of 0 to 7, the clusters are Al-rich and the B-Al binary systems maintain the 3D structure. Whereas, the binary system trends to be quasi-planar structure, and the critical B:Al ratios for the 2D-3D transition are between x=7 and 8. To study the stability of the [BxAl13-x]- clusters, we defined the relative energy (Erel=E([BxAl13-x]-)-xE(B13-)/13–yE(Al13-)/13), where the cluster with a more negative Erel is more stable. At x=1, Erel is the most negative, indicating the highest stability. In order to further understand the stability of clusters, the vertical detachment energies (VDE) and the HOMO-LUMO energy gaps (EH-L) of [BxAl13-x]- (x=0~13) clusters are also calculated. The results show that the energy decreases with the increasing number of B atoms, indicating a lower stability. The largest EH-L of BAl12- cluster indicates that it is the most stable among all the series of this clusters. Molecular orbitals (MO) of BAl12- cluster are analyzed and the result shows that the electronic shells of 1s2 and 1p6 are virtually unchanged when the central Al atom is replaced by the B atom. It also indicates that the electron shell closing model could be regarded as a simple but valid tool for explaining the structures and stabilities of metal clusters. Chemical bonding analysis by Adaptive Natural Density Partitioning (AdNDP) method for the B13- cluster reveals that it is a π-antiaromatic system with 8 delocalized π-electrons.

Key words: density functional theory, [BxAl13-x]- (x=0~13) clusters, optimization, molecular orbitals, adaptive natural density partitioning