非线性光学综述

Mo2Zr的第一性原理研究综述
一、题目研究的背景及意义
钼由于具有熔点高、强度大、硬度高、导电导热性好、热膨胀系数低，耐磨损和抗腐蚀性能强及良好的抗热震性能和耐热疲劳性能等特点，因而被广泛应用于钢铁工业、电子工业、航空航天、原子核能及金属压力加工等领域[1-3]。

然而，由于氧、氮等杂质原子在晶界的富集及其本征脆性，钼及钼合金烧结后脆性非常明显[4-5]。

而Zr具有高的熔点，高温下能与钼产生固溶而低温下溶解度很小，能对钼起到显著的固溶强化及弥散强化效果，且Zr 在加热时能大量吸收氧、氢、氮等元素，可能是改善烧结钼合金室温脆性的有效途径。

所以研究Mo-Zr 合金有着重要的实践意义，Mo-Zr的制备采用粉末冶金方法。

Mo-Zr合金主要应用于工业领域，例如燃气轮机的叶片，加热熔融玻璃时的电极，黄铜及有色金属铸造的模具，涡轮喷气飞机、火箭及核反应堆的某些组件和暴露于腐蚀化学物质的部分。

本文基于密度泛函理论对Mo2Zr的晶体结构作了优化，分别计算了Mo2Zr、Mo、bcc-Zr的晶格常数及力学性质，我们将本文计算值与他人及实验结果作了对比，发现我们计算的结果与他人及实验值吻合的比较好。

二、基本理论方法
第一性原理计算方法，或称Hartree-Fock方法，是First Principles Calculation 的译文，又称从头计算（ab initio calculation）。

第一性原理计算方法的出发点就是将多个原子构成的体系理解为由电子和原子核组成的多粒子系统，然后求解这个多粒子系统的薛定谔方程组，获得描述体系状态的波函数Φ 以及对应的本征能量。

从理论上讲，通过这两项结果就可以推导系统的所有性质。

为了求解多粒子系统的薛定谔方程，在第一性原理计算中采用了三个基本近似，即非相对论近似、Born-Oppenheimer 近似（也称为绝热近似）和Hartree-Fock 近似（单电子近似）。

常用的第一性原理计算方法主要使用的是线性缀加平面波方法（LAPW）和赝势平面波方法。

本课题所使用的是V ASP及PWSCF软件包所采用的赝势平面波方法。

赝势平面波方法基组用平面波展开，晶体势场用赝势代替。

平面波基是一种最简单的正交、完备的函数集合，它有以下优势：（1）具有良好的解析形式，是正交归一的，无需考虑重叠积分。

（2）可以尽可能多的使用平面波作为基底。

（3）平面波是非局域的，它不依赖于原子的位置。

密度泛函理论是多粒子系统基态研究的重要方法，它不但给出了将多电子问题简化为单电子问题的理论基础，也成为分子和固体的电子结构和总能量计算的有力工具。

在第一性原理计算中，要求解系统中原子间的相互作用能，就要求解如下的薛定谔方程：
第一性原理计算的核心是采用合理的近似和简化，利用量子力学求解多体问题。

组成固体的多粒子系统的薛定谔方程：
(,)(,)H H E ψ=ψr R r R （2.1）
如果不考虑其他外场的作用，晶体的哈密顿量应包括原子核和电子的动能以及这些粒子之间的相互作用能，形式上写成
N e N e H H H H -++= (2.2)
我们对研究体系进行简化，把在原子结合中起作用的价电子和内层电子分离，内层电子与原子核一起运动，构成离子实（ion core ），离子实与价电子构成凝聚态体系的基本单元。

晶体哈密顿量可以改写为：
2
222
222,112222i i i j i ij i Z Z e Z e e H m M αβααααβαααβ≠≠⎛⎫⎛⎫=-∇+-∇++- ⎪ ⎪-⎝⎭⎝⎭
∑∑∑∑∑
αr R r R (2.3) 第一项为电子动能，第二项为离子的动能，第三项和第四项是成对离子和电子之
间的静电能，第五项为电子和核之间的吸引作用。

