buct2010

金蟲 (正式寫手)

1ST強帖: 10
應助: 27 (小學生)
金幣: 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 ]

回復此樓

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

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

» 猜你喜歡

211本，11408一志愿中科院277分，曾在中科院自動化所實習已經(jīng)有7人回復
26申博已經(jīng)有3人回復
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以已經(jīng)有3人回復
材料專碩306英一數(shù)二已經(jīng)有7人回復
材料工程專碩274一志愿211求調(diào)劑已經(jīng)有6人回復
一志愿南京大學，080500材料科學與工程，調(diào)劑已經(jīng)有4人回復
332求調(diào)劑已經(jīng)有5人回復
271求調(diào)劑已經(jīng)有14人回復
275求調(diào)劑已經(jīng)有4人回復
302求調(diào)劑已經(jīng)有8人回復

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

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

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

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

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

liuhuimeiwoo

金蟲 (小有名氣)

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

【答案】應助回帖

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

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

贊一下(3人)

回復此樓

paper我要好多好的的paper

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

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

普通回帖

dxcharlary

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

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

【答案】應助回帖

★
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精度，一般的計算，比那個精度高就沒問題。當然算DOS聲子譜啥的需要另說。

贊一下(1人)

回復此樓

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

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

stractor

金蟲 (著名寫手)

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

【答案】應助回帖

★ ★
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人)

回復此樓

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

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

buct2010

金蟲 (正式寫手)

1ST強帖: 10
應助: 27 (小學生)
金幣: 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）一樣的，也不知道對不對

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

贊一下

回復此樓

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

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

buct2010

金蟲 (正式寫手)

1ST強帖: 10
應助: 27 (小學生)
金幣: 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點上除以你晶胞的重復數(shù) 以你的為例就是開始的k點是4*4*6 你生成了一個2*2*2的超 ...

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

贊一下

回復此樓

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

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

xingyexue

鐵桿木蟲 (正式寫手)

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

引用回帖:

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點呢
新手一枚還請不吝賜教~

贊一下

回復此樓

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

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

fanarsenal

木蟲 (正式寫手)

應助: 47 (小學生)
金幣: 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)一下，在按照熟練標準測試選擇即可……

贊一下

回復此樓

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

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

xingyexue

鐵桿木蟲 (正式寫手)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

fanarsenal

木蟲 (正式寫手)

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

你認為要測試的組合全部測試一遍啊，我也沒有什么好辦法，只能一個一個試試，肯定會多出好多，當然可以用多個腳本，這樣可以減少重復，至少根據(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人)

回復此樓

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

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

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

[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 本人考085602 化學工程專碩 +15	不知道叫什么！ 2026-03-15	17/850	2026-03-17 10:13 by Sammy2
[考研] 304求調(diào)劑 +7	小熊joy 2026-03-14	7/350	2026-03-17 08:53 by 霧散后相遇lc
[考研] 藥學383 求調(diào)劑 +3	藥學chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 085601求調(diào)劑 +3	Du.11 2026-03-16	3/150	2026-03-16 20:42 by 無際的草原
[考研] 070300化學學碩求調(diào)劑 +6	太想進步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 297求調(diào)劑 +4	學海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 327求調(diào)劑 +4	Ffff03 2026-03-10	4/200	2026-03-14 00:17 by JourneyLucky
[考研] 復試調(diào)劑 +9	Copy267 2026-03-10	9/450	2026-03-13 23:45 by userper
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 085600材料與化工 309分請求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 26考研求調(diào)劑 +5	丶宏Sir 2026-03-13	5/250	2026-03-13 13:05 by JourneyLucky
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考博] 26讀博 +4	Rui135246 2026-03-12	10/500	2026-03-13 07:15 by gaobiao
[考研] 調(diào)劑 +5	呵唔哦豁 2026-03-10	5/250	2026-03-10 22:00 by 28375m