| 5 | 1/1 | 返回列表 |
| 查看: 1580 | 回復(fù): 9 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】關(guān)于測試SIGMA的問題 已有3人參與
|
|||
|
各位,最近我在測試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ā)表 | |
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +7 | 莫xiao 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +15 | 柳溪邊 2026-03-26 | 16/800 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
[考研] 求生物學(xué)調(diào)劑 +7 | 15172915737 2026-04-01 | 7/350 |
|
|
[考研] 省雙一流重點(diǎn)一本大學(xué)招收調(diào)劑 +4 | wwwwffffff 2026-03-31 | 7/350 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +5 | 大饞小子 2026-03-27 | 7/350 |
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
[考研] 求0861交通運(yùn)輸專碩or材料專碩調(diào)劑 +4 | 勒布朗@ 2026-03-31 | 4/200 |
|
|
[考研] 求調(diào)劑:一志愿:南京大學(xué) 專業(yè):0705 總分320 ,本科985,四六級(jí)已過 +3 | lfy760306 2026-03-31 | 3/150 |
|
|
[考研] 352分-085602-一志愿985 +6 | 海納百川Ly 2026-03-29 | 6/300 |
|
|
[考研] 求調(diào)劑 +8 | 11ggg 2026-03-30 | 8/400 |
|
|
[考研] 262求調(diào)劑 +7 | ZZ..000 2026-03-30 | 8/400 |
|
|
[考研] 282求調(diào)劑 +3 | 呼吸都是減肥 2026-03-31 | 3/150 |
|
|
[考研] 本科211總分289,08工學(xué)真心求調(diào)劑 +3 | utopiaE 2026-03-30 | 3/150 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|