其中，αM 表示第α个离子的质量，相应坐标是{}αR ；m 表示电子质量；i j i j =-r r r 表示电子间距离，
αβαβ=-R R R 表示原子核间距离；αZ 表示第α个原子核的电荷。

2.1 Born-Oppenheimer 绝热近似
原子核质量比电子大得多，运动速度比电子小得多，因此可以把电子的运动与原子核的运动分开考虑。

讨论电子的运动时，离子实是处在它们的瞬时位置上，可以认为离子实始终不动，电子处于固定的离子实产生的势场中，讨论离子实的运动时，不考虑电子在空间的具体分布。

这就是玻恩（M.Born ）和奥本海默（J.E. Oppenheimer ）提出的绝热近似或称为玻恩——奥本海默（Born-Oppenheimer ）近似。

对应本征能量H E 的系统波函数近似为
(,)()(,)χφψ=r R R r R （2.4）
上式表示的就是绝热近似：第一个因子()χR 是描写离子实运动的波函数，第二个因子(,)φr R 是描写多电子体系运动的波函数。

离子实满足如下方程：
()()()2
20
2E M αααχεχ⎛⎫-∇+= ⎪⎝⎭
∑R R R （2.5） 其中，R 表示所有离子实坐标{}αR 集合，()R E 0是电子体系的总能量，以平
均势的身份出现在离子实的动力学方程之中，常被称为Born-Oppenheimer 势能面。

相互作用的电子系统满足如下方程：
()()0(,)(,),BO H E φφ=r R r R R r R （2.6）
其中BO H 被称为Born-Oppenheimer 哈密顿量，表示为：
22
22
2
,1222BO
i
i i j i ij i Z Z Z e e e H m αβααβααβα
≠≠=-∇++--∑∑∑∑r R r R (2.7) 2.2 哈特利-福克（Hartree-Fock ）近似
通过绝热近似，把电子运动与离子实的运动分开，但系统的薛定谔方程仍然
是一个多体方程。

由于电子间存在的库伦相互作用，严格求解这种多电子问题是不可能的。

通过哈特利-福克（Hartree-Fock ）近似,可以将多电子的薛定谔方程简化为单电子有效势方程。

哈特利波函数将多电子波函数表述为每个独立电子波函数的连乘积形式：
()()()()12n n φφφφ=12r r r r （2.8）
哈特利-福克单电子近似方程可以表示为：
()()()()()''2'
',-HF i i i V dr E ρρφφ⎡⎤-⎢⎥∇+-=⎢⎥-⎢⎥⎣⎦
⎰r r r r r r r r （2.9） 哈特利-福克近似虽然包含了电子与电子的交换相互作用，得到了更进一步的结果，却没有考虑电子之间排斥相互作用，因此仍然具有一定的局限性，不能认为是一个严格的单电子理论。

“单电子近似”的近代理论是在密度泛函理论的基础上发展起来的。

2.3 Hohenberg-Kohn 定理
密度泛函理论是一种研究多粒子系统基态的重要方法。

建立在Hohenberg-Kohn 定理上的密度泛函理论不但建立了将多电子问题转化为单电子方程的理论基础，同时也成为分子和固体的电子结构和总能量计算的有力工具。

密度泛函理论可以使复杂的多电子波函数()φr 及其对应的薛定谔方程转化成为简单的粒子数密度函数()ρr 及其对应的计算体系。

密度泛函理论的宗旨是：原子、分子和固体的基态物理性质可以用粒子密度函数()ρr 来描述。

1964年，Hohenberg 和Kohn 利用Thomas-Fermi 模型理论提出了Hohenberg-Kohn （H-K ）第一和第二定理。

这两个定理是密度泛函理论的严格理论基础，可以归结为：
第一定理：不计自旋的全同费密子系统的基态能量是粒子数密度函数()ρr 的唯一泛函。

