| 查看: 2747 | 回复: 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 |
回复此楼» 猜你喜欢
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有22人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
15楼2011-09-14 21:55:54
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
18楼2011-10-08 19:59:39
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
★ ★
zzy870720z(金币+2): 鼓励讨论 2011-10-09 11:30:45
zzy870720z(金币+2): 鼓励讨论 2011-10-09 11:30:45
|
你的赝势截断能设到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 |
19楼2011-10-08 22:11:09
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
20楼2011-10-09 17:25:12
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
23楼2011-11-15 21:52:38