24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2617  |  回复: 22

改王之王

金虫 (正式写手)

[求助] abinit中赝势的含义

最近想算一下声子的raman性质,但是PWSCF里面没有Zn的NCPP,只有USPP赝势,所以自己从abinit里面下了个Zn-LDA.fhi赝势,用pwscf里面的fhi2upf.x来转换,但是转换的时候第一步提示label1,occupancy >
我理解是要电子的占据数吗?看了看下面的赝势文件,我输入了6然后回车,但是没反应…………求指教

赝势文件内容:
zinc, fhi98PP : Trouiller-Martins-type, LDA Ceperley/Alder Perdew/Wang (1992), l= 0 local
30.000 12.000    011001              zatom,zion,pspdat
   6         7    3   0    527     0       pspcod,pspxc,lmax,lloc,mmax,r2well
  0.000  0.000  0.000                 rchrg, fchrg, qchrg
5--- These two lines are available for giving more information, later
6
7-Here follows the cpi file from the fhi98pp code-
0.12000000000000E+02   4
  0.0000    0.0000    0.0000   0.0000
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
527  0.10247000000000E+01
   1 0.20833333333333E-03 0.32462165437455E-04 -.98924907634198E+01
   2 0.21347916666667E-03 0.33263980923760E-04 -.98924908039977E+01
   3 0.21875210208333E-03 0.34085601252577E-04 -.98924908435820E+01
   4 0.22415527900479E-03 0.34927515603515E-04 -.98924908821959E+01
   5 0.22969191439621E-03 0.35790225248781E-04 -.98924909198620E+01
   6 0.23536530468180E-03 0.36674243828480E-04 -.98924909566024E+01
   7 0.24117882770744E-03 0.37580097672687E-04 -.98924909924383E+01
   8 0.24713594475181E-03 0.38508326113058E-04 -.98924910273908E+01
   9 0.25324020258718E-03 0.39459481802581E-04 -.98924910614800E+01
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 Abinit 相关资料

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

发现我输入赝势名称的时候提示
pseudopential without NLCC successfully read
2楼2011-09-13 20:50:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

明白了,是电子轨道占据数的意思,可是根据什么添呢?是根据元素本身电子占据方式吗?和赝势文件内容没关系?
3楼2011-09-13 21:08:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

我是按照元素周期表这样写的
label1, occupancy > 4s 2
label1, occupancy > 3d 10
label1, occupancy > 3p 6
label1, occupancy > 3s 2
4楼2011-09-13 21:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

不知道这样写成不成,我待会儿算个声子的色散关系和以前对比一下…………
5楼2011-09-13 21:30:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

但是O的赝势我这样输入的
label1, occupancy > 2p 2
label1, occupancy > 2s 2
label1, occupancy > 1s 2
label1, occupancy > 1s 2
不知道这样对不对,O只有三个轨道电子,可是要输入四次,最后一次我只能输入1s了,不知道对不对………………
6楼2011-09-13 21:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

我按照以上输入计算自洽后提示这个错误:
Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
     file Zn.LDA.fhi.UPF: wavefunction(s)  3d 3s renormalized
7楼2011-09-13 22:15:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)

完整的错误文件:
Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
     Waiting for input...
     file Zn.LDA.fhi.UPF: wavefunction(s)  3d 3s renormalized

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from readpp : error #         2
     inconsistent DFT read
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
8楼2011-09-13 22:39:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
贺仪(金币+5): 多谢参与! 2011-09-14 09:15:32
引用回帖:
6楼: Originally posted by 改王之王 at 2011-09-13 21:49:26:
但是O的赝势我这样输入的
label1, occupancy > 2p 2
label1, occupancy > 2s 2
label1, occupancy > 1s 2
label1, occupancy > 1s 2
不知道这样对不对,O只有三个轨道电子,可是要输入四次,最后 ...

我只能说FHI赝势的文件头很坑爹~这是氧的:
CODE:
oxygen, fhi98PP : Trouiller-Martins-type, LDA Ceperley/Alder Perdew/Wang (1992), l= 2 local
  8.000  6.000    021003              zatom,zion,pspdat
   6         7    3   2    473     0       pspcod,pspxc,lmax,lloc,mmax,r2well
  0.000  0.000  0.000                 rchrg, fchrg, qchrg
5--- These two lines are available for giving more information, later
6
7-Here follows the cpi file from the fhi98pp code-
6.00000E+00    4
0.0000    0.0000    0.0000   0.0000
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00
  0.0000    .00e+00   .00e+00

我实在很怀疑,对于只有1s2s2p轨道的氧,是怎么搞出来了lmax=3,lloc=2的……要知道,lmax=3已经是f轨道~而lloc=2也都是d轨道了。
根据官方说明:
引用回帖:
the maximal angular momentum of the wavefunctions described in the pseudopotential file (lmax=0 if only s-wavefunctions are present, lmax=1 for s and p wfs, lmax=2 for s, p and d wfs ... ; lmax+1 must be equal to the 'number of pseudopotential components' mentioned in the cpi file, as the second number of the first line)

又要lmax表示最大角动量,又要其等于赝势组分数~真晕~
但是lloc,就实在说不过去了……
引用回帖:
the angular momentum of the potential to be used as local pseudopotential (lloc=0 if s-potential is local; lloc=1 if p-potential is local ; ... )

无论如何,也不能用d轨道作local吧~氧去哪儿找d轨道去啊……
希望上面只是我错误的理解,顺便,鉴于PW的赝势库中,已有赝势的文件头部分都给改得面目全非,反正不影响使用~不排除abinit也这么搞~
思想重于技巧,内涵重于表象
9楼2011-09-14 00:05:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖


zzy870720z(金币+1): 谢谢指点 2011-09-14 17:46:49
引用回帖:
2楼: Originally posted by 改王之王 at 2011-09-13 20:50:36:
发现我输入赝势名称的时候提示
pseudopential without NLCC successfully read

这是因为早期的Abinit不支持读取FHI赝势中的非线性核矫正(不知道是不是这么翻译)部分,所以Abinit中的FHI赝势大多都不带NLCC的~如果需要NLCC,估计要自己生成了……
思想重于技巧,内涵重于表象
10楼2011-09-14 00:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
信息提示
请填处理意见