24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2920  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08600生物与医药-327 +5 18755400796 2026-04-05 5/250 2026-04-05 21:05 by 学员8dgXkO
[考研] 338求调剂 +3 我想上岸ii 2026-04-05 3/150 2026-04-05 19:59 by nepu_uu
[考研] 298求调剂 +7 manman511 2026-04-05 7/350 2026-04-05 10:29 by 唐沐儿
[考研] +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考研] 302分求调剂 一志愿安徽大学085601 +5 zyx上岸! 2026-04-04 5/250 2026-04-05 07:50 by 544594351
[考研] 0703求调剂383分 +8 W55j 2026-04-03 8/400 2026-04-04 20:09 by xhai2011
[考研] 一志愿C9的化学工程(085602) 340分,感觉校内调剂无望,求调剂 +9 万事宜臻 2026-04-04 9/450 2026-04-04 11:49 by 啵啵啵0119
[考研] 278求调剂 +6 Yy7400 2026-04-03 6/300 2026-04-04 09:53 by zhangdingwa
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 一志愿深大085601材料工程专业(专硕)300分可以调剂去哪 +8 10160315 2026-04-02 8/400 2026-04-03 09:36 by hypershenger
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 085600,320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 262求调剂 +9 励志一定发文章 2026-03-31 10/500 2026-04-01 12:22 by sunshine0013
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
信息提示
请填处理意见