24小时热门版块排行榜    

查看: 2874  |  回复: 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
回复此楼
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[文学芳草园] 伙伴们,祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
信息提示
请填处理意见