其核心是：粒子数密度函数()ρr 是一个决定系统基态物理性质的基本变量。

()ρr 定义为
()()()ρφφ+≡ψψr r r （2.10）
其中φ是基态波函数。

第二定理：能量泛函[]ρE 在粒子数不变的条件下对正确的粒子数密度泛函
()ρr 取极小值，并等于基态能量。

其核心是：在粒子数不变条件下能量泛函对密度函数的变分就得到系统基态的能量[]G E ρ。

获得基态，可以预测很多性质。

例如，分子的键长，振动频率，固体的晶胞边长、弹性系数张量，甚至是化学键的断裂或是生成，对电子而言都
是基态的性质。

多电子系统的哈密顿量为：
H T U V =++ （2.11）
T 为电子动能，U 为库仑排斥项，V 为由对所有粒子都相同的局域势()υr 表示的外场的影响。

对于给定的()υr ，能量泛函[]ρE 定义为：
()()()E d T U ρυρφφ≡++⎰r r r （2.12）
定义一未知的、与外场无关的泛函[]F ρ：
[]F T U ρφφ≡+ （2.13）
[][]()()[]1
2XC F T d d E ρρρρρ''=++'
-⎰⎰r r r r r r (2.14)
其中第一项是具有粒子数密度()ρr 的非相互作用电子系统的动能，第二项与无相互作用粒子模型的库仑排斥项相对应，最后一项[]XC E ρ表示交换- 关联能，代表了所有未包含在无相互作用粒子模型中的相互作用项，仍然是未知的。

目前，对于问题的求解仍然存在三方面的困难：一是如何确定粒子数密度函数()ρr ，二是如何确定动能泛函[]T ρ，这两个问题可以由Kohn 和Sham （沈吕九）在1965年提出的Kohn-Sham 方程解决。

三是如何确定交换关联能泛函[]XC E ρ，这个问题一般通过采用局域密度近似（local-density approximation, LDA ）解决。

2.4 Kohn-Sham 方程
Kohn 和Sham 构造了用无相互作用粒子模型代替有相互作用粒子哈密顿量中的相应项，而将有相互作用粒子的全部复杂性归入交换关联相互作用泛函
[]XC E ρ中去，从而转化为单电子图像。

即：
()()2
1
N
i i ρφ==∑r r (2.15)
用无相互作用粒子的动能泛函[]s T ρ代替动能泛函[]T ρ，
[]()()()21N
s i i i T d ρφφ*==-∇∑⎰r r r （2.16）
对ρ的变分可用对()i φr 的变分代替，有
()()()()110N
i i i i i E E d δρφφδφ*
=⎧⎫⎡⎤--=⎡⎤⎨⎬⎣⎦⎣⎦⎩⎭
∑⎰r r r r r （2.17） 于是可得
()()[]()()()2-+XC i i i E d E δρρνφφδρ⎧⎫'⎪⎪'∇++=⎨⎬'-⎪⎪⎩⎭
⎰r r r r r r r r （2.18）
其中
()()()()()
XC eff E V d δρρνδρ⎡⎤'⎣⎦'=++
'-⎰r r r r r r r r （2.19） （2.15），（2.18）和（2.19）一起称为Kohn-Sham 方程。

有效势()eff V r 是第一项外势()υr 、第二项库仑势和第三项交换关联势之和。

基态密度函数可从求解式（2.18）得到()i φr 后根据式(2.15)构成。

实际求解Kohn-Sham 方程的过程是一个自洽循环的过程，当输出的粒子数密度()ρr 与上一次计算结果的差值小于一定的收敛精度时，即可认为得到自洽的结果。

2.5 常用的交换关联函数
在Hohenberg-Kohn-Sham 方程的框架下，多电子薛定谔方程简化为有效的单电子Kohn-Sham 方程，这种计算方案是完全严格的，其中唯一的近似包含在交换关联能[]XC E ρ一项中，所以找出合理的交换关联能泛函形式是求解Kohn-Sham 方程的关键。

