24小时热门版块排行榜    

查看: 1270  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

冻冻2007

铜虫 (小有名气)

[求助] siesta 计算电子占据数的问题

我用siesta计算石墨烯纳米带的H和C原子的占据数发现,用SZ基组算出来的和用DZP计算出来的结果有很大的区别,DZP计算的结果很不合理。接着我用水分子做了一下测试,结果如下:
SZ:Species: H
Atom  Qatom  Qorb
               1s
   1  0.840   0.840
   2  0.840   0.840
   3  0.840   0.840
   4  0.840   0.840

Species: C
Atom  Qatom  Qorb
               2s      2py     2pz     2px
   5  4.639   1.216   1.141   1.141   1.141

DZP:
Atom  Qatom  Qorb
               1s      1s      1Ppy    1Ppz    1Ppx
   1  1.203   0.831   0.272   0.034   0.033   0.033
   2  1.203   0.831   0.272   0.033   0.034   0.033
   3  1.204   0.833   0.271   0.033   0.033   0.033
   4  1.204   0.833   0.271   0.033   0.033   0.034
             2s      2s      2py     2pz     2px     2py     2pz     2px
               2Pdxy   2Pdyz   2Pdz2   2Pdxz   2Pdx2-y2
   5     3.187  -0.047   0.750   0.194   0.194   0.192   0.622   0.622   0.622
              0.004   0.008   0.009   0.004   0.011
很明显DZP 的结果不对,是不是siesta程序本身的原因?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

冻冻2007

铜虫 (小有名气)

我的甲烷的input文件为:
SystemName          ch4
SystemLabel         Ch4
NumberOfAtoms      5
NumberOfSpecies     2

MeshCutoff  50 Ry

%block ChemicalSpeciesLabel
1  1  H      # Species index, atomic number, species label
2  6  C
%endblock ChemicalSpeciesLabel

PAO.BasisSize   DZP
XC.functional GGA
XC.authors    PBE
WriteMullikenPop   1

AtomicCoordinatesFormat  Ang
%block AtomicCoordinatesAndAtomicSpecies
-0.380000123   -1.038179245    0.278179290  1
-0.379999607    0.278179338   -1.038179422  1
-0.379999607    0.760000084    0.760000084  1
1.139999981    0.000000000    0.000000000  1
0.000000000    0.000000000    0.000000000  2
%endblock AtomicCoordinatesAndAtomicSpecies

save-rho T
save-delta-rho T
save-total-potential T
save-neutral-atom-potential T
save-hs T
2楼2011-12-08 10:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冻冻2007

铜虫 (小有名气)

有没有人发现这个问题啊?可以用个简单的甲烷测试一下。帮帮忙也
3楼2011-12-09 14:31:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 冻冻2007 的主题更新
信息提示
请填处理意见