版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

Vaucanson

金蟲 (小有名氣)

應(yīng)助: 16 (小學(xué)生)
金幣: 812.7
紅花: 2
帖子: 283
在線: 86.6小時
蟲號: 2852894
注冊: 2013-12-06
性別: GG
專業(yè): 金屬結(jié)構(gòu)材料

[求助] 關(guān)于用VASP計算彈性常數(shù)的兩個問題已有2人參與

最近在學(xué)習(xí)用VASP 5.3.5計算彈性常數(shù)，設(shè)置IBRION = 6, ISIF = 4計算。發(fā)現(xiàn)有兩個問題想請教大家：
1. 當(dāng)我用16核計算時，設(shè)置NCORE = 4，走完第一個離子步后報錯，停止計算：

這是否意味著要直接取消掉NCORE一類的相關(guān)參數(shù)？取消之后是warning，提示我設(shè)置NCORE。

2. 由于計算時晶胞有一定的變形，導(dǎo)致K點網(wǎng)格和晶格不匹配，所以往往產(chǎn)生以下的提示，但是會繼續(xù)計算：

VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful... 96

不知道是否對結(jié)果有什么影響？

謝謝大家了。

回復(fù)此樓

1樓 2015-05-17 16:23:23

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

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

Vaucanson

金蟲 (小有名氣)

應(yīng)助: 16 (小學(xué)生)
金幣: 812.7
紅花: 2
帖子: 283
在線: 86.6小時
蟲號: 2852894
注冊: 2013-12-06
性別: GG
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

8樓: Originally posted by vasp001 at 2015-06-14 16:04:43
你好，我最近使用VASP計算D3h晶體的彈性常數(shù)。然后輸出的結(jié)果不是很懂。
不知道如何從VASP的OUTCAR中獲得 Cij, 請多指教。

TOTAL ELASTIC MODULI (kBar)
Direction XX YY ZZ X ...

確定Cij的獨立變量看文獻(xiàn)就行，或者彈塑性力學(xué)上有，這個還是比較容易查找到的。

OUTCAR中給出的是一個二維六階矩陣，XX-XX就是C11，XX-YY是C22，以此類推，不同的是XY-XY對應(yīng)的是C66，我看你在另一個帖子中回復(fù)了，相信你也已經(jīng)知道啦。

QQ什么的還是算啦，因為我自己就是一個菜鳥。

贊一下(1人)

回復(fù)此樓

11樓2015-06-14 20:20:17

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

普通回帖

awdawdawdd

金蟲 (小有名氣)

應(yīng)助: 24 (小學(xué)生)
金幣: 679
散金: 200
紅花: 6
帖子: 111
在線: 67小時
蟲號: 2237233
注冊: 2013-01-12
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Vaucanson: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-17 19:02:39

IBRION=6的時候有關(guān)parallel的開關(guān)只能選擇默認(rèn)值（NPAR NCORE之類的）

我覺得彈性常數(shù)最好還是手動算，尤其是對于對稱性不高的情況。通過施加strain，每一個類型的strain多弄幾個點，在ISIF=2下固定cell，relax所有原子，從OUTCAR中獲得stress，然后線性擬合一下就可以獲得彈性常數(shù)。這樣一般都可以獲得比較準(zhǔn)確的值。利用stress-strain法的話，你的相關(guān)參數(shù)都要調(diào)高一點（ENCUT KPOINTS EDIFFG之類的），否則OUTCAR里面的stress不準(zhǔn)確。

或者，也可以利用energy-strain的方法計算，參數(shù)要求沒有stress-strain那么高，但是有多少個獨立的彈性常數(shù)就要對cell進(jìn)行多少組變形，計算量很大。

ref：
stress-strain法：
http://journals.aps.org/prb/abstract/10.1103/PhysRevB.65.104104
http://scitation.aip.org/content ... 0/10.1063/1.2711762
energy-strain法：
http://scitation.aip.org/content ... /9/10.1063/1.368733 (D9的第二行分母我覺得應(yīng)該是減號，否則不volume conserving)

贊一下

回復(fù)此樓

2樓2015-05-17 18:49:02

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

Vaucanson

金蟲 (小有名氣)

應(yīng)助: 16 (小學(xué)生)
金幣: 812.7
紅花: 2
帖子: 283
在線: 86.6小時
蟲號: 2852894
注冊: 2013-12-06
性別: GG
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

2樓: Originally posted by awdawdawdd at 2015-05-17 18:49:02
IBRION=6的時候有關(guān)parallel的開關(guān)只能選擇默認(rèn)值（NPAR NCORE之類的）

