应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于测试SIGMA的问题
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 1591  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lgxyz

铜虫 (正式写手)

[交流] 【求助】关于测试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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-15 09:03:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
不知道你列的结果是不是平均到一个原子后的能量
赞一下(1人) 回复此楼
2楼2010-10-15 10:27:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by quby7503 at 2010-10-15 10:27:34:
不知道你列的结果是不是平均到一个原子后的能量

多谢！我是直接 grep “entropy T*S” 的能量，这个应该是总的吧，然后还需要除以整个体系的原子数以得到ev/atom，是吗？
赞一下 回复此楼
3楼2010-10-15 10:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢交流 2010-10-15 10:49:30
恩，是单个原子达到那个标准。多谢你的数据，ISMEAR=1没用过，原来算金属用ISMEAR=1同样的SIGMA它的entropy要小得多
赞一下(2人) 回复此楼
心之所向无惧无悔
4楼2010-10-15 10:47:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by quby7503 at 2010-10-15 10:27:34:
不知道你列的结果是不是平均到一个原子后的能量

多谢！我是直接 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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2010-10-15 10:47:15:
恩，是单个原子达到那个标准。多谢你的数据，ISMEAR=1没用过，原来算金属用ISMEAR=1同样的SIGMA它的entropy要小得多

恩 应该是这样的吧 manual上也说算金属ISMEAR=1要比ISMEAR=0好一些。
赞一下 回复此楼
6楼2010-10-15 10:50:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ISMEAER=0用的比较多，手册上的说法不管什么计算取ISMEAR=0然后找一个合适的SIGMA就可以。另外ISMEAR=1或2时可能会有负占据，在计算态密度时不能用
赞一下(1人) 回复此楼
心之所向无惧无悔
7楼2010-10-15 10:57:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2010-10-15 10:57:41:
ISMEAER=0用的比较多，手册上的说法不管什么计算取ISMEAR=0然后找一个合适的SIGMA就可以。另外ISMEAR=1或2时可能会有负占据，在计算态密度时不能用

多谢！这个合适的SIGMA(ISMEAR=0)值的标准也是entropy/atom要小于1meV左右吗？那上面的例子岂不是得取SIGMA=0.03，好像有点小吧，你们一般取多少呢?
赞一下 回复此楼
8楼2010-10-15 11:07:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
在保证收敛的前提下，SIGMA越小越好，没有所谓有点小的问题。既然测试了，用测试结果是最可信的
赞一下(2人) 回复此楼
心之所向无惧无悔
9楼2010-10-15 11:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)
引用回帖:
Originally posted by sunyang1988 at 2010-10-15 11:12:59:
在保证收敛的前提下，SIGMA越小越好，没有所谓有点小的问题。既然测试了，用测试结果是最可信的

好的 受教了 多谢！
回复此楼
10楼2010-10-15 11:46:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297分083200求助 +3 aekx 2026-04-05 3/150 2026-04-06 07:57 by 小小树2024
[考研] 求调剂 +8 xzghyuj 2026-04-04 8/400 2026-04-06 06:47 by houyaoxu
[考研] 化学0703-一志愿211-338分求调剂 +8 vants 2026-04-05 8/400 2026-04-06 06:17 by houyaoxu
[考研] 331求调剂 +8 于征yz 2026-04-05 8/400 2026-04-06 00:54 by fmesaito
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 296求调剂 +3 汪！？！ 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 085410人工智能 初试316分 求调剂 +6 残星拂曙 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[有机交流] 甲醇/二氯 1:15过柱子 5+3 a哎y呦w喂 2026-03-31 3/150 2026-04-05 10:42 by 88817753
[考研] 290求调剂085701 +10 1314捧花 2026-04-02 10/500 2026-04-05 10:19 by Sealedwind
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 生物学求调剂 +3 15064154688 2026-04-03 3/150 2026-04-03 10:28 by macy2011
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 一志愿山东大学化学与化工学院材料与化工专硕，360分求调剂 +4 不愿透露姓名的� 2026-04-02 4/200 2026-04-03 09:29 by 遗忘消失的灆
[考研] 调剂 +3 osbbx 2026-04-02 3/150 2026-04-03 07:47 by cc8418
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
信息提示
关闭
请填处理意见
关闭