| 查看: 1552 | 回复: 3 | ||
[求助]
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程序本身的原因？ |
回复此楼» 猜你喜欢
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有13人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 用siesta画声子谱的问题
已经有12人回复
- siesta计算dipole moment（电偶极矩）
已经有7人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】siesta 计算DOS应该注意什么问题
已经有3人回复
- 【求助】用siesta计算出现的问题
已经有3人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】关于siesta中赝势价电子轨道的问题
已经有10人回复
|
我的甲烷的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
3楼2011-12-09 14:31:06
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
冻冻2007: 金币+5, ★★★很有帮助 2012-05-04 15:52:07
冻冻2007: 金币+5 2012-05-15 22:28:57
感谢参与，应助指数 +1
冻冻2007: 金币+5, ★★★很有帮助 2012-05-04 15:52:07
冻冻2007: 金币+5 2012-05-15 22:28:57
|
首先，需要声明的是对于采用原子轨道基组计算，随着基组增大，计算结果并不一定要更好，但是趋势是好的，这一点和平面波基组是不同的。 第二，关于你的计算的问题。存在至少两个疑问。一个是你的MeshCutoff 50 Ry我认为是小的，建议使用150Ry甚至200Ry。另外，一个问题是你的基组的问题，建议加入这个选项：PAO.EnergyShift 100 meV，并且使用我的这个值，默认是272 meV。 估计这两个选项可以改善你的问题。 再有就是你的赝势的问题。可以最后考虑它的因素。你先用官方网站上的赝势计算试试。 |
4楼2012-02-28 15:03:10
50