| 10 | 1/1 | 返回列表 |
| 查看: 1550 | 回復(fù): 9 | |||
[交流]
【求助】關(guān)于測試SIGMA的問題 已有3人參與
|
|
各位,最近我在測試Pt(111)面的SIGMA時出現(xiàn)下面的結(jié)果,KPOINTS為9x9x1, INCAR為靜態(tài)計算(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的結(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同樣取上面的一系列值時,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左右,然而二者都沒有達(dá)到這個要求,不知各位有什么看法;對比ISMEAR=1 和0, 發(fā)現(xiàn)同樣的SIGMA值,ISMEAR=1比ISMEAR=0好一些,Pt是金屬體系,但是看到文獻(xiàn)里好像還是用ISMEAR=0的多一些。另外,好像大家的SIGMA值一般去0.1或0.2, 但顯然二者對應(yīng)的entropy值是比較大的。敬請大家賜教,非常感謝! [ Last edited by lgxyz on 2010-10-15 at 09:06 ] |
至尊木蟲 (著名寫手)

|
多謝!我是直接 grep “entropy T*S” 的能量,這個應(yīng)該是總的吧,然后還需要除以整個體系的原子數(shù)以得到ev/atom,是嗎? 如果除以原子個數(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 這樣來看,對于ISMEAR=1來說,SIGMA=0.1就足夠了,而對于ISMEAR=0,SIGMA至少得取到0.03,這個怕是有點太小吧。 請教各位!多謝! |
至尊木蟲 (著名寫手)

至尊木蟲 (著名寫手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +9 | 紙魚ly 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
[基金申請] 山東省面上項目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|