| 查看: 3208 | 回复: 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 |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
改王之王
金虫 (正式写手)
- 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 ：电子结构
2楼2011-09-13 20:50:36
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
3楼2011-09-13 21:08:33
改王之王
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1120.7
- 散金: 61
- 红花: 3
- 帖子: 817
- 在线: 113.5小时
- 虫号: 968169
- 注册: 2010-03-11
- 专业: 凝聚态物性 II ：电子结构
4楼2011-09-13 21:29:19
25