精华区文章阅读

发信人: zjliu (秋天的萝卜), 信区: Chemistry
标  题: 初学VASP(五) 赝势·选择
发信站: BBS 哈工大紫丁香站 (Thu Jun 17 20:47:26 2004)

发信站: 日月光华

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
V_RHFIN file
Sc: 6p d2 s1
    8  21.   .002000  44.95590 125.   .25E-05 .300  200FCA 12.00000
   .7  1.0   0
  1.0   .0   .5 -320.8847 2.0000
  2.0   .0   .5  -34.4217 2.0000
  2.0  1.0  1.5  -28.2366 6.0000
  3.0   .0   .5   -3.7944 2.0000
  3.0  1.0  1.5   -2.2591 6.0000
  3.0  2.0  2.5    -.1113 2.0000
  4.0   .0   .5    -.2699 1.0000
  4.0  3.0  2.5    -.1000  .0000
第一行是注释行给出基本的信息
第二行是最重要的控制行
    8  21.   .00 2000  44.95590 125.   .25E-05 .300  200 F CA 12.00000
    J   Z   XION    N  AM         H    DELRVR   PHI  NC1 | CH QCOR
                                                         |
                                                     GREEN

J - 轨道数  Z - 原子序数    XION - 离子化程度一般设为0 N - 格点数
AM - 原子质量   H - 决定格点间距    DELRVR - 自洽收敛标准
PHI - 线性拟合参数 NC1 - 最大自洽循环次数 GREEN - 是否存在初始的势
CH - 交换相关能(XC)类型
CH - 交换相关能(XC)类型
      Slater-XC
    HL Hedin Lundquist (1971)
    CA Ceperly and Alder parameterized by J.Perdew and Zunger
    WI Wigner interpolation
    PB Perdew -Becke
    PW Perdew -Wang 86
    LM Langreth-Mehl-Hu
    91 Perdew -Wang 91
QCOR - 非价键电子数（core electrons）
第三行开始是每个轨道的具体参数，依次为
n l j(=l±1/2) 原子轨道能占有率

PSCTR file of LDA/H1.25
   TITEL  = US H
   LULTRA =    T    use ultrasoft PP ?
   RWIGS  =    0.57  nn distance ! Wigner-Seitz radius

   RCLOC  =   .65
   NE     =   100
   LCOR   =   .TRUE.
   QCUT   =   -1
   RMAX   = 3.0 ! core radius for proj-oper
   RMAX   = 3.0 ! core radius for proj-oper

   Description
l     E      TYP RCUT    TYP  RCUT(cutoff radius)
0   0         15   0.80  23   1.25
0   0.5       15   0.80  23   1.25
1  -0.2       15   0.80  23   1.25
最重要的地方上面已经用颜色标出来啦：）
说明一下，TYP是指赝势的类型,RCUT是半径截距，TYP可取的值如下：
正则
1 BHS
2 TM
3 VAN
6 XNC
7 RRKJ wave function possibly with node
15 RRKJ wave function strictly no node
非正则 +8

最后一个问题是LDA or GGA。貌似没有定论目前。
这个最好是两个一起做做看啦。或者看文献别人验证过哪个数据好。
其实据说目前最好的是PAW(P.E.Blochl,Phys.Rev.B 50,17953(1994).,Phys.Rev.B 59,
1758(1999).)，PP已经落伍了，不过好像我们用的VASP不带这个，就不展开讨论了。
--
╔═══════════════════╗
║★★★★★友谊第一  比赛第二★★★★★║
╚═══════════════════╝

※ 来源:·哈工大紫丁香 http://bbs.hit.edu.cn·[FROM: 202.118.229.*]

Chemistry 版 (精华区)