我覺得彈性常數(shù)最好還是手動算，尤其是對于對稱性不高的情況。通過施加strain，每一個類型的strain多弄幾個點，在ISIF=2下固定cell，rela ...

非常感謝你的回答，您說的手動算確實是一個選擇，我現(xiàn)在嘗試的是面心立方的鋁的彈性常數(shù)的計算，像這樣高對稱的結(jié)構(gòu)如果用VASP自帶的功能計算而出現(xiàn)我說的提示不知道是否要緊？因為它會繼續(xù)計算。
另外，不太明白您說的“parallel的開關(guān)只能選擇默認(rèn)值”，好像NCORE不存在默認(rèn)值吧？是NCORE =1的意思？
最后，非常感謝你附的參考文獻(xiàn)。

贊一下

回復(fù)此樓

3樓2015-05-17 19:08:45

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

awdawdawdd

金蟲 (小有名氣)

應(yīng)助: 24 (小學(xué)生)
金幣: 679
散金: 200
紅花: 6
帖子: 111
在線: 67小時
蟲號: 2237233
注冊: 2013-01-12
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Vaucanson: 金幣+20, ★★★很有幫助, 多謝回答 2015-05-18 21:07:14
franch: 金幣+4, 應(yīng)助指數(shù)+1, 謝謝回帖交流。。。。 2015-05-20 11:16:54

引用回帖:

3樓: Originally posted by Vaucanson at 2015-05-17 19:08:45
非常感謝你的回答，您說的手動算確實是一個選擇，我現(xiàn)在嘗試的是面心立方的鋁的彈性常數(shù)的計算，像這樣高對稱的結(jié)構(gòu)如果用VASP自帶的功能計算而出現(xiàn)我說的提示不知道是否要緊？因為它會繼續(xù)計算。
另外，不太明白 ...

有關(guān)parallelisation，主要是關(guān)系到multi-core multi-node的運算的效率
NCORE和NPAR是差不多的，INCAR內(nèi)聲明一個即可。NCORE默認(rèn)值為1。parallelisation只會影響運算效率，一般不影響運算的結(jié)果。

但I(xiàn)BRION=5，6貌似都只能用默認(rèn)值算(也就是你不在INCAR聲明NPAR或者NCORE，或者把這些注釋掉)，然后你用了默認(rèn)值，輸出文件里面肯定會提示你NPAR或者NCORE使用默認(rèn)值計算的效率不高，此時忽略警告即可（因為你只能用默認(rèn)值）。

還有，IBRION=6默認(rèn)應(yīng)該只取兩個點（NFREE=2），加上strain=0的情況總共3個點，個人覺得不是很準(zhǔn)確，因為擬合直線的話一般多取一些比較好（例如NFREE=4），然后調(diào)整好合適的strain的大�。≒OTIM默認(rèn)為0.015Angst，根據(jù)你的體系調(diào)整）

ref：
parallelisation http://cms.mpi.univie.ac.at/vasp ... PLANE_KPAR_tag.html
NCORE http://cms.mpi.univie.ac.at/wiki/index.php/NCORE
IBRION=6、NFREE、POTIM： http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html

贊一下

回復(fù)此樓

4樓2015-05-17 19:24:54

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

Li小喵

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 91
帖子: 39
在線: 29.3小時
蟲號: 3480917
注冊: 2014-10-17
性別: MM
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

你好，我最近也在嘗試著計算彈性常數(shù)，也出現(xiàn)了相同的情況，請問你解決了嗎？我是VASP的初學(xué)者，菜鳥一枚，不知道該如何解決這個問題

贊一下

回復(fù)此樓

5樓2015-05-19 15:15:49

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

huangsunchao

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1437
散金: 144
帖子: 70
在線: 75.1小時
蟲號: 1306932
注冊: 2011-05-26
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Vaucanson: 金幣+10, ★★★很有幫助, 謝謝回帖，我照你的方法試一下 2015-06-13 09:24:28

第一個問題沒有遇見,第二個問題可以通過修改KPOINTS來解決.K點設(shè)置準(zhǔn)則:a1*K1 =a3*K2 =a3*K3 =N,其中N為70或者80.這樣做的目的有兩個:1. 盡可能的均勻劃分K空間,2. 在計算彈性常數(shù)時,晶格體積發(fā)生形變,對應(yīng)K空間劃分也會發(fā)生變化,所以要讓a1*K1 盡可能是10的倍數(shù).

贊一下

回復(fù)此樓

6樓2015-06-12 09:27:32

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

Vaucanson

金蟲 (小有名氣)

