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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】关于测试SIGMA的问题
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 1587  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 一志愿青科085500，初试295分，公共课213分 +3 遇到的人愿望都� 2026-04-05 3/150 2026-04-05 18:45 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666！
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 270分求调剂 +4 maxjxbsk 2026-04-01 4/200 2026-04-05 17:04 by yulian1987
[考研] 085410人工智能 初试316分 求调剂 +6 残星拂曙 2026-03-31 6/300 2026-04-05 12:15 by rainbow11
[考研] 331求调剂 +3 niby 2026-04-02 3/150 2026-04-04 19:56 by 蓝云思雨
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 285求调剂 +6 FZAC123 2026-03-30 6/300 2026-04-03 12:22 by xingguangj
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +5 厦大化工 2026-04-01 5/250 2026-04-02 10:03 by jp9609
[考研] 【求调剂】新能源材料本科，一志愿211，初试321 +6 求调剂学校， 2026-04-02 6/300 2026-04-02 09:41 by 晴空210210
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 318求调剂 +8 七忆77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
信息提示
关闭
请填处理意见
关闭