应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » abinit中赝势的含义
北京石油化工学院2026年研究生招生接收调剂公告
231/3返回列表123下一页
查看: 2922  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

改王之王

金虫 (正式写手)

[求助] 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 相关资料

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-13 20:34:55
已阅   回复此楼   关注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也这么搞~
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
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，估计要自己生成了……
赞一下(1人) 回复此楼
思想重于技巧，内涵重于表象
10楼2011-09-14 00:08:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 改王之王 的主题更新
231/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 chenxrlkx 2026-04-05 6/300 2026-04-05 22:05 by imissbao
[考研] 285求调剂 +4 恶法大二的气味� 2026-04-05 5/250 2026-04-05 20:32 by 286640313
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 调剂 +5 好好读书。 2026-04-01 5/250 2026-04-05 17:54 by liucky
[考研] 385分 生物学（071000）求调剂 +11 qf626 2026-04-01 11/550 2026-04-05 17:35 by Ecowxq666！
[考研] 304求调剂（085602，过四级，一志愿985） +15 化工人999 2026-04-04 15/750 2026-04-05 16:25 by 我是电风扇r
[考研] 085410人工智能 初试316分 求调剂 +6 残星拂曙 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 考研调剂 +6 15615482637 2026-04-04 6/300 2026-04-04 22:43 by yu221
[考研] 一志愿武理材料工程302调剂环化或化工 +19 Doleres 2026-03-31 20/1000 2026-04-04 16:44 by 啊俊！
[考研] 求调剂，一志愿北京中医药大学 +3 小小达不溜 2026-04-02 3/150 2026-04-03 22:55 by 冲矢昴星团
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 机械专硕297 +3 Afksy 2026-04-03 3/150 2026-04-03 14:24 by 1753564080
[考研] 工科341分调剂 +3 洛多罗 2026-04-03 3/150 2026-04-03 14:20 by 1753564080
[考研] 302求调剂 +9 zyx上岸！ 2026-04-02 9/450 2026-04-02 23:07 by 马儿快快地跑
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业 +6 what张 2026-04-01 7/350 2026-04-02 16:48 by zzsw+
[考研] 275学硕081000服从调剂到其他专业，保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
[考研] 070300化学279求调剂 +15 哈哈哈^_^ 2026-03-31 17/850 2026-04-01 21:37 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