Fe3C化合物电子结构磁性和硬度的第一性原理研究.pdfVIP

Fe C 化合物电子结构、磁性和硬度的第一性原理研究1 3 闫牧夫，吴业琼 哈尔滨工业大学材料科学与工程学院，哈尔滨（150001 ） E-mail ：yanmufu@ 摘 要：本文应用第一性原理研究了正交和六方结构Fe C 化合物的电子结构、磁性和硬度。 3 计算结果表明，正交结构的Fe C 化合物较六方结构的Fe C 化合物易形成且稳定；两种结构 3 3 中C 和Fe 原子都成共价键，正交结构的共价作用强于六方结构；两种结构中C 原子具有相 同的磁矩，但六方结构中的Fe 原子磁矩大于正交结构中的Fe 原子磁矩；文中两种硬度计算 方法获得的正交结构Fe3C 化合物的硬度值都高于六方结构的硬度值，且前一种计算方法获 得的两种结构的硬度值都大于后一种计算方法得到的硬度值。 关键词：Fe3C；电子结构；磁性；硬度；第一性原理 0. 引言 渗碳体(Fe C)是钢中常见且非常重要的铁碳化合物，特别是在钢和铁的渗碳过程中往往 3 伴随着渗碳体的析出，其数量、形态、尺寸和分布影响钢和渗层的性能。渗碳体是一种具有 复杂晶体结构的间隙化合物，具有正交和六方两种晶体结构，硬度较高，塑韧性较差。对于 [1~3] 正交结构Fe C 的弹性性能和电子结构的研究有见报道 ，而对于六方结构Fe C 的研究为 3 3 之甚少。纯Fe3C 化合物硬度的定量化研究，目前还是一片空白。因此，本文应用第一性原 理计算的方法研究了两种结构Fe3C 化合物的能量、电子结构及磁性和硬度等性能，对这两 种结构的稳定性和性能等进行了比较。 1. 计算方法 第一性原理计算采用剑桥大学研发的CASTEP (Cambridge Serial Total Energy Package) [4] [5] [6] 软件包 ，交换关联势选用广义梯度近似(GGA) ，离子势选用超软(ultra soft)势 。为了获 得计算精度与收敛速度的平衡，平面波截断能量Eout =300eV ，能量收敛精度为1 × 10-5eV， 作用在每个原子上的力不大于0.03eV/Å，内应力不大于0.05GPa，原子位移不大于0.001Å 。 布里渊区积分采用Monkhorst-pack 特殊k 点[7]对第一布里渊区求和，正交和六方结构的k 点 网格划分分别为5×4×6 和6×6×6 。为研究Fe3C 的磁性性能，采用了自选极化近似。 2. 计算结果与讨论 正交结构的Fe C 晶胞中每个碳原子被6 个铁原子所包围，构成八面体结构，每个铁原 3 子又为两个八面体共用，即碳原子配位数为6，铁原子配位数为2 。其空间群为PNMA(62) ， 晶格常数为a=5.09Å，b=6.748Å ，c=4.523Å，α=β=γ=90° ，处于8d 和4c 位置的铁原子分别 命名为Fe1 和Fe2 ，其晶胞结构示意图如图 1(a)所示。六方结构的Fe3C 晶胞中碳原子也处 在6 个铁原子所构成的八面体间隙中，碳原子和铁原子的配位数分别为6 和2 。其空间群为 P6322(182) ，晶格常数为a=b=4.767Å，c=4.354Å，α=β=90 ，γ=120° ，其晶胞结构示意图如图 1(b)所示。 1本课题得到高等学校博士学科点专项科研基金项目(20060213017) 资助。 - 1 -