侯柱峰vasp教程

VASP

Cij

§1§2§3

σ

[1,2]Voigt

(zfhou@)

2006

8

21

VASP

Cij

§1

[1,2]

Cij

xx→1,yy→2,zz→3,yz→4,xz→5

e1

2e5

=

1

12e41

2e4

.e3 σi=

1

,

i

=0

Cij=

1

,

i j

=0

1

123

xy→6

(1)

(2)

(3)

2

i,j=1

6

Cijeiej+...

→′

R

(4)

E(V0,0) →R

V0

→′ →

R=R (I+ )

(5)

I

a

=e=(e1,e2,e3,e4,e5,e6)

(△E=E(V, ) E(V,0))

C11,C12,C13,C33

C44

§2

√ →a1

2

a

1

R= 3

2a00

0c

c

a2a3

= .

e=(δ,δ,0,0,0,0)

C11+C12[3]

E

V=

1

V=

1

V=C0

44δ2

e=(δ,δ,δ,0,0,0)

C11C12C13

C33

E

(6)

(10)

AlN2

3.11553

1.0000000.0000000.000000-0.5000000.8660250.0000000.0000000.0000001.60500022Direct

0.000000000.000000000.000000000.3333333

VASP计算的理论及实践总结

一、赝势的选取 二、收敛测试

1、VASP测试截断能和K 点 2、MS测试

三、结构弛豫

四、VASP的使用流程（计算性质）

1、VASP的四个输入文件的设置 2、输出文件的查看及指令 3、计算单电能

(1) 测试截断能 (2) 测试K点 4、进行结构优化 5、计算弹性常数 6、一些常用指令

一、赝势的选取

VASP赝势库中分为：PP和PAW两种势，PP又分为SP（标准）和USPP（超软）。

交换关联函数分为：LDA（局域密度近似）和GGA（广义梯度近似）。GGA又分为PW91和PBE。

在VASP中，其中pot ,pot-gga是属于超软势（使用较少）。Paw, paw-pbe ,和paw-gga是属于PAW。采用较多的是PAW-pbe 和PAW-gga。

此外vasp 中的赝势分为几种，包扩标准赝势(没有下标的)、还有硬（harder）赝势(_h)、软（softer）赝势(_s), 所谓的硬（难以赝化），就是指该元素原子的截断动能比较大，假想的势能与实际比较接近，计算得到的结果准确，但比较耗时，难以收敛。软（容易赝化），表示该元素原子的截断动能比较小，赝势模型比较粗糙，但相对简单，可以使计算很快收

打包压缩命令：tar zcvf 文件名.tar.gz 源文件名

采用link方式避免重复的文件浪费内存：ln –s 源文件 命名

Eg：ln –s ../optic/MME ./EuO.mme

交互式绘图工具gnuplot: 命令行打:gnuplot进入 格式:plot “文件名” 退出:quit

求磁矩:getmag

如果由于节点掉线在提交任务后秒退,不输出outcar可以指定节点提交任务 LJRS -l nodes=c0104:ppn=4

chmod +x 名称——使脚本可以执行 构造potcar，以A和B元素为例：

如果是以Z为拓展名的文件： zcat A/POTCAR.Z B/POTCAR.Z > POTCAR 如果是解压后的potcar文件： cat A/POTCAR B/POTCAR > POTCAR

INCAR中的RWIGS通过POTCAR文件获得

单位晶胞体积：grep “vol” OUTCAR

自动计算加应力情况下的最优化情况：

