buct2010

金蟲 (正式寫手)

1ST強帖: 10
應(yīng)助: 27 (小學(xué)生)
金幣: 2328.6
散金: 306
紅花: 11
帖子: 351
在線: 231.2小時
蟲號: 1129027
注冊: 2010-10-22
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

[求助] VASP計算時的幾個細節(jié)設(shè)置問題，希望有經(jīng)驗的人過來指點一下！

上手VASP 有一段時間了，參數(shù)一般看文獻取的。但是有兩個基本問題不明白：
VASP計算時，都必須測試 ENCUT 、KPOINTS 和 SIGMA ，使得總能量收斂嗎？還有就是：如果我算原胞時的 KPOINTS 是4*4*6，如果算2*2的超晶胞，KPOINTS 是該設(shè)置為 8*8*12 呢還是保持不變？

[ Last edited by buct2010 on 2011-5-8 at 16:15 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

材料計算模擬實用技巧	VASP	第一性原理計算-tg	計算-vasp
第一性原理計算方法、理論、參數(shù)設(shè)置	科研

» 猜你喜歡

085601材料工程專碩求調(diào)劑已經(jīng)有3人回復(fù)
293求調(diào)劑已經(jīng)有10人回復(fù)
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有6人回復(fù)
085600材料與化工已經(jīng)有3人回復(fù)
材料與化工專碩調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有3人回復(fù)
0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 已經(jīng)有3人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有4人回復(fù)
268求調(diào)劑已經(jīng)有8人回復(fù)
26申博已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

摻雜求高手指點：vasp 不同摻雜濃度模型設(shè)置。已經(jīng)有8人回復(fù)
新手編譯vasp出錯，求高人指點已經(jīng)有3人回復(fù)
vasp自洽計算K點不匹配問題，急��！已經(jīng)有15人回復(fù)
VASP計算：LORBIT設(shè)置出現(xiàn)問題已經(jīng)有5人回復(fù)
vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂已經(jīng)有3人回復(fù)
VASP計算表面吸附非常耗時，請大家指點一下！已經(jīng)有6人回復(fù)
【求助完結(jié)】Linux系統(tǒng)下MS5.5的并行計算調(diào)度問題已經(jīng)有12人回復(fù)
【求助】用vasp跑TS的設(shè)置，大家看看有什么問題嗎？跑了五天了，沒結(jié)束！著急！已經(jīng)有18人回復(fù)
【求助】使用vasp進行帶心凍結(jié)聲子計算的幾個問題~~~ 已經(jīng)有9人回復(fù)
【求助】求VASP計算ZnO能帶和態(tài)密度的詳細輸入設(shè)置已經(jīng)有15人回復(fù)
【求助】新蟲求助！哪位前輩指點一下用vasp做空位缺陷的流程已經(jīng)有7人回復(fù)

1樓 2011-05-08 16:13:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

liuhuimeiwoo

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1469.1
散金: 10
紅花: 3
帖子: 235
在線: 40.2小時
蟲號: 911880
注冊: 2009-11-24
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
youzhizhe(金幣+2): 謝謝交流。 2011-05-30 16:43:24

主要就是測試ECUT和k點一般人家有文獻中會有一些現(xiàn)成的數(shù)值一般來講夠大就行了至于k點那是到空間的 2*2*2的你只要在原來的k點上除以你晶胞的重復(fù)數(shù) 以你的為例就是開始的k點是4*4*6 你生成了一個2*2*2的超晶格所以超胞的k點設(shè)置就應(yīng)該是(4/2)*(4/2)*(6/2) 就是2*2*3

贊一下(3人)

回復(fù)此樓

paper我要好多好的的paper

5樓2011-05-09 10:27:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

dxcharlary

專家顧問 (職業(yè)作家)

專家經(jīng)驗: +83
1ST強帖: 2
應(yīng)助: 251 (大學(xué)生)
金幣: 10417.6
散金: 155
紅花: 65
帖子: 3137
在線: 757.4小時
蟲號: 936718
注冊: 2010-01-04
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

【答案】應(yīng)助回帖

★
buct2010(金幣+2): 很感謝你 2011-05-08 19:35:29
youzhizhe(金幣+1, 1ST強帖+1): 謝謝交流。 2011-05-08 21:07:29

KPOINT是在倒格矢空間中劃分的。。。
如果算2*2*2的超晶胞，KPOINTS 是該設(shè)置為 2*2*3
收斂性測試最好是做，不過一般只要把ecut的數(shù)值設(shè)的大一點都會收斂的。另外ms有一個ultrafine精度，一般的計算，比那個精度高就沒問題。當(dāng)然算DOS聲子譜啥的需要另說。

贊一下(1人)

回復(fù)此樓

2樓2011-05-08 19:21:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

stractor

金蟲 (著名寫手)

1ST強帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
散金: 666
紅花: 30
沙發(fā): 1
帖子: 2224
在線: 426.8小時
蟲號: 394087
注冊: 2007-06-06
專業(yè): 半導(dǎo)體材料

【答案】應(yīng)助回帖

★ ★
buct2010(金幣+3): 謝謝你，知道了。 2011-05-08 19:35:16
youzhizhe(金幣+2, 1ST強帖+1): 謝謝交流。 2011-05-08 21:06:53

"VASP計算時，都必須測試 ENCUT 、KPOINTS 和 SIGMA ，使得總能量收斂嗎？"
主要測試ENCUT和KPOINTS，這兩個必須測試，可以看總能，晶格常數(shù)或者體積模量是否收斂。
“如果我算原胞時的 KPOINTS 是4*4*6，如果算2*2的超晶胞，KPOINTS 是該設(shè)置為 8*8*12 呢還是保持不變？”
采用2*2*2的超胞時，K點選擇2*2*3的采樣方式。

贊一下(1人)

回復(fù)此樓

3樓2011-05-08 19:28:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

buct2010

金蟲 (正式寫手)

1ST強帖: 10
應(yīng)助: 27 (小學(xué)生)
金幣: 2328.6
散金: 306
紅花: 11
帖子: 351
在線: 231.2小時
蟲號: 1129027
注冊: 2010-10-22
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

Originally posted by stractor at 2011-05-08 19:28:57:
"VASP計算時，都必須測試 ENCUT 、KPOINTS 和 SIGMA ，使得總能量收斂嗎？"
主要測試ENCUT和KPOINTS，這兩個必須測試，可以看總能，晶格常數(shù)或者體積模量是否收斂。
“如果我算原胞時的 KPOINTS 是4* ...

總算是弄明白了。之前算的贗勢、ENCUT、KPOINTS都是和文獻（JPCC）一樣的，也不知道對不對

，以后算得自己測試了。謝謝你的解答

贊一下

回復(fù)此樓

4樓2011-05-08 19:44:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

buct2010

金蟲 (正式寫手)

1ST強帖: 10
應(yīng)助: 27 (小學(xué)生)
金幣: 2328.6
散金: 306
紅花: 11
帖子: 351
在線: 231.2小時
蟲號: 1129027
注冊: 2010-10-22
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
管轄: 第一性原理

引用回帖:

Originally posted by liuhuimeiwoo at 2011-05-09 10:27:09:
主要就是測試ECUT和k點一般人家有文獻中會有一些現(xiàn)成的數(shù)值一般來講夠大就行了至于k點那是到空間的 2*2*2的你只要在原來的k點上除以你晶胞的重復(fù)數(shù) 以你的為例就是開始的k點是4*4*6 你生成了一個2*2*2的超 ...

你說的對，因為k空間是倒格式，與原胞的基矢成反比。謝謝你哦

贊一下

回復(fù)此樓

6樓2011-05-09 14:56:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xingyexue

鐵桿木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 5026.4
散金: 362
紅花: 1
帖子: 678
在線: 79.9小時
蟲號: 1165188
注冊: 2010-12-08
專業(yè): 理論和計算化學(xué)

引用回帖:

3樓: Originally posted by stractor at 2011-05-08 19:28:57:
"VASP計算時，都必須測試 ENCUT 、KPOINTS 和 SIGMA ，使得總能量收斂嗎？"
主要測試ENCUT和KPOINTS，這兩個必須測試，可以看總能，晶格常數(shù)或者體積模量是否收斂。
“如果我算原胞時的 KPOINTS 是4* ...

之前看過侯老師的例子發(fā)現(xiàn)例子上的k點都是8 8 8這樣相等的
那么如果k點事像 5 3 2這樣的又該怎么測試k點呢
新手一枚還請不吝賜教~

贊一下

回復(fù)此樓

7樓2012-02-22 10:46:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fanarsenal

木蟲 (正式寫手)

應(yīng)助: 47 (小學(xué)生)
金幣: 2023.2
紅花: 23
帖子: 515
在線: 659.5小時
蟲號: 1069215
注冊: 2010-08-04
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

7樓: Originally posted by xingyexue at 2012-02-22 10:46:57:
之前看過侯老師的例子發(fā)現(xiàn)例子上的k點都是8 8 8這樣相等的
那么如果k點事像 5 3 2這樣的又該怎么測試k點呢
新手一枚還請不吝賜教~

寫個shell腳本，將你想測試的k點組合全部循環(huán)一下，在按照熟練標準測試選擇即可……

贊一下

回復(fù)此樓

8樓2012-02-22 10:54:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xingyexue

鐵桿木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 5026.4
散金: 362
紅花: 1
帖子: 678
在線: 79.9小時
蟲號: 1165188
注冊: 2010-12-08
專業(yè): 理論和計算化學(xué)

引用回帖:

8樓: Originally posted by fanarsenal at 2012-02-22 10:54:19:
寫個shell腳本，將你想測試的k點組合全部循環(huán)一下，在按照熟練標準測試選擇即可……

不是很明白能說的詳細一點么
是將每一個方向的k點分別循環(huán)么
請問你有這樣檢測腳本的例子么？

贊一下

回復(fù)此樓

9樓2012-02-22 14:20:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fanarsenal

木蟲 (正式寫手)

應(yīng)助: 47 (小學(xué)生)
金幣: 2023.2
紅花: 23
帖子: 515
在線: 659.5小時
蟲號: 1069215
注冊: 2010-08-04
專業(yè): 金屬結(jié)構(gòu)材料

你認為要測試的組合全部測試一遍啊，我也沒有什么好辦法，只能一個一個試試，肯定會多出好多，當(dāng)然可以用多個腳本，這樣可以減少重復(fù)，至少根據(jù)supercell的基本信息，abc三個軸的大小，可以排定kakbkc的大小
for ka in
do
for kb in
do
for kc in
do
cat > KPOINTS< A
0
G
$ka $kb $kc
0 0 0
!
done
done
done

贊一下(1人)

回復(fù)此樓

10樓2012-02-23 11:26:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +3	張vvvv 2026-03-15	5/250	2026-03-16 20:25 by 張vvvv
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 材料080500調(diào)劑求收留 +3	一顆meteor 2026-03-13	3/150	2026-03-14 10:54 by peike
[考研] 285 求調(diào)劑資源與環(huán)境一志愿北京化工大學(xué) +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 290求調(diào)劑 +7	ADT 2026-03-12	7/350	2026-03-13 15:17 by JourneyLucky
[考研] 求調(diào)劑資源與環(huán)境 285 +3	未名考生 2026-03-10	3/150	2026-03-13 10:31 by houyaoxu
[考研] 296求調(diào)劑 +3	大口吃飯身體健 2026-03-13	3/150	2026-03-13 10:31 by 學(xué)員8dgXkO
[考研] 調(diào)劑 +5	呵唔哦豁 2026-03-10	5/250	2026-03-10 22:00 by 28375m