 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5413 | 回復(fù): 12 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
關(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一類的相關(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ù)此樓» 收錄本帖的淘帖專輯推薦
 彈性計(jì)算 |
» 猜你喜歡
求b區(qū)院校調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
347求調(diào)劑
已經(jīng)有3人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有4人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
-
359求調(diào)劑
已經(jīng)有3人回復(fù)
-
各位老師您好：本人初試372分
已經(jīng)有5人回復(fù)
-
285求調(diào)劑
已經(jīng)有3人回復(fù)
-
282求調(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ù)
- 救命��！VASP計(jì)算磁性出錯(cuò)
已經(jīng)有14人回復(fù)
- 感覺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è)維度放開經(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ù)
awdawdawdd
金蟲 (小有名氣)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時(shí)
- 蟲號(hào): 2237233
- 注冊(cè): 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
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è)人覺得不是很準(zhǔn)確，因?yàn)閿M合直線的話一般多取一些比較好（例如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
awdawdawdd
金蟲 (小有名氣)
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時(shí)
- 蟲號(hào): 2237233
- 注冊(cè): 2013-01-12
- 專業(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的開關(guān)只能選擇默認(rèn)值（NPAR NCORE之類的） 我覺得彈性常數(shù)最好還是手動(dòng)算，尤其是對(duì)于對(duì)稱性不高的情況。通過(guò)施加strain，每一個(gè)類型的strain多弄幾個(gè)點(diǎn)，在ISIF=2下固定cell，relax所有原子，從OUTCAR中獲得stress，然后線性擬合一下就可以獲得彈性常數(shù)。這樣一般都可以獲得比較準(zhǔn)確的值。利用stress-strain法的話，你的相關(guān)參數(shù)都要調(diào)高一點(diǎn)（ENCUT KPOINTS EDIFFG之類的），否則OUTCAR里面的stress不準(zhǔn)確。 或者，也可以利用energy-strain的方法計(jì)算，參數(shù)要求沒有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的第二行分母我覺得應(yīng)該是減號(hào)，否則不volume conserving) |
2樓2015-05-17 18:49:02
3樓2015-05-17 19:08:45
Li小喵
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 91
- 帖子: 39
- 在線: 29.3小時(shí)
- 蟲號(hào): 3480917
- 注冊(cè): 2014-10-17
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2015-05-19 15:15:49
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求b區(qū)院校調(diào)劑 +3 | 周56 2026-03-24 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +6 | 顧九笙要謙虛 2026-03-24 | 6/300 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +4 | 1124361 2026-03-24 | 4/200 |
|
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +15 | 吃吃吃才有意義 2026-03-19 | 16/800 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿國(guó)科過(guò)程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +5 | prayer13 2026-03-23 | 5/250 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 生物學(xué)調(diào)劑 +5 | Surekei 2026-03-21 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研]
320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|