目前，局域密度近似泛函(Local Density Approximation ，LDA)和广义梯度近似泛函(Generalized Gradient Approximation ，GGA)被广泛使用。

局域密度近似泛函（LDA ）是一个简单可行而又富有实效的近似。

其基本思想是利用一个局域的均匀系统来代替非均匀系统。

也就是说，原本需要知道整个 ()ρr 函数分布才能确定空间中各点的()()XC r ερ的大小，
LDA 近似成只要给定位置r ，代入()ρr 得出该位置的ρ，就可以得到该位置的LDA
XC ε值。

即()()XC r ερ大小
只跟那个位置的电荷密度大小有关，称作局域密度近似（Local Density Approximation ），简称LDA 。

实践证明，在密度泛函理论、局域密度近似（LDFT ）框架下的计算都能得
到合理的结果。

在局域密度近似下，交换关联能[]LDA XC E ρ可以写为：
[]()()()3LDA
XC XC E d ρρερ=⎰r r r （2.20）
[]LDA
XC E ρ满足：
()()()()()()
LDA
XC XC E ρερδρδρρ⎡⎤∂⎡⎤⎣⎦⎣⎦=∂r r r r r （2.21）
其中()()XC r ερ是密度为ρ的均匀系统中每个粒子的交换-关联能。

LDA 方法普遍高估了结合能，特别是对于结合较弱的体系，误差较大。

为了弥补LDA 在计算过程当中的缺陷，不断发展出新的修正方法。

广义梯度近似（Generalized Gradient Approximation ，GGA ）是在LDA 的基础上引入了电荷密度的梯度修正，以考虑电荷分布的不均匀性，更适合处理密度的非均匀性。

与LDA 相比，GGA 在很大程度上改进了原子的交换能和相关能的计算结果，但并不总是优于LDA 。

大量计算表明，GGA 会高估晶格常数，而LDA 会低估晶格常数， LDA 和GGA 都会低估半导体带隙。

GGA 近似下，交换关联能与密度及其梯度有关：
()3,GGA XC XC E f d ρρ=∇⎰r (2.22)
通常将XC E 分为交换和相关两个部分各自寻找合适的泛函。

对于GGA ，常用的交换关联能有三种形式有：PW91，PBE ，RPBE 等形式。

2.6 布洛赫（Bloch ）定理
固体是存在大量原子核和电子的多粒子系统，需要太多的时间完成计算。

无穷大的体系是不可能精确求解的，在Born-Oppenheimer 近似和Kohn-Sham 定理基础上，需要进一步的近似才能求解多粒子体系的能级和本征函数问题。

为了求解问题，我们通常构造一个有限大小的基元，形成具有宏观周期性的物质体系。

体系中所有离子势场和其它电子的平均场被看作是周期性势场,体系中电子的运动可以简化成求解周期场作用下的单电子薛定谔方程。

由于具有空间周期性，作为反映电子状态的相空间和能量空间必然具有相同的周期性，可以在有限大小的相空间（第一布里渊区）中进行久期方程的求解，这就是周期性繁衍思想和周期性边界条件的应用价值所在。

对于周期体系当中电子行为的描述，布洛赫波是极为重要的一个概念。

布洛赫定理：周期势场中的电子波函数必定是按晶格周期函数调幅的平面波。

在周期场中，波动方程的解φ具有如下性质：
()()n
i n e φφ⋅+=k R r R r （2.23）
式(2.23)就是布洛赫定理。

它确定了周期势场中波动方程解的基本特征。

根据布洛赫定理可以把波函数写为：
()()i e u φ⋅=K r r r (2.24)
上式表达的波函数称为布洛赫函数，它是平面波与周期函数的乘积。

其中
()u r 具有与晶格同样地周期性，即
()()n u u +=r R r (2.25)
k 为一矢量，称为简约波矢，是对应于平移操作本征值的量子数。

它的物理
意义是：表示原胞之间电子波函数位相的变化。

波函数和能量本征值都和k 值有关，不同的k 值表示电子不同的状态。

