化学学报 ›› 2008, Vol. 66 ›› Issue (8): 914-922. 上一篇    下一篇

研究论文

单质Si晶格稳定性的第一原理研究

陶辉锦*,1,2,尹健3,尹志民1,张传福2,谢佑卿1,3   

  1. (1中南大学材料科学与工程学院 长沙 410083)
    (2中南大学冶金工程博士后流动站 长沙 410083)
    (3中南大学粉末冶金国家重点实验室 长沙 410083)
  • 投稿日期:2007-09-19 修回日期:2007-12-06 发布日期:2008-04-28
  • 通讯作者: 陶辉锦

First Principles Study of the Lattice Stabilities of Elemental Si

TAO Hui-Jin*,1,2 YIN Jian3 YIN Zhi-Min1 ZHANG Chuan-Fu2 XIE You-Qing1,3   

  1. (1 School of Materials Science and Engineering, Central South University, Changsha 410083)
    (2 Postdoctoral Mobile Station of Metallurgy Engineering, Central South University, Changsha 410083)
    (3 State Key Laboratory of Powder Metallurgy, Central South University, Changsha 410083)
  • Received:2007-09-19 Revised:2007-12-06 Published:2008-04-28
  • Contact: TAO Hui-Jin

采用SGTE纯单质数据库中Gibbs能的表达式, 结合JANAF热力学数据, 采用最小二乘法对SGTE纯单质数据库中元素Si的Gibbs能表达式进行了重新评估, 得到了比SGTE数据库更精确的结果; 同时, 将SGTE数据库中CALPHAD方法得到的晶格稳定参数外推至0 K, 与第一原理总能赝势平面波和投影缀加波方法的结果进行了对比, 发现第一原理总能赝势平面波方法得到的晶格稳定参数结果为 > > >0, 与CALPHAD方法外推结果一致. 同时, 研究发现: 第一原理总能平面波赝势方法计算的diamond-Si的晶格常数和原子体积比实验值大, fcc-, hcp-和bcc-Si的结果比第一原理投影缀加波方法的小, 总能绝对值出现了类似的结果; 四种结构中所有的价态电子密度分布到0至-15 eV区间, 部分s态电子转化为p态电子, 其中diamond-Si的转化数目最多, 结构最稳定, 这与电子态密度的计算结果一致.

关键词: Si, Gibbs能, 最小二乘法, 晶格稳定性, 第一原理

Taking the form of Gibbs energy functions in the SGTE (Scientific Group Thermodata Europe) database of pure elements, adopting the newly available thermochemical reference data of JANAF (Joint Army-Navy-Air Force) database, the correction of transition data and reassessment of the parameters of Gibbs energy of elemental Si have been performed with the least-square method. The results agree more accurately with JANAF data than with the SGTE data. At the same time, the lattice stability parameters obtained by the CALPHAD method in the SGTE database have been extrapolated to 0 K, and these results have been compared with those of the total energy plane wave pseudopotential method and projector augmented wave method in first principles.It was found that the result of the total energy plane wave pseudopotential method agrees with that of CALPHAD, > > >0. Besides, it was found that the results of lattice constants and atomic volumes of diamond-Si calculated by the total energy plane wave pseudopotential method were much larger than experimental data and those of fcc-, hcp- and bcc-Si were smaller than those by the projector augmented wave method in first principles. The same case occurred in total energy calculations. All the electrons in valence distribute in the energy zone from -15 to 0 eV, part of s state electrons per atom in diamond-, fcc-, hcp-, and bcc-Si are changed into p state electrons, and diamond-Si has the largest number of transition and is the most stable phase, which were further proved by the results of density of state calculations.

Key words: Si, Gibbs energy, least-square method, lattice stability, first principles