| 查看: 1431 | 回复: 9 | |||
[交流]
【求助】关于测试SIGMA的问题 已有3人参与
|
|
各位，最近我在测试Pt（111）面的SIGMA时出现下面的结果，KPOINTS为9x9x1, INCAR为静态计算（ENCUT=396 PREC=Accurate），请大家指点： ISMEAR=1 SIGMA=0.01 0.03 0.05 0.07 0.09 0.1 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30时，entropy T*S的结果如下： 0.01 0.000393 0.03 0.002586 0.05 -0.007696 0.07 -0.026601 0.09 -0.040747 0.1 -0.045824 0.12 -0.054866 0.14 -0.065127 0.16 -0.079274 0.18 -0.099725 0.20 -0.127807 0.22 -0.163586 0.24 -0.206818 0.26 -0.255509 0.28 -0.310302 0.30 -0.370317 当ISMEAR=0 SIGMA同样取上面的一系列值时，entropy T*S的值如下： 0.01 -0.003791 0.03 -0.043966 0.05 -0.121563 0.07 -0.223301 0.09 -0.347799 0.1 -0.419281 0.12 -0.581029 0.14 -0.767504 0.16 -0.977665 0.18 -1.210150 0.20 -1.463341 0.22 -1.735587 0.24 -2.025388 0.26 -2.331446 0.28 -2.652650 0.30 -2.988023 mannual上说，SIGMA要尽量小要使得entropy T*S项在0.001~0.002eV/atom左右，然而二者都没有达到这个要求，不知各位有什么看法；对比ISMEAR=1 和0, 发现同样的SIGMA值，ISMEAR=1比ISMEAR=0好一些，Pt是金属体系，但是看到文献里好像还是用ISMEAR=0的多一些。另外，好像大家的SIGMA值一般去0.1或0.2, 但显然二者对应的entropy值是比较大的。敬请大家赐教，非常感谢！ [ Last edited by lgxyz on 2010-10-15 at 09:06 ] |
回复此楼» 猜你喜欢
2025年遐想
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有12人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 金币请教达人一求方差的问题，十分着急，要考试了！
已经有16人回复
- VASP的ISMEAR=0，测试SIGMA的时候，TS为0是怎么回事？
已经有5人回复
- 【EBSD】关于实验测得的轴角对取向差与准确的CSL取向差之间的偏差角度怎么计算？
已经有13人回复
- 水中总氮测试
已经有32人回复
- 【求助】各位大侠！小女子求助钙钛矿结构计算磁性的问题？
已经有36人回复
- 【求助】关于核磁的重水溶剂的问题
已经有3人回复
- 【求助】测试sigma取值时，T*S项的值是负数合理吗？
已经有6人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有51人回复
2楼2010-10-15 10:27:34
3楼2010-10-15 10:34:06
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2010-10-15 10:47:15
|
多谢！我是直接 grep “entropy T*S” 的能量，这个应该是总的吧，然后还需要除以整个体系的原子数以得到ev/atom，是吗？ 如果除以原子个数36得到一下数据： ISMEAR=1 0.01 0.000393 1.09167E-05 0.03 0.002586 7.18333E-05 0.05 -0.007696 -0.000213778 0.07 -0.026601 -0.000738917 0.09 -0.040747 -0.001131861 0.1 -0.045824 -0.001272889 0.12 -0.054866 -0.001524056 0.14 -0.065127 -0.001809083 0.16 -0.079274 -0.002202056 0.18 -0.099725 -0.002770139 0.2 -0.127807 -0.003550194 0.22 -0.163586 -0.004544056 0.24 -0.206818 -0.005744944 0.26 -0.255509 -0.007097472 0.28 -0.310302 -0.0086195 0.3 -0.370317 -0.010286583 ISMEAR=0 0.01 -0.003791 -0.000105306 0.03 -0.043966 -0.001221278 0.05 -0.121563 -0.00337675 0.07 -0.223301 -0.006202806 0.09 -0.347799 -0.009661083 0.1 -0.419281 -0.011646694 0.12 -0.581029 -0.016139694 0.14 -0.767504 -0.021319556 0.16 -0.977665 -0.027157361 0.18 -1.21015 -0.033615278 0.2 -1.463341 -0.040648361 0.22 -1.735587 -0.04821075 0.24 -2.025388 -0.056260778 0.26 -2.331446 -0.064762389 0.28 -2.65265 -0.073684722 0.3 -2.988023 -0.083000639 这样来看，对于ISMEAR=1来说，SIGMA=0.1就足够了，而对于ISMEAR=0，SIGMA至少得取到0.03，这个怕是有点太小吧。 请教各位！多谢！ |
5楼2010-10-15 10:47:50
6楼2010-10-15 10:50:41
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2010-10-15 10:57:41
8楼2010-10-15 11:07:35
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
9楼2010-10-15 11:12:59
10楼2010-10-15 11:46:45