晶格周期性和周期性边界条件确定了k 只能在第一布里渊区内取N （晶体原胞数目）个值。

借助于布洛赫定理，可以逐个k 求解Kohn-Sham 方程：
()()()22-2eff n n n V m φεφ⎡⎤∇+=⎢⎥⎣⎦
k k k r r r (2.26) 计算之前，首先要在布里渊选择一定数目的k 点，然后逐个k 求解本征方程。

通过运用布洛赫定理，电子无穷数量的问题就化为了在周期晶胞的第一个布里渊区内用无穷倒空间向量k 表示波函数的问题了。

三、 工作内容
1、Mo2Zr结构的搭建及优化
Mo2Zr晶体结构（C15）Mo2Zr优化后的结构Mo2Zr优化后的晶格常数为7.628 Å，与实验值7.586 Å [6]接近，他人用ASW的方法计算出结果为7.505 Å [6]，用FLG计算出的结果为7.479 Å [6]，用的方法是the augmented spherical wave (ASW) method with the von Barth-Hedin local density approximation (LDA) parametrization。

We still use the LDA (full potential LMTO called FPL in this paper，full potential LMTO with gradient correction，called FLG)。

用的方法是线性能带计算方法，采用全势LDA 加梯度修正（文献[6]中称为FLG方法）。

2、Mo2Zr的力学性质
我们利用应变-总能变化的关系式，计算出Mo2Zr的弹性常数C ij
，然后利用V oigt- Reuss-Hill (VRH)近似计算出体弾模量B和剪切模量G。

下表为目前初步计算的Mo2Zr的弹性常数及体弾模量与文献和实验值的对比结果。

剪切模量G V=159.505，G R=82.887。

3、利用优化后的结构，改变体积求对应的能量，得到Mo2Zr的E-V曲线；改变体积，求不同体积下的体弹模量，得到B-V曲线。

E
V
Mo2Zr 的E-V 曲线
b
v
Mo2Zr 的B-V 曲线
4、计算了Mo 、bcc-Zr 的晶格常数、弹性常数，为了对比，我们列出了实验值和他人的计算结构，图一为晶格常数与他人结果和实验值的对比，图二为弹性常数、体弹模量计算值与他人计算值和实验值的对比。

图一：Mo、bcc-Zr的晶格常数
图二：Mo、bcc-Zr的弹性常数、体弹模量计算值与实验值和他人结果的对比从对比结果可以看出，我们的计算值与他人及实验值吻合的很好，从而验证了本文计算的合理性。

参考文献：
[1]Moerkurova H H. 钼冶金[M]. 北京: 冶金工业出版社, 1984: 4−18.
[2]曹维成, 刘静, 任宜霞. 掺杂不同微量元素对钼材性能的影响[J]. 稀
有金属快报, 2006, 25(8): 29−31.
[3]张久兴, 刘燕琴, 刘丹敏, 等. 微量La2O3 对钼的韧化作用.中国有
色金属学报, 2004, 14(1): 13−17.
[4]Mrotzek T, Hoffmann A, Martin U. Hardening mechanisms and
recrystallization behaviour of several molybdenum
alloys[J].International Journal of Refractory Metals & Hard
Materials,2006, 24: 298−305.
[5]王德志, 刘心宇, 周美玲. MoLa2O3烧结坯间隙杂质的研究[J]. 稀有
金属材料与工程, 2000, 29(14): 266−268.
[6]H. Anton, P.C. Schmidt. Theoretical investigations of the elastic
constants in Laves phases[J].Intermetallics ,1997:449-465.
[7]Hideaki Ikehata, Naoyuki Nagasako, Tadahiko Furuta, Atsuo Fukumoto,
Kazutoshi Miwa, Takashi Saito. First-principles calculations for
development of low elastic modulus Ti alloys[J].
Phys.Rev.B,2004,70(17):174113.
[8]王保田,重金属及其氧化物高压相变的第一性原理研究[D].山西大
学,2011.。