24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

[求助] 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

回复此楼

» 猜你喜欢

Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
2025冷门绝学什么时候出结果已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2011-09-13 20:34:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
散金: 61
红花: 3
帖子: 817
在线: 113.5小时
虫号: 968169
注册: 2010-03-11
专业: 凝聚态物性 II ：电子结构

★ ★
zzy870720z(金币+2): 鼓励讨论 2011-10-09 11:30:45

引用回帖:

17楼: Originally posted by souledge at 2011-10-08 17:32:39:

做了新的Zn和O的PBE的模守赝势，要不要来试试，是不是比你自己转化的那个好啊，话说，100Ry的截断才能勉强的收敛，实在有点过吧～～～～～

你的赝势截断能设到47ry就收敛的比较好了，和我算的87ry的结果相当，可能是我的赝势转化的有错误，我是按下面的方法做的
Zn
label1, occupancy > 3d 10
label1, occupancy > 4s 2
label1, occupancy > 4p 0
label1, occupancy > 4f 0

O
label1, occupancy > 2s2
label1, occupancy > 2p 4
label1, occupancy > 3s 0
label1, occupancy > 3p 0

然后我去maillist搜了一下看到paolo说的这句话
On Tuesday 14 April 2009 05:19, shangy wrote:

> When I run the fhi2upf.x, it asked the "Wavefunction label and occupancy"
> [...] I tried it like this
> ---------------------
> Wavefunction # 1: label, occupancy > 4s 2
> Wavefunction # 2: label, occupancy > 4p 2
> Wavefunction # 3: label, occupancy > 4d 0
> Wavefunction # 4: label, occupancy > 4f 0
> ---------------------
> but I don't know if it is right?

the FHI format doesn't contain - to the best of my very limited
knowledge - any information on the electronic configuration
used to generate the pseudopotential. Atomic orbitals are
actually used only to produce the starting wavefunctions
and charge density for self-consistency. You should specify
the labels of l=0 to lmax valence orbitals (in the order:s,p,d,f,..)
and a reasonable value for their occupancy. You have Ge with
4 valence electrons, so your choice seems to be correct.

我感觉是不是我Zn的设错了？4s层电子是在3d的内部，
所以Zn应该这样排
label1, occupancy > 4s 2
label1, occupancy > 3d 10
label1, occupancy > 4p 0
label1, occupancy > 4f 0