POTCAR 赝势产生理论.docx

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 fileSc: 6p d2 s18 21..002000 44.95590 125. .25E-05 .300 200FCA 12.00000.7 1.0 01.0 .0 .5 -320.8847 2.00002.0 .0 .5 -34.4217 2.00002.0 1.0 1.5 -28.2366 6.00003.0 .0 .5 -3.7944 2.00003.0 1.0 1.5 -2.2591 6.00003.0 2.0 2.5 -.1113 2.00004.0 .0 .5 -.2699 1.00004.0 3.0 2.5 -.1000 .0000第一行是注释行给出基本的信息 第二行是最重要的控制行8 21..00 2000 44.95590 125. .25E-05 .300 200 F CA 12.00000J Z XION N AM H DELRVR PHI NC1 | CH QCOR|GREENJ - 轨道数 Z - 原子序数 XION - 离子化程度一般设为 0 N - 格点数AM - 原子质量 H - 决定格点间距 DELRVR - 自洽收敛标准PHI - 线性拟合参数 NC1 - 最大自洽循环次数 GREEN - 是否存在初始的势CH -交换相关能(XC)类型Slater-XCHL Hedin Lundquist (1971)CA Ceperly and Alder parameterized by J.Perdew and ZungerWI Wigner interpolationPB Perdew -BeckePW Perdew -Wang 86LM Langreth-Mehl-Hu91 Perdew -Wang 91QCOR -非价键电子数(core electro ns)第三行开始是每个轨道的具体参数，依次为n l j(=l ±1/2) 原子轨道能占有率PSCTR file of LDA/H1.25TITEL = US HLULTRA = T use ultrasoftPP ?RWIGS = 0.57 nndistance ! Wigner-Seitz radiusRCLOC = .65NE = 100LCOR = .TRUE.QCUT = -1RMAX = 3.0 !core radius for proj-operDescriptionl E TYP RCUT TYP RCUT(cutoff radius)0 0 15 0.80 23 1.250 0.5 15 0.80 23 1.251 -0.2 15 0.80 23 1.25 最重要的地方上面已经用颜色标出来啦：)说明一下，TYP是指贋势的类型,RCUT是半径截距，TYP可取的值如下： 正则1 BHS2 TM3 VAN6 XNC7 RRKJ wave function possibly with node15 RRKJ wave function strictly no node 非正则 +8最后一个问题是LDA or GGA。

貌似没有定论目前 这个最好是两个一起做做看啦或者看文献别人验证过哪个数据好其实据说目前最好的是 PAW(P.E.Blochl,Phys.Rev.B 50,17953(1994).,Phys.Rev.B59,1758(1999).)。