vasp 赝势选取

http://blog.sina.com.cn/s/blog_17589e0920102xinl.html
http://blog.sina.com.cn/s/blog_17589e0920102xinl.html
1.VASP赝势按产生方法可以分为PP (standard pesudopotential，其中大部分是USPP, ultrasoft pesudopotential) 和 PAW (projector augmented wave method)。按交换关联函数的不同又可以有 LDA (local density approximation) 和 GGA (generalized gradient approximation)，其中GGA之下又可以再分为PW91和PBE。
vasp赝势包里有 pot， pot_GGA， potpaw， potpaw_GGA ， potpaw_PBE。
以上各个目录对应起来分别是 pot ==> PP, LDA ; pot_GGA ==> PP, GGA ; potpaw ==> PAW, LDA ; potpaw_GGA ==> PAW, GGA, PW91 ; potpaw_PBE ==> PAW , GGA, PBE。选择某个目录进去，我们还会发现对应每种元素往往还会有多种赝势存在。标准的没有后缀；根据对截断能量的选取不同还可以分为Ga,Ga_s,Ga_h， _h 硬一点 _s 软一点；或者根据半芯态的不同还可以分为Ga,Ga_sv,Ga_pv的不同， _pv,_sv,_d 就是说semi-core的p,s或者d也当做价态处理了。
对于H元素，赝势库 potpaw—PBE 里面 除了H， H_s， H_h 之外 还有H1.25 H1.33 H1.5 H1.66 H1.75 H.25 H.33 H.42 H.5 H.58 H.66 H.75 这个只有H原子有这些 其它元素都没有。H后面的数字代表电荷数目。普通计算选1就可以，其它几个带数字的分数氢是表面计算用的，中和体系电荷用的。H_h是高压用的，如果是常压下的话，推荐用H，H_s有时候算出来结果不准。
一般推荐选取PAW_PBE。
2.POTCAR的建立：
选好哪一种赝势之后，进入对应的目录，你会看到里边有这么几个文件，POTCAR.Z PSCTR.Z V_RHFIN.Z WS_FTP.LOG 。我们需要的是第一个。把它解压，如 zcat POTCAR.Z > Ga 。对As元素我们也可以类似得到一个As文件。用 cp 命令或者 mv 命令把这两个文件都移到我们的工作目录里。然后再用 cat 命令把这两个文件合并在一起，如 cat Ga As > POTCAR ，这样就得到了我们需要的POTCAR。同理，有多个元素的POTCAR也可以这样产生。这里需要注意的是，记住元素的排列顺序，以后在POSCAR里各个元素的排列就是按着这里来的。
查看POTCAR里的截断能量为多大，用 grep -in “enmax” POTCAR 命令。
3.官方推荐赝势
VASP官方是有一些推荐选择的赝势。详情可以访问VASP的官方manual，
http://cms.mpi.univie.ac.at/vasp/vasp/Recommended_PAW_potentials_DFT_calculations_using_vasp_5_2.html 其中加粗显示的即为推荐的。下面这个表将其中推荐的挑出来了。
Number Element ENMAX (eV) valency
1 H 250 1
2 He 479 2
3 Li sv 499 3
4 Be 248 2
5 B 319 3
6 C 400 4
7 N 400 5
8 O 400 6
9 F 400 7
10 Ne 344 8
11 Na pv 260 7
12 Mg 200 2
13 Al 240 3
14 Si 245 4
15 P 255 5
16 S 259 6
17 Cl 262 7
18 Ar 266 8
19 K sv 259 9
20 Ca sv 267 10
21 Sc sv 223 11
22 Ti sv 275 12
23 V sv 264 13
24 Cr pv 266 12
25 Mn pv 270 13
26 Fe 268 8
27 Co 268 9
28 Ni 270 10
29 Cu 295 11
30 Zn 277 12
31 Ga d 283 13
32 Ge d 310 14
33 As 209 5
34 Se 212 6
35 Br 216 7
36 Kr 185 8
37 Rb sv 220 9
38 Sr sv 229 10
39 Y sv 203 11
40 Zr sv 230 12
41 Nb sv 293 13
42 Mo sv 243 14
43 Tc pv 264 13
44 Ru pv 240 14
45 Rh pv 247 15
46 Pd 251 10
47 Ag 250 11
48 Cd 274 12
49 In d 239 13
50 Sn d 241 14
51 Sb 172 5
52 Te 175 6
53 I 176 7
54 Xe 153 8
55 Cs sv 220 9
56 Ba sv 187 10
57 La 219 11
58 Ce 273 12
59 Pr 3 182 11
60 Nd 3 183 11
61 Pm 3 177 11
62 Sm 3 177 11
63 Eu 2 99 8
64 Gd 3 154 9
65 Tb 3 156 9
66 Dy 3 156 9
67 Ho 3 154 9
68 Er 3 155 9
69 Tm 3 149 9
70 Yb 2 113 8
71 Lu 3 155 9
72 Hf pv 220 10
73 Ta pv 224 11
74 W pv 223 12
75 Re 226 7
76 Os 228 8
77 Ir 211 9
78 Pt 230 10
79 Au 230 11
80 Hg 233 12
81 Tl d 237 13
82 Pb d 238 14
83 Bi d 243 15
84 Po d 265 16
85 At d 266 17
86 Rn 152 8
87 Fr sv 215 9
88 Ra sv 237 10
89 Ac 172 11
90 Th 247 12
91 Pa 252 13
92 U 253 14
93 Np 254 15
94 Pu 254 16
95 Am 256 17
96 Cm 258 18

分享：