源资培训班课程VASP软件理论基础.pdf

1、VASP量子力学计算软件包VASP软件总体介绍软件总体介绍VASP软件理论基础软件理论基础VASP软件发展史软件发展史目录目录一、一、VASP软件软件总体介绍总体介绍目前材料模拟和计算物质科学研究中最流行的商用软件之一，广泛应用于固体，表面、界面和分子等各种体系的研究原子尺度进行电子结构计算和量子力学-分子动力学模拟软件包既可以在密度泛函理论框架内求解Kohn-Sham方程，也可以在Hartree-Fock的近似下求解Roothaan方程支持格林函数方法和微扰理论使用平面波基组，电子与离子间的相互作用使用缀加投影波贋势（PAW）描述计算效率高，可实现大规模高效并行计算由维也纳大学开发，Mate

2、rials design公司代理，目前最新版本5.3VASP基本基本介绍介绍采用周期性边界条件采用周期性边界条件(或超原胞模型或超原胞模型)处理原子、分子、团簇、纳米处理原子、分子、团簇、纳米线线(或管或管)、薄膜、晶体、准晶和无定性材料，以及表面体系和固体。、薄膜、晶体、准晶和无定性材料，以及表面体系和固体。几何性质结构参数（键长、键角、晶格参数、原子位置），稳定构型等 电子性质电子态密度，能带结构，电荷密度分布，电子局域化函数（ELF）等 状态方程和力学性质体弹性模量和弹性常数等 表面性质表面重构、缺陷等结构，表面能量，表面吸附能，反应机理等 光学性质介电常数，吸收光谱，折射率等 分子动力

3、学模拟扩散系数和粘性系数 磁学性质共线和非共线性磁性、自旋轨道耦合VASP计算性质计算性质VASP材料材料表面表面催化催化凝聚凝聚态态VASP应用领域应用领域VASP发表文章发表文章Images from patents citing VASPVASP文献文献Science（48篇）31.027Nature Materials（36篇）35.749J.Am.Chem.Soc.（450篇）10.667Angew.Chem.（115篇）10.667Advanced materials（71篇）14.829Nano letters（288篇）13.025ACS Nano（246篇）12.062Phy

4、sical Review Letters（880篇）7.943Physical Review B（3946篇）3.767J.Phys.Chem.C（2170篇）4.835J.Phys.Chem.B（330篇）3.337Phys.Chem.Chem.Phys.（789篇）4.198截止到截止到2014年年VASP并行速度并行速度VASP支持支持CPU、GPU并行并行VASP并行速度并行速度VASP的并行运算的并行运算速度：速度：GPUCPUhttp:/cms.mpi.univie.ac.at/vasp-forum/forum.phpVASP技术论坛技术论坛二二、VASP软件软件理论基础理论基础根

5、据原子核和电子互相作用的原理及其基本运动规律，运用量子力学根据原子核和电子互相作用的原理及其基本运动规律，运用量子力学原理，从具体要求出发，经过一些近似处理后直接求解薛定谔方程的原理，从具体要求出发，经过一些近似处理后直接求解薛定谔方程的算法。算法。第一性原理通常是跟计算联系在一起的，是指在进行计算的时候除了第一性原理通常是跟计算联系在一起的，是指在进行计算的时候除了告诉程序你所使用的原子和他们的位置外，没有其他的实验的，经验告诉程序你所使用的原子和他们的位置外，没有其他的实验的，经验的或者半经验的参量，且具有很好的移植性。的或者半经验的参量，且具有很好的移植性。因此在材料科学领域里，有着其他

6、计算方法无法比拟的优势，其不但因此在材料科学领域里，有着其他计算方法无法比拟的优势，其不但能够使被模拟的体系性质与特征更接近真实情况，而且与真实的实验能够使被模拟的体系性质与特征更接近真实情况，而且与真实的实验相比，第一性原理计算能更快地设计出符合要求的实验，很多时候甚相比，第一性原理计算能更快地设计出符合要求的实验，很多时候甚至可以模拟实验无法实现的工作。至可以模拟实验无法实现的工作。VASP理论背景理论背景VASP理论基础理论基础薛定谔薛定谔方程方程Hartree-FockThomas-Fermipost-HFHohenberg-KohnKohn-ShamBO近似近似密度泛密度泛函理论函理

