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

分享：