立方尖晶石氮化物电子结构的第一性原理研究.docVIP

立方尖晶石氮化物电子结构的第一性原理研究 　　摘 要：该文利用第一原理膺势平面波方法对立方尖晶石氮化物GexSi3-xN4（x=0，1，2，3）晶体结构进行自洽计算，并分析了其稳定性和电子结构，Ge3N4的能量最低，结构最稳定，其次结构比较稳定的就是GeSi2N4，Si3N4的结构式最不稳定，利用最优化后的晶格常数计算出能带和态密度，多Ge化合物（?＼-Ge3N4，?＼-Ge2SiN4）的带隙明显小于多Si（?＼-GeSi2N4，?＼-Si3N4）化合物，但是带隙却没有随Ge含量的减小而同时变大，与?＼-Ge3N4相比，?＼-Ge2SiN4表现出明显的带隙变小的行为。计算结果表明这几种物质也都具有直接带隙半导体的特征。目前对高温高压下体系稳定性存在分歧，故需要进一步去研究此体系，该文通过广义梯度近似来详细讨论氮化物这个体系就有非常大的意义和参考价值。 　　关键词：第一原理计算，电子结构，最优化 　　硅基材料是应用广泛、研究比较深入的半导体材料，由于其高热导率，高电子饱和漂移速率，和大临界击穿电压等特点，成为研制半导体微电子器件和电路的理想材料，在通信，航空，航天等领域有广阔的应用前景[1]。但是单纯的硅材料已经不能满足日益发展的半导体行业的需求，这需要研究其相关体系以发现具有更实际意义的材料。人们对此系列氮化物从理论、实验方面都进行了大量研究[2]。2001年，研究了Ge3N4用Si替代可能形成?＼-Ge2SiN4，?＼-GeSi2N4，并对他们稳定性，电子结构进行研究，认为?＼-GeSi2N4是一个不稳定的结构，而?＼-Ge2SiN4是稳定的结构[3]。此后不久通过实验来验证这个结论，在这个实验中发现了一种新的氮化物（SixGe1-x）3N4，x～0.6。对x射线图案的分析表明更小的离子Si+4离子倾向于占据八面体体心位置，而更大的离子Ge+4倾向于占据四面体位置。这个实验结果明显与理论计算得到的结果相悖。采用局域密度近似（LDA）密度泛函理论方法来研究了在Si3N4-Ge3N4体系中新生成的氮化物的稳定性，得到了与实验比较一致的结果，即最稳定的物质是?＼-GeSi2N4[4]。本文通过广义梯度近似来讨论氮化物体系有非常大的意义和参考价值。 　　一、计算过程 　　本文采用基于密度泛函理论的WIEN2k软件包研究了GexSi3-xN4（x=0，1，2，3）这4种氮化物的稳定性与电子结构。计算中的交换关联势选取Perdew Bruke Ernzerhof 广义梯度近似 （GGA-PBE） [5]。价态和芯态的分离能是-6.0 By且Rmin ??max = 7.0， where Rmin 是原子芯球半径， Kmax 平面波截止参数， Si、Ge原子芯球半径2.5， N 的为2.2。系统总能量和电荷密度在Brillouin区的积分计算使用Monkhorst-Pack方案[6]来选择k 网格点都为6??？进行自洽迭代。4种化合物的结构模型都选用立方晶系，其中以立方的Ge3N4为基础，通过不同比例的Si取代得到初始计算结构模型。然后对初始结构进行优化，直至能量变化小于10-6eV达到收敛标准并停止运算。 　　二、计算结果和讨论 　　1、结构优化结果 　　在结构优化计算中，采用广义梯度近似对晶格常数和原子位置进行了优化，得到能量最低的局域稳定状态。优化曲线如图1所示，从图中可以看出，Si的取代引入对于整体结构的稳定性有非常大的影响。而且，经过对比不难发现，随着结构中Si的取代增多，结构变得更加不稳定，即Ge3N4的能量最低，结构最稳定，其次结构比较稳定的就是GeSi2N4，Si3N4的结构式最不稳定。这和之前的研究结果是一致的，同时也是和实验中的结论也是一致的[7]。 　　2、能带结构 　　GexSi3-xN4（x=0，1，2，3）的能带结构显示 4种物质均为宽带隙直接带隙半导体，导带低和价带顶都在?点，?＼-Ge3N4，?＼-Ge2SiN4，?＼-GeSi2N4，?＼-Si3N4禁带宽度分布为2.21ev，1.95ev ，2.85ev ，3.62ev。由于GGA通常会低估电子与电子之间的交互关联作用[8]，所以计算结果很有可能比实验值偏低，这并不影响对于该类体系的定性分析。多Ge化合物（?＼-Ge3N4，?＼-Ge2SiN4）的带隙明显小于多Si（?＼-GeSi2N4，?＼-Si3N4）化合物，带隙却没有随Ge含量的减小而同时变大，带隙的大小顺序为?＼-Ge2SiN4?＼-Ge3N4?＼-GeSi2N4?＼-Si3N4。与?＼-Ge3N4相比，?＼-Ge2SiN4表现出明显的带隙变小的行为。对于这种行为的出现，需要对其态密度进行详细的分析。在?＼-Si3N4，?＼-Ge3N4，?＼-GeSi2N4的能带图中