7、论CI、CC、MPx量子化学量子化学量子力学量子力学密度泛函理论密度泛函理论Schr dinger equation在量子力学中，粒子的状态用波函数描述。在量子力学中，粒子的状态用波函数描述。薛定谔方程（定态）：薛定谔方程（定态）：要准确描述一个体系，就必须通过求解薛定要准确描述一个体系，就必须通过求解薛定谔方程来得到谔方程来得到 波函数波函数n 薛定谔方程能够准确的求解出仅含一个电子薛定谔方程能够准确的求解出仅含一个电子的的H原子的波函数以及原子的波函数以及H的基态能量。的基态能量。但是，如果要处理多电子的复杂体系时，但是，如果要处理多电子的复杂体系时，求求解薛定谔方程就显得尤为复杂！解薛定

8、谔方程就显得尤为复杂！E.Schr dingerBorn-Oppenheimer Approximation玻恩玻恩-奥本海默近似奥本海默近似(BO近似近似)原子核质量远大于电子质量，一般要大3-4个数量级，因而在同样的相互作用下，根据动量守恒，原子核运动比电子运动慢得多。因此，通常在研究某一瞬间电子结构时可以忽略原子核的速度，假定它和原子核长期固定在某位置时的电子结构一样，这就是BO近似近似。其核心思想核心思想是：把核的运动和电子的运动分开处理 处理电子运动时，认为核是固定不动的 处理核运动时，认为快速运动的电子建立一个平均化了的负电荷分布，核在电子的负电荷平均场中运动量子化学量子化学：Ha

9、rtree-Fock含时薛定谔方程含时薛定谔方程定态薛定谔方程定态薛定谔方程玻恩玻恩-奥本海默近似奥本海默近似假设电子之间没有任假设电子之间没有任何相互作用，那么，何相互作用，那么，体系内多电子波函数体系内多电子波函数就可以转化为单电子就可以转化为单电子波函数的简单乘积。波函数的简单乘积。此为单电子近似。此为单电子近似。单电子哈密顿算符单电子哈密顿算符中的第三项可以看中的第三项可以看成体系中其他电子成体系中其他电子的库仑相互作用的的库仑相互作用的一个平均场近似。一个平均场近似。引入引入Slater行列式行列式量子化学量子化学：电子关联与：电子关联与post-HFHF理论采用单行列式多体波函数，

10、它假定一个电子在由原子核和其它电子形成的平均势场中独立运动，没有考虑电子间的瞬时相关作用，也就是电子关联。但是，HF理论中采用的Slater行列式波函数满足Pauli不相容原理，所以它在一定程度上包含了同自旋电子的关联作用，这部分电子关联被单独命名为交换相互作用。所以，通常所说电子关联能是指精确的基态能量与HF 能量之差。关联能一般只占体系电子总能量的0.32%，但是它可以与化学反应的反应热或者活化能相当甚至更大。一般说来，化学反应途径上电子关联程度不同。所以，在考虑化学反应和电子激发等问题时电子关联十分重要。由电子瞬时作用产生的电子关联称为动态关联。此外，还有非动态关联，它指体系几个对称性相

11、同的电子组态函数（单行列式）接近简并时导致HF方法给出很大的误差。常用的电子关联计算方法有：组态相互作用（CI）：Configuration Interaction 耦合簇方法（CC）：Coupledcluster method 多体微扰理论（MPx）：MollerPlesset perturbation theory1927年，Thomas和Fermi提出了建立在均匀电子气模型上的Thomas-Fermi模型。在均匀电子气模型中，电子不受任何外力，彼此没有相互作用。此时，描述电子运动的薛定谔方程就被简化为波动方程。Thomas-Fermi模型是一个很粗糙的模型，它是以均匀电子气的密度得到动能