應(yīng)助: 16 (小學(xué)生)
金幣: 812.7
紅花: 2
帖子: 283
在線: 86.6小時
蟲號: 2852894
注冊: 2013-12-06
性別: GG
專業(yè): 金屬結(jié)構(gòu)材料

引用回帖:

5樓: Originally posted by Li小喵 at 2015-05-19 15:15:49
你好，我最近也在嘗試著計算彈性常數(shù)，也出現(xiàn)了相同的情況，請問你解決了嗎？我是VASP的初學(xué)者，菜鳥一枚，不知道該如何解決這個問題

第一個問題我覺得不能設(shè)置并行參數(shù)，只能把NCORE注釋掉；第二個暫時還沒有什么好方法，而且最近算的彈性常數(shù)老是錯的，全是負(fù)數(shù)

贊一下

回復(fù)此樓

7樓2015-06-13 09:46:51

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

vasp001

金蟲 (正式寫手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線: 326小時
蟲號: 727265
注冊: 2009-03-20
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

你好，我最近使用VASP計算D3h晶體的彈性常數(shù)。然后輸出的結(jié)果不是很懂。
不知道如何從VASP的OUTCAR中獲得 Cij, 請多指教。

TOTAL ELASTIC MODULI (kBar)
Direction XX       YY       ZZ       XY       YZ       ZX
--------------------------------------------------------------------------------
XX       614.9821 162.6485    -1.1175    -0.5210    0.0000    0.0000
YY       162.6485 615.4451    -1.0843    0.9220    0.0000    0.0000
ZZ       -1.1175    -1.0843    -1.7237    0.0287    0.0000    0.0000
XY       -0.5210    0.9220    0.0287 227.1156    0.0000    0.0000
YZ          0.0000    0.0000    0.0000    0.0000    -0.4616    0.0064
ZX          0.0000    0.0000    0.0000    0.0000    0.0064    -0.4690

另外如何根據(jù)晶體對稱性來確定Cij中獨立變量，
如何從VASP的OUTCAR文件中獲得Cij數(shù)值？

還請不吝指教。

另外你有即時溝通工具嗎？有的話，方便提供一下嗎？

贊一下

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

8樓2015-06-14 16:04:43

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

vasp001

金蟲 (正式寫手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線: 326小時
蟲號: 727265
注冊: 2009-03-20
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by awdawdawdd at 2015-05-17 19:24:54
有關(guān)parallelisation，主要是關(guān)系到multi-core multi-node的運算的效率
NCORE和NPAR是差不多的，INCAR內(nèi)聲明一個即可。NCORE默認(rèn)值為1。parallelisation只會影響運算效率，一般不影響運算的結(jié)果。

但I(xiàn)BRION=5， ...

贊一下

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

9樓2015-06-14 16:04:56

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

vasp001

金蟲 (正式寫手)

應(yīng)助: 69 (初中生)
金幣: 1513.7
散金: 134
紅花: 24
帖子: 702
在線: 326小時
蟲號: 727265
注冊: 2009-03-20
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by Vaucanson at 2015-06-13 09:46:51
第一個問題我覺得不能設(shè)置并行參數(shù)，只能把NCORE注釋掉；第二個暫時還沒有什么好方法，而且最近算的彈性常數(shù)老是錯的，全是負(fù)數(shù)

...

贊一下

回復(fù)此樓

一定要精通MaterialsStudios,VASP.

10樓2015-06-14 16:05:03

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 Vaucanson 的主題更新

返回列表

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +9	Ymlll 2026-03-18	12/600	2026-03-20 00:01 by 23Postgrad
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	3/150	2026-03-19 23:58 by 23Postgrad
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +3	sbdksD 2026-03-19	3/150	2026-03-19 23:21 by fmesaito
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +18	rare12345 2026-03-18	18/900	2026-03-19 22:36 by 學(xué)員8dgXkO
[考研] 085600材料與化工求調(diào)劑 +14	enenenhui 2026-03-13	15/750	2026-03-19 21:53 by allen-yin
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3	一瓶苯甲酸 2026-03-14	3/150	2026-03-19 15:17 by 盡舜堯1
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11	yangfz 2026-03-17	11/550	2026-03-19 15:06 by houyaoxu
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 0817調(diào)劑 +3	沒有答案_ 2026-03-14	3/150	2026-03-19 09:51 by Xu de nuo
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 268求調(diào)劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥于飛Chao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏

24小時熱門版塊排行榜

[求助] 關(guān)于用VASP計算彈性常數(shù)的兩個問題 已有2人參與

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

» 猜你喜歡

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

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

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

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

[求助] 關(guān)于用VASP計算彈性常數(shù)的兩個問題已有2人參與

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