 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1943 | 回復(fù): 11 | ||
Gina88木蟲(chóng) (正式寫(xiě)手)
|
[求助]
ISMEAR=1的測(cè)試問(wèn)題
|
|
|
一個(gè)體系有70多個(gè)原子， 我計(jì)算的時(shí)候ISMEAR=1，然后測(cè)試sigma，當(dāng)sigma=0.2的時(shí)候，得到TS為-0.036318 但是幾何優(yōu)化的時(shí)候sigma=0.2時(shí)，得到TS=-0.089262，然后采用sigma=0.15，得到TS=-0.063869 那我sigma應(yīng)該選0.2呢還是應(yīng)該選0.15呢？ 大家有遇到靜態(tài)測(cè)試盒幾何優(yōu)化時(shí)，ISMEAR和sigma取相同的值，卻得到不同數(shù)值TS的情況嗎？ 還有大家測(cè)試ENCUT和Kpoints的時(shí)候，最終精度是不是都達(dá)到了delta-total energy小于1meV呢？我覺(jué)得我自己測(cè)試的時(shí)候好像很難到達(dá)1meV的精度啊，即使到達(dá)1meV的精度，那ENCUT和KPOINTs取值也比較大，因而計(jì)算量很大。但是手冊(cè)上及侯老師的確這么說(shuō)的。 然后查看文獻(xiàn)的時(shí)候，發(fā)現(xiàn)有的人用的delta-total energy的精度是10meV。 所以對(duì)著這個(gè)參數(shù)測(cè)試的精度犯糊涂了， 不知道諸位在具體計(jì)算的時(shí)候都是怎么測(cè)試的，還望不吝賜教！ 先謝謝各位大俠了！ |
回復(fù)此樓» 猜你喜歡
本人考085602 化學(xué)工程 專(zhuān)碩
已經(jīng)有23人回復(fù)
-
求調(diào)劑院校信息
已經(jīng)有4人回復(fù)
-
085600材料與化工306
已經(jīng)有4人回復(fù)
-
286求調(diào)劑
已經(jīng)有10人回復(fù)
-
328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑
已經(jīng)有4人回復(fù)
-
286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有8人回復(fù)
-
資源與環(huán)境 調(diào)劑申請(qǐng)(333分)
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有12人回復(fù)
-
269專(zhuān)碩求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 用VASP-DFPT+phonopy計(jì)算石墨烯聲子譜求助
已經(jīng)有37人回復(fù)
- HSE計(jì)算的一些問(wèn)題，希望大家?guī)臀医鉀Q謝謝
已經(jīng)有13人回復(fù)
- 關(guān)于Si算例結(jié)果分析
已經(jīng)有5人回復(fù)
- 超晶格結(jié)構(gòu)原子弛豫優(yōu)化出問(wèn)題
已經(jīng)有5人回復(fù)
- 求助一下如何超算中心提交vasp任務(wù)
已經(jīng)有8人回復(fù)
- 求大神幫忙解釋截?cái)嗄軠y(cè)試中的問(wèn)題！
已經(jīng)有7人回復(fù)
- sigma測(cè)試問(wèn)題
已經(jīng)有5人回復(fù)
- NEB過(guò)渡態(tài)問(wèn)題，搜索的結(jié)果為什么和插入結(jié)果差距這么大？
已經(jīng)有6人回復(fù)
- 關(guān)于VASP中對(duì)ENCUT收斂性測(cè)試的疑問(wèn)
已經(jīng)有5人回復(fù)
- vasp階段能測(cè)試求助，謝謝�。〖鼻�。
已經(jīng)有15人回復(fù)
- NEB的計(jì)算出現(xiàn)負(fù)值，測(cè)試的其他能量都比初始的00和終止的10低
已經(jīng)有17人回復(fù)
- VASP計(jì)算ZnO的性質(zhì)，為什么ENCUT的測(cè)試曲線是這個(gè)樣子的？
已經(jīng)有10人回復(fù)
- 【重新開(kāi)帖】KPOINTS的測(cè)試腳本問(wèn)題
已經(jīng)有2人回復(fù)
- ENCUT測(cè)試
已經(jīng)有12人回復(fù)
- 求助：帶電缺陷HSE06總能計(jì)算結(jié)果匪夷所思，請(qǐng)賜教！
已經(jīng)有11人回復(fù)
- [求助成功]VASP 計(jì)算 graphene 的錯(cuò)誤
已經(jīng)有15人回復(fù)
- 緊急求助vasp5.2 計(jì)算氧化鈰表面吸附遇到的若干問(wèn)題
已經(jīng)有4人回復(fù)
- 利用GGA+U計(jì)算ZnO帶隙的問(wèn)題
已經(jīng)有3人回復(fù)
- vasp4.6問(wèn)題
已經(jīng)有2人回復(fù)
- 為什么我進(jìn)行ENCUT測(cè)試時(shí)是這個(gè)樣子的？問(wèn)題出在哪里？
已經(jīng)有8人回復(fù)
- 【求助】關(guān)于測(cè)試SIGMA的問(wèn)題
已經(jīng)有9人回復(fù)
wzbhit
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 84 (初中生)
- 金幣: 11959.7
- 散金: 26
- 紅花: 36
- 帖子: 1162
- 在線: 256.7小時(shí)
- 蟲(chóng)號(hào): 1128043
- 注冊(cè): 2010-10-21
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+2, 謝謝提示 2012-11-04 19:47:13
Gina88: 金幣+5, ★★★★★最佳答案 2012-11-04 23:34:25
感謝參與，應(yīng)助指數(shù) +1
liliangfang: 金幣+2, 謝謝提示 2012-11-04 19:47:13
Gina88: 金幣+5, ★★★★★最佳答案 2012-11-04 23:34:25
| 收斂標(biāo)準(zhǔn)是1meV/atom。sigma=0.2時(shí)，得到TS=-0.089262 除以70=-0.001275已經(jīng)在可接受的范圍了。cutoff和kpoint測(cè)試一樣的收斂收斂標(biāo)準(zhǔn)，沒(méi)必要超胞能量變化在1meV，精度太高，因?yàn)槟愕捏w系70atom，屬于比較大的體系。 |
3樓2012-11-04 11:54:45
colinchen2
金蟲(chóng) (小有名氣)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 1652
- 帖子: 59
- 在線: 33.7小時(shí)
- 蟲(chóng)號(hào): 1536228
- 注冊(cè): 2011-12-13
- 性別: GG
- 專(zhuān)業(yè): 半導(dǎo)體物理
2樓2012-11-04 10:21:38
Gina88
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
4樓2012-11-04 23:33:41
Gina88
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
|
您說(shuō)的意思是不一定需要TS小于1meV/atom，只需要量級(jí)在1meV/atom？ 還有，如果按照測(cè)試，sigma要比0.2大好多，如果后面的計(jì)算sigma也設(shè)置得比0.2大好多，得到的TS卻比1meV/atom大好多啊，遇到這種情況該怎么辦呢？ 還有，您說(shuō)cutoff和kpoint測(cè)試一樣的收斂標(biāo)準(zhǔn)，指的是delta-total energy達(dá)到1meV/atom即可認(rèn)為收斂，不必一定達(dá)到手冊(cè)上說(shuō)的delta-total energy達(dá)到1meV？比如：對(duì)于我的體系，delta-total energy只要小于70meV即可認(rèn)為收斂，而不必小于1meV？ |
5樓2012-11-04 23:43:04
colinchen2
金蟲(chóng) (小有名氣)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 1652
- 帖子: 59
- 在線: 33.7小時(shí)
- 蟲(chóng)號(hào): 1536228
- 注冊(cè): 2011-12-13
- 性別: GG
- 專(zhuān)業(yè): 半導(dǎo)體物理
6樓2012-11-05 07:26:43
wzbhit
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 84 (初中生)
- 金幣: 11959.7
- 散金: 26
- 紅花: 36
- 帖子: 1162
- 在線: 256.7小時(shí)
- 蟲(chóng)號(hào): 1128043
- 注冊(cè): 2010-10-21
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
7樓2012-11-05 09:20:31
Gina88
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
|
非常感謝！ 那這么說(shuō)來(lái)，前面測(cè)試就沒(méi)什么用啦？ 還是測(cè)試標(biāo)準(zhǔn)更高，以保證后面的計(jì)算滿足TS在1meV/atom量級(jí)？ 大俠幫我看看我的sigma測(cè)試結(jié)果，怎么不是sigma越大，TS（EENTRO）越大啊？ 0.10 -0.021321 0.12 -0.036190 0.14 -0.048559 0.16 -0.052802 0.18 -0.047653 0.20 -0.036058 0.22 -0.022676 0.24 -0.011552 0.26 -0.005187 |
8樓2012-11-05 09:48:31
wzbhit
至尊木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 84 (初中生)
- 金幣: 11959.7
- 散金: 26
- 紅花: 36
- 帖子: 1162
- 在線: 256.7小時(shí)
- 蟲(chóng)號(hào): 1128043
- 注冊(cè): 2010-10-21
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
9樓2012-11-05 11:44:51
Gina88
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
|
for i in 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 do cat >INCAR< SYSTEM=Si111_relax ISTART=0;ICHARG=2 PREC=Accurate ENCUT=300 GGA=91;VOSKOWN=1 ISMEAR=1;SIGMA=$i EDIFF=1E-6 LREAL=A ALGO=Veryfast;IALGO=48 NGX=170 NGY=170 NGZ=80 NBANDS=200 !ISIF=4 !ISYM=2 !NSW=200 !IBRION=2 !POTIM=0.5 LWAVE=.FALSE. LCHARG=.FALSE. ! |
10樓2012-11-05 21:34:45
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Gina88 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專(zhuān)材355 2026-03-19 | 5/250 |
|
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽(tīng)風(fēng)就是雨； 2026-03-16 | 3/150 |
|