12、的表达式，同时忽略了电子间的交换相关作用，直接使用比较困难。Thomas-Fermi模型存在的问题：虽然提出了可以密度为处理变量简化薛定谔方程，但是基于均匀电子气，这样子是不合理的。1964年，Hohenberg和Kohn以这个模型为基础，同时打破其能量泛函形式的束缚，创立了严格的密度泛函理论。Thomas-Fermi Model汤姆斯汤姆斯-费米模型费米模型(TF)1964年，现代密度泛函理论解开序幕：HK 第一定理证明了体系的基态能量只与电荷密度有关，是电荷密度的泛函。HK 第二定理证明了以基态密度为变量，将体系能量最小化之后就得到了基态能量。Hohenberg-Kohn(HK)Theor

13、em=|eeexteeextET VVTVV 计算的复杂转移到了如何精确计算的复杂转移到了如何精确交换相关泛函交换相关泛函Exc的表达形式的表达形式动能项动能项原子核与电原子核与电子之间的相子之间的相互作用势互作用势哈特里哈特里-库仑能库仑能交换相关能交换相关能Kohn-Sham 定理（KS）：实际体系动能的主要部分可以用无相互作用体系的动能来描述。其中动能的误差部分和相互作用能与库伦作用能之差归入交换相关泛函。若交换相关泛函严格精确，则无相互作用体系的基态密度将严格地等价于真实体系的基态密度，也将得到精确的基态能量。事实上，寻找越来越精确和普适的交换相关泛函代表了现代密度泛函方法发展的主线。

14、Kohn-Sham Theorem交换相关泛函：交换相关泛函：LDA：Local Density Approximations,LDA,局域密度近似将均匀电子气的交换关联能密度代替非均匀电子气的交换关联密度。GGA：Generalized Gradient Approximations,GGA,广义梯度近似考虑了电子密度的非均匀性，其中GGA-PBE应用最广。还有经过修正的rPBE、PBEsol等。meta-GGA：此类泛函比GGA多的变量是动能密度(kinetic energy density)或是局域自旋密度的二阶导数。主要包括：TPSS、revTPSS、M06-L等。hybrid Fun

15、ctional：综合了Hartree-Fock和DFT。对于固体的电子结构及性质计算更加准确。常用的杂化密度泛函包括：B3LYP、PBE0、HSE06等。vdW-DFT：考虑了范德华相互作用。选用标准：选用标准：1.具体问题具体分析，多看文献很重要。具体问题具体分析，多看文献很重要。2.meta-GGA和杂化泛函在很多方面更精确，但是也同样计算量更巨大。和杂化泛函在很多方面更精确，但是也同样计算量更巨大。Density Functional Theory(DFT)LDA局域（自旋）密度近似(L(S)DA)是对于寻找密度泛函理论中交换关联泛函的的第一台阶。Kohn和Sham提出把具有相同密度的均

16、匀电子的交换相关作用泛函近似的作为非均匀体系的交换相关泛函，得到自平均电子气模型（uniform electron gas,UEG）。这样，对于实际体系电子的近似交换相关泛函可表示为：(r)dr (r)drLDALDALDAxcxcxcxcEEEE 在体系空间中的每一点电子密度决定了该点的交换能密度和相关能密度，与其它的电子密度无关。将Kohn-Sham方程作LDA近似后，使得DFT可以应用到大部分的金属和半导体体系的基态性质，比如结合能，晶格常数，晶体力学性质；对原子，分子的键角、键长、电子的自旋密度、声子的振动频率也能做出较为准确的计算。但是，由于LDA是对于多电子体系作了均匀电子气的近似处理，因此LDA对于自旋非极化的系统给出能量的全局最小值，LDA能够较为准确的处理，而对于磁性材料,电子能量会有多个局部最小值，因此LDA精度并不理想。同时，由局域密度近似计算得到的分子和固体的结合能偏高，金属d带宽度及半导体的禁带宽度的计算结果与实验相比总是偏小。因此，在之后的发展中，又对LDA作了进一步的修正，发展了广义梯度近似、在位库伦修正、自相互作用修正等方法。GGA由于现代密度泛函理论的