24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2630  |  回复: 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的回帖

souledge

专家顾问 (著名写手)

★ ★
zzy870720z(金币+2): 鼓励讨论 2011-10-10 09:23:14
引用回帖:
20楼: Originally posted by 改王之王 at 2011-10-09 17:25:12:
昨天试了试这个…………还是不太好…………
label1, occupancy > 4s 2
label1, occupancy > 3d 10
label1, occupancy > 4p 0
label1, occupancy > 4d 0

别试了,P.Giannozzi已经说过了,那个fhi2upf.x的代码中,很多转换部分都是猜测级的~虽然他认为那些猜测都还reasonable,但是不敢保证一定可以,除非你能问到初试的生成参数。
另外,大多数fhi的PP都是multi-projector(不知道是不是这么个名称,就是考虑到很多不同的激发态之类的)的,所以偶尔参照态就算猜准了,也不见得转换出来一定“源汁源味”~
思想重于技巧,内涵重于表象
21楼2011-10-09 18:33:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

改王之王

金虫 (正式写手)

发现我输入赝势名称的时候提示
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的回帖
信息提示
请填处理意见