| 查看: 2468 | 回复: 2 | ||
fangyongxinxi新虫 (正式写手)
|
[求助]
求助：详细解释vasp势函数库中，各文件的含义。
|
|
自己对MS有一定的操作经验，在MS里面对选择GGA,LDA之类的东东，觉得很方便，刚学vasp，但在vasp中就不知道怎么办了。 如题，自己把vasp的势函数打开，如附件所示，出现 pot pot_GGA potpaw potpaw_GGA potpaw_PBE （问题1. lda 在什么地方 ？） 进入pot---element 会出现各个元素的势文件 如 Ag Al C C_s ( 问题2. 这个C_s 是什么意思 ？什么时候用C 什么时候用 C_s ？有时候还会有_pv, _d , _s, _200eV, 不知道什么意思) 进入Ag 出现 DDE POTCAR PSCTR V_RHFIN V_TABIN( 问题3. 除了POTCAR知道有用外，其他几个文件不知道什么意思，有何用处)  |
回复此楼» 猜你喜欢
材料专硕283求调剂
已经有7人回复
-
一志愿北京2，材料与化工308求调剂
已经有4人回复
-
求调剂
已经有5人回复
-
一志愿9材料学硕297已过六级求调剂推荐
已经有8人回复
-
278求调剂
已经有9人回复
-
285求调剂
已经有5人回复
-
一志愿武理材料工程302调剂环化或化工
已经有20人回复
-
348求调剂
已经有3人回复
-
本科211 分数293请求调剂
已经有3人回复
-
一志愿南农090401，268，求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何设定vasp中交换势
已经有5人回复
- 求助： vasp中化合物势函数的选取。
已经有6人回复
- 【求助】VASP中H2O的势函数［求助成功］
已经有4人回复
- 【求助】VASP中Zn的化学势是多少？
已经有15人回复
- 【求助】vasp中COOP如何理解及输出文件如何处理？
已经有12人回复
- 【求助】怎么建立POTCAR，vasp中找不到赝势文件
已经有4人回复
- 【求助】VASP中K点设置
已经有14人回复
- 【求助】求教VASP并行版的含义？
已经有6人回复
- 【求助】vasp检验Ti的赝势kpoint如何设置
已经有13人回复
- 【求助】vasp的计算时间含义
已经有12人回复
- 【求助】求VASP的赝势库。。
已经有11人回复
- 【求助】请问vasp有B3LYP赝势吗？
已经有7人回复
漂泊四方
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 4 (幼儿园)
- 金币: 1321.2
- 散金: 50
- 红花: 6
- 帖子: 198
- 在线: 115.2小时
- 虫号: 897145
- 注册: 2009-11-08
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
franch(金币+3): 谢谢回帖交流 2011-05-31 09:32:20
franch(金币+3): 谢谢回帖交流 2011-05-31 09:32:20
|
纯复制，精华贴里面有很多类似的东西，可以找找看，有不少惊喜。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 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. .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)类型 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 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 |
2楼2011-05-30 22:53:52
fangyongxinxi
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 41
- 散金: 69
- 红花: 3
- 帖子: 876
- 在线: 138.5小时
- 虫号: 437418
- 注册: 2007-08-26
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
3楼2011-05-30 23:12:12
