首页     专家     文章          软件     小组     视频    

注册  |  登录
首页 > 文章 > VASP入门系列文章[5]初识POTCAR

VASP入门系列文章[5]初识POTCAR

作者:obaica   (在线)   [VASP ]   [POTCAR ]   时间:2016-09-04 22:48:33  向他请教

工具/原料:VASP,POTCAR

步骤/方法

上一讲链接http://www.goscience.cn/article/articleInfo?type=2&id=465

今天继续讲解POTCAR输入文件,

生成POTCAR的命令是:

cat POTCAR_file1 POTCAR_file2 POTCAR_file3 > POTCAR

需要注意的是:

(1) 这里元素的顺序要和POSCAR里的顺序一一对应起来,

(2) 在选择赝势的时候,注意文件末尾字符的区别,标准的没有后缀 _h 硬一点 _s 软一点_pv,_sv,_d 就是说semi-core的p,s或者d也当做价态处理了,如果是数字的话,表示的可能是不同的半径截距。

———————————————————————————————————

 PAW_PBE Ag 02Apr2005

   11.0000000000000

 parameters from PSCTR are:

   VRHFIN =Ag : s1 d10

   LEXCH  = PE

   EATOM  =  1037.2675 eV,   76.2370 Ry


   TITEL  = PAW_PBE Ag 02Apr2005

   LULTRA =        F    use ultrasoft PP ?

   IUNSCR =        1    unscreen: 0-lin 1-nonlin 2-no

   RPACOR =    2.000    partial core radius

   POMASS =  107.868; ZVAL   =   11.000    mass and valenz

   RCORE  =    2.500    outmost cutoff radius

   RWIGS  =    2.840; RWIGS  =    1.503    wigner-seitz radius (au A)

   ENMAX  =  249.844; ENMIN  =  187.383 eV

   ICORE  =        3    local potential

   LCOR   =        T    correct aug charges

   LPAW   =        T    paw PP

   EAUG   =  412.496

   DEXC   =    0.000

   RMAX   =    2.560    core radius for proj-oper

   RAUG   =    1.300    factor for augmentation sphere

   RDEP   =    2.576    radius for radial grids

   RDEPT  =    2.072    core radius for aug-charge


   Atomic configuration

   12 entries

     n  l   j            E        occ.

     1  0  0.50    -25324.9161   2.0000

     2  0  0.50     -3728.5341   2.0000

     2  1  1.50     -3344.0404   6.0000

     3  0  0.50      -683.5581   2.0000

     3  1  1.50      -557.3225   6.0000

     3  2  2.50      -358.3853  10.0000

     4  0  0.50       -95.0341   2.0000

     4  1  1.50       -58.5511   6.0000

     4  2  2.50        -7.4440  10.0000

     5  0  0.50        -4.4045   1.0000

     5  1  0.50        -1.3606   0.0000

     4  3  2.50        -1.3606   0.0000

   Description

     l       E           TYP  RCUT    TYP  RCUT

后面略

———————————————————————————————————

VRHFIN =Ag : s1 d10   外层价电子,可以参考网站http://www.ptable.com/#Orbital

LEXCH  = PE 一般是对应的PBE赝势,如果LEXCH  = CA对应的是LDA赝势。

注意ENMAX  =  249.844; ENMIN  =  187.383 eV的取值


关于每个参数的详细含义可以参考VASP手册

http://cms.mpi.univie.ac.at/wiki/index.php/The_VASP_Manual

———————————————————————————————————

有兴趣的可以加入我的小组进行交流和分享经验: 

http://www.goscience.cn/group/groupInfo?id=124 (第一性原理计算交流小组)

———————————————————————————————————

好了,下一节内容说一下其他输入文件~欢迎各位批评指正和点赞!

发表提问

您有0条新消息