C/a ：vaspcaopt softmode-e -=*(画曲线）

如果要看某种材料是FM还是AFM

实例1

VASP赝势文件POTCAR详解 VASP赝势文件POTCAR详解

POTCAR 赝势文件

可以理解为分子力学模拟中的力场文件 但包括的信息更多 VASP4.6将各元素优化的INCAR里的参数也包括在这里了，作为支持PREC的缺省选择 通常各元素的POTCAR已经包括在软件包里了

我们只需要按照POSCAR里的顺序，将各元素的POTCAR按顺序连接起来就可以了 如以下命令：

cat file1 file2 file3 > POTCAR

软件包自带的绝大多数赝势是超软赝势（US-PP）了，但不少元素有两个版本，如何 选取呢？

一个简单的办法是看后缀

标准的没有后缀 _h 硬一点 _s 软一点

_pv,_sv,_d 就是说semi-core的p,s或者d也当做价态处理了 如果是数字的话，表示的可能是不同的半径截距

也可以参考各版本同目录下的V_RHFIN file ，PSCTR file 这两个文件告知该版本的赝势是如何生成的。比如： V_RHFIN file Sc: 6p d2 s1

8 21. .002000 44.95590 125. .2

VASP和CASTEP的参数设置区别

请问大家：再用VASP计算的时候，K点应该怎么选取呢？我用CASTEP默认的那个设置算的，有的会提示说 VERY BAD NEWS! internal error in subroutine IBZKPT: Tetrahedron method fails for NKPT<4. NKPT =2：貌似是说我的K点设置的不合理吧，那我应该用什么办法选取K点呢？在文献上找不到跟我的模型一样的，没办法参考别人的。 第二个问题是计算能带结构的时候高对称点有应该怎么选取呢， 建立不同大小的超晶胞选用的截断能和高对称点应该一样吗？

我一直搞不清楚这个问题，就是VASP计算的时候选用的势函数、截断能可以按照CASTEP默认的来吗？我自己测试的截断能和CASTEP默认的怎么相差特别大呢？

1. VASP中采用Tetrahedron method 时，k点不能少于4个。解决方法：可以选ISMEAR=0;SIGMA=0.1

2. 计算能带时的高堆成点可用CASTEP提供的，即CASTEP计算能带时用的k特殊点。不同超晶胞选用的截断能不能一样啊，要不然超晶胞太大导致无法运算，不是吗

3. vasp计算的时候选用的势函数、截断能可以用CASTEP默认的来，但是不一定精确，所以一般还是要自己测试之后才用。以测试得到的结果为准。

谢谢你的回答，我还想在请问一下，你的意思就是说：VASP就直接采用CASTEP默认的K点吗？我看好多人用的跟默认的还不一样，这是怎么选的呢？

还有

实例1

VASP赝势文件POTCAR详解 VASP赝势文件POTCAR详解

POTCAR 赝势文件

可以理解为分子力学模拟中的力场文件 但包括的信息更多 VASP4.6将各元素优化的INCAR里的参数也包括在这里了，作为支持PREC的缺省选择 通常各元素的POTCAR已经包括在软件包里了

我们只需要按照POSCAR里的顺序，将各元素的POTCAR按顺序连接起来就可以了 如以下命令：

cat file1 file2 file3 > POTCAR

软件包自带的绝大多数赝势是超软赝势（US-PP）了，但不少元素有两个版本，如何 选取呢？

一个简单的办法是看后缀

标准的没有后缀 _h 硬一点 _s 软一点

_pv,_sv,_d 就是说semi-core的p,s或者d也当做价态处理了 如果是数字的话，表示的可能是不同的半径截距

也可以参考各版本同目录下的V_RHFIN file ，PSCTR file 这两个文件告知该版本的赝势是如何生成的。比如： V_RHFIN file Sc: 6p d2 s1

8 21. .002000 44.95590 125. .2

VASP计算的理论及实践总结

一、赝势的选取 二、收敛测试

1、VASP测试截断能和K 点 2、MS测试

三、结构弛豫

四、VASP的使用流程（计算性质）

1、VASP的四个输入文件的设置 2、输出文件的查看及指令 3、计算单电能

(1) 测试截断能 (2) 测试K点 4、进行结构优化 5、计算弹性常数 6、一些常用指令

一、赝势的选取

VASP赝势库中分为：PP和PAW两种势，PP又分为SP（标准）和USPP（超软）。

交换关联函数分为：LDA（局域密度近似）和GGA（广义梯度近似）。GGA又分为PW91和PBE。

在VASP中，其中pot ,pot-gga是属于超软势（使用较少）。Paw, paw-pbe ,和paw-gga是属于PAW。采用较多的是PAW-pbe 和PAW-gga。

此外vasp 中的赝势分为几种，包扩标准赝势(没有下标的)、还有硬（harder）赝势(_h)、软（softer）赝势(_s), 所谓的硬（难以赝化），就是指该元素原子的截断动能比较大，假想的势能与实际比较接近，计算得到的结果准确，但比较耗时，难以收敛。软（容易赝化），表示该元素原子的截断动能比较小，赝势模型比较粗糙，但相对简单，可以使计算很快收

VASP5.2安装心得

2014-05-07 来源：小木虫作者： yysskk

花了五天时间终于学会怎么装VASP了，在此写下心得体会，供后人参考。个人觉得最难的一步就是makefile文件，网上流传着各种各样的版本，每个人都说自己编译成功了，却又各不相同，也说不清为什么，给新手极大的困扰。在此会详细介绍makefile的文件结构。其余大部分内容都是参考前人的，就不一一注明出处了。

一、系统、编译程序及准备工作

我用的是centos6.5+icc2011+ifort2011+openmpi1.6.5 1.1编译器安装

系统安装不说了，网上教程多得是。Icc和ifort可以申请免费非商业版本，icc和ifort都各自带了一个MKL，使用的时候别搞混了。装2011的时候会缺组件，用yum都可以免费下载。装编译器的时候会要求关闭selinux，按照给出的步骤关闭即可。之后会说系统不兼容，但是可以继续装，默认安装路径是在/opt下面。装完之后会有提示，把安装目录/bin/ifortvars.sh 写到环境变量中，注意32/64位系统的参数不一样。C语言编译器建议用icc，毕竟是intel出品，针对自家cpu肯定有大量优化，效率上高于gcc是肯定的。

VASP计算的理论及实践总结

一、赝势的选取 二、收敛测试

1、VASP测试截断能和K 点 2、MS测试

三、结构弛豫

四、VASP的使用流程（计算性质）

1、VASP的四个输入文件的设置 2、输出文件的查看及指令 3、计算单电能

(1) 测试截断能 (2) 测试K点 4、进行结构优化 5、计算弹性常数 6、一些常用指令

一、赝势的选取

VASP赝势库中分为：PP和PAW两种势，PP又分为SP（标准）和USPP（超软）。

交换关联函数分为：LDA（局域密度近似）和GGA（广义梯度近似）。GGA又分为PW91和PBE。

在VASP中，其中pot ,pot-gga是属于超软势（使用较少）。Paw, paw-pbe ,和paw-gga是属于PAW。采用较多的是PAW-pbe 和PAW-gga。

此外vasp 中的赝势分为几种，包扩标准赝势(没有下标的)、还有硬（harder）赝势(_h)、软（softer）赝势(_s), 所谓的硬（难以赝化），就是指该元素原子的截断动能比较大，假想的势能与实际比较接近，计算得到的结果准确，但比较耗时，难以收敛。软（容易赝化），表示该元素原子的截断动能比较小，赝势模型比较粗糙，但相对简单，可以使计算很快收

用VASP计算Ｈ原子的能量

氢原子的能量为 。在这一节中，我们用VASP计算H原子的能量。对于原子计算，我们可以采用如下的INCAR文件 PREC=ACCURATE

NELMDL = 5 make five delays till charge mixing ISMEAR = 0; SIGMA=0.05 use smearing method 采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。 Monkhorst Pack 0 Monkhorst Pack 1 1 1 0 0 0

采用如下的POSCAR文件 atom 1

15.00000 .00000 .00000 .00000 15.00000 .00000 .00000 .00000 15.00000 1 cart

0 0 0

采用标准的H的POTCAR 得到结果如下：

k-point 1 : 0.0000 0.0000 0.0