24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2637  |  回复: 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(金币+1): 鼓励交流 2011-09-14 17:47:01
改王之王(金币+10): 感谢回复………… 2011-09-14 19:19:13
  分析对了,果然是fhi的赝势本身提供信息太少,所以转化起来很……这封P.Giannozzi的回信估计会比较有用吧:
引用回帖:
shudong wang wrote:

> I want to use fhi2upf.x to convert 27-Co.LDA.fhi to upf format, but I
> have some questions:
> 1.What the wavefunction1-4 should I write?

you should ask this to whoever generated that pseudopotential.
The FHI format does not contain all the info needed to be used
with other codes.

You should in any case know how many valence electrons there
are. Occupancies are used only for starting the calculation
and are not critical. The correct angular momentum is important.
I think it is s,p,d,f, in this order, but no warranty.

> If my converted pp is not precise, how can I solve this? Which
> parameters should I change to enhance the precision??

this is a big question that doesn't have a quick answer.
Please have a look at the atomic_doc/pseudo-gen.pdf file

P.
--
Paolo Giannozzi, Democritos and University of Udine, Italy

思想重于技巧,内涵重于表象
11楼2011-09-14 00:38:10
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见