密度泛函理论.doc

摘要：介绍了密度泛函理论的发展与完善，运用密度泛函理论研究了钒(Vanadium)在高压下的结构相变。通过计算体心立方结构的钒在不同压强下剪切弹性系数C44,发现当压强约95GPa时C440,说明体心立方结构的钒在此条件下是不稳定的。进一步计算分析得到钒在高压下发生了从体心立方到菱面体的结构相变,相变压强约70GPa,这一结果与实验结果符合。还首次发现当压强约380GPa时,将会发生菱面体到体心立方的结构相变,这有待实验的验证。

引言：相变的研究受到广泛重视,通过相变研究可以认识物质的内部结构,可以了解原子核的内部性质。尤其是极端条件下—高温、高压下相变的研究一直是人们关注的热点,能量很高的重离子反应能形成高温、高密的区域,在这种条件下会出现许多奇异现象[1]。原子在高压下也会出现许多新的特征,如发生结构相变。过渡金属钒由于有较高的超导转变温度Tc,最近成为实验和理论研究的主题[2—8]。Ishizuka等[2]对钒的实验研究发现:常压下钒的转变温度Tc为5.3K,并随压强成线性增长的关系,当压强为120GPa时Tc=17.2K(迄今是金属中最大的Tc),但压强大于120GPa,Tc出现了反常,即不再随压强成线性增长而保持不变。Takemura等[8]对高压下的钒进行了X射线衍射实验,结果显示状态方程并没有奇异性,体心立方结构的钒在压强达到154GPa时仍是稳定的。Suzuki和Ostani利用第一性原理对进行了计算,发现横向声子模在加压下有明显的软化,当压强约130GPa时变成虚的,能说明可能发生了结构相变,但并未给出相变细节[3]。Nirmal等[4]理论计算表明,压强约140GPa时会发生体心立方到简立方(sc)的结构相变。Landa等[5,6]计算了体心立方结构的钒在加压下剪切弹性系数C44的大小,发现压强约200GPa时会出现力学不稳定,并用费米面嵌套解释了不稳定的原因,但并没有给出相变后的结构。最近Ding等[7]在常温下首次从实验上得到当准静压约63GPa时钒会发生从体心立方到菱面体的结构相变,并分析了产生结构相变的原因。他们认为,排除传统的s-d电子跃迁的驱动,相变可能与来自于费米面嵌套、带的Jahn-Teller扭曲以及电子拓扑跃迁等因素有关。

基于如上原因,本文运用密度泛函理论研究钒在高压下的结构相变,即通过计算体心立方结构的钒在不同压强下的剪切弹性系数C44得到结构相变的信息。研究得到钒在高压下发生从体心立方到菱面体的结构相变,相变压强为70GPa;当压强约380GPa时,又会发生菱面体到体心立方的结构相变,这需要实验的论证。

密度泛函理论是凝聚态物理和计算化学领域最常用的方法之一,是一种研究多电子体系电子结构的量子力学方法,理论用只有三个变量的电子密度函数代替有多个变量的电子波函数,使处理问题更加方便[34]。密度泛函理论的概念以Thomas-Fermi模型为出发点,在Hohenberg-Kohn定理提出后得到了发展[36,37]。

绝热近似

通过建立多电子粒子系统的定态薛定愕方程,计算固体的电子能级。价电子对物质的结构和性质上起重要作用,将内层电子和原子核看成一个整体叫做离子实,多电子粒子系统由价电子和离子实组成,求解多电子粒子系统的定态薛定谔得方程,即求解离子实和价电子组成的系统的定态薛定愕方程:

电子坐标用表示;离子实坐标用表示。若不考虑其他外场的作用,哈密顿量表示如下:

其中和分别表示电子体系和离子实的哈密顿量,为电子与离子实的相互作用能。电子体系哈密顿量为:

其中表示电子动能项;表示库企相互作用势能项。对除外的所有电子进行求和。离子实的哈密顿量为:

第一项表示离子实的动能;第二项表示离子实间的相互作用能。对除外的所有的离子实求和,表示第个离子实的质量。这里假定离子实间的作用能与离子实之间的位矢差有关。电子与离子实的相互作用能形式上为:

式(2.1.1)~(2.1.5)中每立方米对/和/求和的数量级为1029,所以需要进一步的近似。

在电子体系的哈密顿量中,只出现电子坐标;在离子实的哈密顿量中,只出现离子实坐标;电子坐标和离子实坐标及同时出现在离子实与电子相互作用项中。由于他们的作用

得到对于单粒子波函数的单电子方程。以上三式统称Kohn-Sham方程。根据Hohenberg-Kohn定理,将该系统的波函数、基态的能量以及各物理算符的期待值等由求出的粒子数密度函数准确地得出,单粒子能量为得到的本征值。

交换关联近似

应用Kohn-Sham方程,能用有效的单电子问题描述有效的多电子系统基态特性问题,但是由于在交换关联能泛函仍然包含着多粒子系统相互作用的所有复杂性,因此解K