| 5 | 1/1 | 返回列表 |
| 查看: 1572 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】關(guān)于測(cè)試SIGMA的問題 已有3人參與
|
|||
|
各位,最近我在測(cè)試Pt(111)面的SIGMA時(shí)出現(xiàn)下面的結(jié)果,KPOINTS為9x9x1, INCAR為靜態(tài)計(jì)算(ENCUT=396 PREC=Accurate),請(qǐng)大家指點(diǎn): 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時(shí),entropy T*S的結(jié)果如下: 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 當(dāng)ISMEAR=0 SIGMA同樣取上面的一系列值時(shí),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項(xiàng)在0.001~0.002eV/atom左右,然而二者都沒有達(dá)到這個(gè)要求,不知各位有什么看法;對(duì)比ISMEAR=1 和0, 發(fā)現(xiàn)同樣的SIGMA值,ISMEAR=1比ISMEAR=0好一些,Pt是金屬體系,但是看到文獻(xiàn)里好像還是用ISMEAR=0的多一些。另外,好像大家的SIGMA值一般去0.1或0.2, 但顯然二者對(duì)應(yīng)的entropy值是比較大的。敬請(qǐng)大家賜教,非常感謝! [ Last edited by lgxyz on 2010-10-15 at 09:06 ] |
至尊木蟲 (著名寫手)

|
多謝!我是直接 grep “entropy T*S” 的能量,這個(gè)應(yīng)該是總的吧,然后還需要除以整個(gè)體系的原子數(shù)以得到ev/atom,是嗎? 如果除以原子個(gè)數(shù)36得到一下數(shù)據(jù): 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 這樣來看,對(duì)于ISMEAR=1來說,SIGMA=0.1就足夠了,而對(duì)于ISMEAR=0,SIGMA至少得取到0.03,這個(gè)怕是有點(diǎn)太小吧。 請(qǐng)教各位!多謝! |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)專業(yè)279調(diào)劑 +9 | 哈哈哈^_^ 2026-03-31 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 086000調(diào)劑 +5 | 7901117076 2026-03-26 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +6 | qzenlc 2026-03-29 | 6/300 |
|
|
[考研] 0856 335分 | 封裝or激光加工 老師看看我! +9 | cccchenso 2026-03-29 | 9/450 |
|
|
[考研] 一志愿:西北大學(xué),英一數(shù)一408-284分求調(diào)劑 +6 | 12.27 2026-03-27 | 6/300 |
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +10 | 26要上岸 2026-03-28 | 10/500 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +6 | Ong3 2026-03-27 | 6/300 |
|
|
[考研] 367求調(diào)劑 +7 | 芋泥啵! 2026-03-28 | 7/350 |
|
|
[考研] 生物考研337分求調(diào)劑 +4 | cgxin 2026-03-30 | 6/300 |
|
|
[考研] 322求調(diào)劑 +10 | 舊吢 2026-03-24 | 10/500 |
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考博] 26申博 +3 | 加油沖! 2026-03-26 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|