 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5410 | 回復(fù): 12 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
[求助]
關(guān)于用VASP計(jì)算彈性常數(shù)的兩個(gè)問(wèn)題 已有2人參與
|
|||
|
最近在學(xué)習(xí)用VASP 5.3.5計(jì)算彈性常數(shù)，設(shè)置IBRION = 6, ISIF = 4計(jì)算。發(fā)現(xiàn)有兩個(gè)問(wèn)題想請(qǐng)教大家： 1. 當(dāng)我用16核計(jì)算時(shí)，設(shè)置NCORE = 4，走完第一個(gè)離子步后報(bào)錯(cuò)，停止計(jì)算： 這是否意味著要直接取消掉NCORE一類(lèi)的相關(guān)參數(shù)？取消之后是warning，提示我設(shè)置NCORE。 2. 由于計(jì)算時(shí)晶胞有一定的變形，導(dǎo)致K點(diǎn)網(wǎng)格和晶格不匹配，所以往往產(chǎn)生以下的提示，但是會(huì)繼續(xù)計(jì)算： 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 不知道是否對(duì)結(jié)果有什么影響？ 謝謝大家了。 |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 彈性計(jì)算 |
» 猜你喜歡
材料與化工328分調(diào)劑
已經(jīng)有3人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
293求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有9人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
344求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿哈工大，085400，320，求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于BCC-Fe的彈性常數(shù)的能量-應(yīng)變計(jì)算
已經(jīng)有47人回復(fù)
- 采用VASP計(jì)算得到的禁帶寬度為什么比實(shí)驗(yàn)值低
已經(jīng)有11人回復(fù)
- 救命�。ASP計(jì)算磁性出錯(cuò)
已經(jīng)有14人回復(fù)
- 感覺(jué)VASP計(jì)算缺陷存在很大問(wèn)題
已經(jīng)有17人回復(fù)
- VASP算DFT+U求助
已經(jīng)有7人回復(fù)
- VASP做結(jié)構(gòu)優(yōu)化時(shí)用ISIF=3只想對(duì)晶格常數(shù)中的一個(gè)維度放開(kāi)經(jīng)行優(yōu)化
已經(jīng)有7人回復(fù)
- 彈性常數(shù)的計(jì)算
已經(jīng)有6人回復(fù)
- VASP計(jì)算TiO2表面能問(wèn)題 求助
已經(jīng)有8人回復(fù)
- vasp計(jì)算彈性時(shí)的錯(cuò)誤（使用defvector方法）
已經(jīng)有25人回復(fù)
- vasp自旋極化計(jì)算時(shí)，初始設(shè)置自旋能否不被優(yōu)化
已經(jīng)有6人回復(fù)
- vasp5.2如何計(jì)算DOS和band structure？
已經(jīng)有10人回復(fù)
- 計(jì)算彈性常數(shù)時(shí)的問(wèn)題
已經(jīng)有10人回復(fù)
- 彈性常數(shù)的計(jì)算
已經(jīng)有12人回復(fù)
- VASP5.2 直接計(jì)算彈性常數(shù)
已經(jīng)有52人回復(fù)
- vasp計(jì)算中，如何選擇正確的K點(diǎn)
已經(jīng)有7人回復(fù)
- 【求助】vasp 計(jì)算 INCAR的設(shè)置
已經(jīng)有9人回復(fù)
- 【求助】一定壓強(qiáng)下的彈性常數(shù)計(jì)算方法
已經(jīng)有25人回復(fù)
- 【求助】算VASP彈性常數(shù)時(shí)的能量-應(yīng)變出錯(cuò)
已經(jīng)有18人回復(fù)
- 【討論】用vasp計(jì)算過(guò)彈性常數(shù)的請(qǐng)進(jìn)��！更新中��！
已經(jīng)有52人回復(fù)
|
非常感謝你的回答，您說(shuō)的手動(dòng)算確實(shí)是一個(gè)選擇，我現(xiàn)在嘗試的是面心立方的鋁的彈性常數(shù)的計(jì)算，像這樣高對(duì)稱(chēng)的結(jié)構(gòu)如果用VASP自帶的功能計(jì)算而出現(xiàn)我說(shuō)的提示不知道是否要緊？因?yàn)樗鼤?huì)繼續(xù)計(jì)算。 另外，不太明白您說(shuō)的“parallel的開(kāi)關(guān)只能選擇默認(rèn)值”，好像NCORE不存在默認(rèn)值吧？是NCORE =1的意思？ 最后，非常感謝你附的參考文獻(xiàn)。  |
3樓2015-05-17 19:08:45
awdawdawdd
金蟲(chóng) (小有名氣)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時(shí)
- 蟲(chóng)號(hào): 2237233
- 注冊(cè): 2013-01-12
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
Vaucanson: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-17 19:02:39
感謝參與，應(yīng)助指數(shù) +1
Vaucanson: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-17 19:02:39
|
IBRION=6的時(shí)候 有關(guān)parallel的開(kāi)關(guān)只能選擇默認(rèn)值（NPAR NCORE之類(lèi)的） 我覺(jué)得彈性常數(shù)最好還是手動(dòng)算，尤其是對(duì)于對(duì)稱(chēng)性不高的情況。通過(guò)施加strain，每一個(gè)類(lèi)型的strain多弄幾個(gè)點(diǎn)，在ISIF=2下固定cell，relax所有原子，從OUTCAR中獲得stress，然后線性擬合一下就可以獲得彈性常數(shù)。這樣一般都可以獲得比較準(zhǔn)確的值。利用stress-strain法的話(huà)，你的相關(guān)參數(shù)都要調(diào)高一點(diǎn)（ENCUT KPOINTS EDIFFG之類(lèi)的），否則OUTCAR里面的stress不準(zhǔn)確。 或者，也可以利用energy-strain的方法計(jì)算，參數(shù)要求沒(méi)有stress-strain那么高，但是有多少個(gè)獨(dú)立的彈性常數(shù)就要對(duì)cell進(jìn)行多少組變形，計(jì)算量很大。 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的第二行分母我覺(jué)得應(yīng)該是減號(hào)，否則不volume conserving) |
2樓2015-05-17 18:49:02
awdawdawdd
金蟲(chóng) (小有名氣)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時(shí)
- 蟲(chóng)號(hào): 2237233
- 注冊(cè): 2013-01-12
- 專(zhuān)業(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
Vaucanson: 金幣+20, ★★★很有幫助, 多謝回答 2015-05-18 21:07:14
franch: 金幣+4, 應(yīng)助指數(shù)+1, 謝謝回帖交流。。。。 2015-05-20 11:16:54
|
有關(guān)parallelisation，主要是關(guān)系到multi-core multi-node的運(yùn)算的效率 NCORE和NPAR是差不多的，INCAR內(nèi)聲明一個(gè)即可。NCORE默認(rèn)值為1。parallelisation只會(huì)影響運(yùn)算效率，一般不影響運(yùn)算的結(jié)果。 但I(xiàn)BRION=5，6貌似都只能用默認(rèn)值算(也就是你不在INCAR聲明NPAR或者NCORE，或者把這些注釋掉)，然后你用了默認(rèn)值，輸出文件里面肯定會(huì)提示你NPAR或者NCORE使用默認(rèn)值計(jì)算的效率不高，此時(shí)忽略警告即可（因?yàn)槟阒荒苡媚J(rèn)值）。 還有，IBRION=6默認(rèn)應(yīng)該只取兩個(gè)點(diǎn)（NFREE=2），加上strain=0的情況總共3個(gè)點(diǎn)，個(gè)人覺(jué)得不是很準(zhǔn)確，因?yàn)閿M合直線的話(huà)一般多取一些比較好（例如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 |
4樓2015-05-17 19:24:54
Li小喵
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 91
- 帖子: 39
- 在線: 29.3小時(shí)
- 蟲(chóng)號(hào): 3480917
- 注冊(cè): 2014-10-17
- 性別: MM
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2015-05-19 15:15:49
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專(zhuān)碩找調(diào)劑 +3 | 哈哈哈吼吼吼哈 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專(zhuān)碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運(yùn)的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | www_q 2026-03-20 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|