| 10 | 1/1 | 返回列表 |
| 查看: 1551 | 回復(fù): 9 | |||
[交流]
【求助】關(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)教各位!多謝! |
至尊木蟲 (著名寫手)

至尊木蟲 (著名寫手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +4 | Ymlll 2026-03-24 | 4/200 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 材料學(xué)碩297已過四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|