DyCu电子结构的第一性原理计算.pdfVIP

1 DyCu 电子结构的第一性原理计算 史耀君，杜宇雷，陈国良，陈光 南京理工大学材料科学与工程系，南京 （210094 ） E-mail ：shiyaojun@163.com 摘 要：本文利用基于密度泛函理论（DFT ）的全势线性缀加平面波法（FP_LAPW ）计算 了DyCu 的电子结构。利用体积优化比较了广义梯度近似法（GGA ）和LDA+U 法的适用性， 结果证明了LDA+U 法的正确性。分析比较了电子自旋向上和自旋向下的电子态密度和能带 结构的异同，结果表明DyCu 的电子相互作用的共价键较弱，具有典型的金属键特征，预测 其在室温下具有明显的室温塑性。 关键词： 第一性原理；DyCu；电子结构；LDA+U 法 中图分类号： TG146.4 1. 引言 基于第一性原理的局域自旋密度近似（LSDA ）或广义梯度近似（GGA ）法对于弱关联 系统（一般为主族元素）的计算是非常精确的，取得了巨大的成功，但是对于强关联系统（含 有局域化d 电子或f 电子的过渡族或稀土族元素）的计算，往往导致错误而荒谬的结果。具 有5f 电子系统的UO2 用LSDA 计算的电子态密度（DOS ）的结果是5f 态电子被牵制在费米 能级（EF ）处，从而得出的结果是UO2 是金属性的导体，这显然是错误的，而使用LDA+U 法后得到的DOS 图中，在费米能级处是明显的能隙，UO2 是绝缘体，这才是正确的计算结 果[1] 。在强关联系统中，由于非常局域化的d 态电子或f 态电子之间的库仑排斥力较强，甚 至比自身的动能都突出，因而使它们的运动相互关联。LSDA 和GGA 法是建立在单电子模 型的基础上，而单电子模型对强关联系统的局域化电子是不适用的，所以会导致错误的结果。 LDA+U 法就是在LSDA 法的基础上对其进行了修正，把强关联系统的电子分成两个子系统， 局域化的d 态或f 态电子的能量用Hubbard 项U 值描述，其它巡回电子仍用LSDA 描述[2]， Hubbard 项 U 值描述的是在同一地点两个电子之间的库仑力作用能。对于本实验计算的 DyCu 的电子结构，文献中尚未见报道，但根据上述分析，Dy 原子含有4f 电子，属于强关 联系统，简单用第一性原理的LSDA 或GGA 法肯定会导致错误的结果，所以采用LDA+U 法可以解决DyCu 电子结构的计算难题。 2. 计算方法 [3] 本实验计算利用最新的 WIEN2k 程序包中的全势线形缀加平面波法（FP_LAPW ） 。 DyCu 是简单立方格子（Pm-3m ），具有B2 （CsCl 型）结构，属于空间群221[4, 5] 。晶格参 数由X-ray 粉末衍射确定，其值为a = b = c = 3.460 Å ，α = β = γ = 90° 。图1 是DyCu 晶体 0 0 0 结构的图例，从图中可以知道Dy 原子占据着立方单胞八个顶角的位置（0, 0, 0 ），Cu 原子 占据着立方单胞的中心位置（0.5, 0.5, 0.5 ）。 在计算中，首先对DyCu 晶体进行了体积优化，比较了广义梯度近似法（GGA ）和LDA+U 法，在确定了 LDA+U 法更适合本体系后，计算了 DyCu 的电子结构，交换关联势取 [6] Perdew-Burke-Ernzerhof （PBE ）形式 。Dy 和Cu 原子的球半径分别选取为2.7a.u.和2.2a.u. ， 价态与束缚态的分界能量选取－6.0 Ryd （－81.6 eV），自洽循环计算中总能量的收敛标准 是10 µRy/f.u.。在LDA+U 的计算过程中，要考虑自旋极化（spin polarized ）和 局域轨道势 1 本课题得到教育部高等学校博士学科点专项科研基金(项目编号：20040288018)资助