第一性原理计算方法在新材料研究中的应用
Pd基Heusler合金Pd2CrGa晶体结构预测 采用基于基因遗传算法软件包USPEX对Pd2CrGa进行晶体结构预测。 步骤:首先建立Pd2CrGa的L21结构模型,并采用自选极化方式对Pd2CrGa晶胞进行结构优化,体积优化,通过对一系列不同体积和能量的拟合,确定基态能量,进而获得结构的晶格常数。 通过基于密度泛函理论(DFT)的VASP软件包,采用经相对论校正的投影缀加波PAW(Projector Augmented Wave)方法, 交换关联能采用广义梯度近似(GGA), 平面波截断能选取500eV, K点网格采用12 ×12 × 12, 计算过程均采用自旋极化的处理方式对Pd2CrGa晶体结构的磁性、态密度、弹性常数进行了计算。 Pd基Heusler合金Pd2CrGa晶体结构预测 Pd2MGa(M=Cr, Fe)的结构 (a)L21结构;(b)四方结构( -Pd原子 -M原子 - Ga原子) Pd基Heusler合金Pd2CrGa晶体结构预测 * 磁性 Pd2CrGa呈现铁磁性,Cr原子是Pd2CrGa总磁矩的主要来源。 Pd基Heusler合金Pd2CrGa晶体结构预测 从图中可以看出,Pd2CrGa在两种状态下均表现处很强的自旋极化,而且在费米面附近尤为明显Cr. 原子的态密度差异是Pd2CrGa总态密度差异的主要来源,而Pd原子、Ga原子的自旋向上和自旋向下态密度的对称性较高,对总磁矩的贡献有限. 所以,Cr原子是Pd2CrGa磁性的主要贡献者 态密度 Pd基Heusler合金Pd2CrGa晶体结构预测 弹性常数 对于Pd2CrGa四方结构包含六个独立的弹性常数C11,C12,C13,C33,C44和C66. 四方结构的体积模量可由公式: B=1/9(2 C11 +2 C12 +4 C13 + C33)得到。 结果 C11 C12 C13 C33 C44 C66 B Pd2CrGa 218.07 70.89 124.36 184.99 23.74 75.53 140.04 四方结构的弹性常数满足四方相稳定判据: Pd2MnGa的四方结构满足稳定性判据.能稳定存在。 第一性原理计算实例三 FeZn13弹性常数的第一原理计算 自20世纪90年代起,热镀锌在汽车上、大型建筑上的应用迅速发展,新应用对镀锌板质量提出了更高的要求。热镀锌的几个相层中,含铁量为5%—6%的ζ相[1](FeZn13相)是热镀锌过程中一个非常重要的相层。提高镀锌钢材的力学性能,改善锌镀层粘附性能,都需要深入细致的研究FeZn13相。 应用第一原理方法计算材料的电荷分布和弹性常数,可以获得材料的特征参数,从而能够表征、预测甚至设计材料的结构与性能,为以后对FeZn13的验研究和理论研究提供理论依据。 FeZn13弹性常数的第一原理计算 材料弹性是人们选择和使用材料的依据之一;作为减振元件或结构则要求材料应变能要高。近代航空、航天、无线电及其精密仪器仪表工业对材料的弹性有更高的要求,不仅要有高的弹性模量,而且还要恒定,同时,准确测定材料的弹性常数,对于研究材料间原子的相互作用和相变都具有工程和理论意义。 FeZn13弹性常数的第一原理计算 FeZn13晶胞模型 FeZn13晶胞的拓展图 FeZn13 晶格常数a=13.394?,b=7.598?,c=5.066?,α=γ=90o, β=127.23o 单斜立方结构,(空间群12:C2/m),晶胞中一共有28个原子,其中Z方向的四条棱的中间和单斜立方的中心各有一个Fe原子,晶胞中和面心位置一共有26个Zn原子。 FeZn13弹性常数的第一原理计算 Fe原子和Zn原子连接的4种键长 x y z Fe1 0.0000 0.0000 ,0.5000 Zn2 0.0000 0.0000 0.0000 Zn3 0.2227 0.0000 0.0693 Zn4 0.1144 0.5000 0.2929 Zn5 0.1764 0.1792 0.5517 Zn6 0.0753 0.2898 0.8360 a=1.32313nm, b=0.76326nm, c=0.50567nm β=126.966o, Vol=0.2041nm3 晶胞中个原子的位置 FeZn13弹性常数的第一原理计算 弹性常数是表征材料弹性的量。对大多数晶体来说,弹性常数有21个非零独立分量。而高对称性的晶体弹性常数还可以进一步减少。对于单斜晶体来说,一共有13个非零独立分量,要想得到这些弹性常数,只需对未变形的晶胞结构施加微小的变形,得到变形后晶胞的受力情况。 弹性常数的第一性原理计算可以分为两个阶段:(1)对晶体施加弹性极限范围内的应力。在计算中以矩阵的形