 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 5395 | 回復(fù): 12 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
關(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ù)此樓» 收錄本帖的淘帖專輯推薦
 彈性計算 |
» 猜你喜歡
268求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑
已經(jīng)有4人回復(fù)
-
302求調(diào)劑
已經(jīng)有11人回復(fù)
-
306求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué)
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
313求調(diào)劑
已經(jīng)有3人回復(fù)
-
南昌大學(xué)材料專碩311分求調(diào)劑
已經(jīng)有6人回復(fù)
-
316求調(diào)劑
已經(jīng)有6人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于BCC-Fe的彈性常數(shù)的能量-應(yīng)變計算
已經(jīng)有47人回復(fù)
- 采用VASP計算得到的禁帶寬度為什么比實驗值低
已經(jīng)有11人回復(fù)
- 救命��！VASP計算磁性出錯
已經(jīng)有14人回復(fù)
- 感覺VASP計算缺陷存在很大問題
已經(jīng)有17人回復(fù)
- VASP算DFT+U求助
已經(jīng)有7人回復(fù)
- VASP做結(jié)構(gòu)優(yōu)化時用ISIF=3只想對晶格常數(shù)中的一個維度放開經(jīng)行優(yōu)化
已經(jīng)有7人回復(fù)
- 彈性常數(shù)的計算
已經(jīng)有6人回復(fù)
- VASP計算TiO2表面能問題 求助
已經(jīng)有8人回復(fù)
- vasp計算彈性時的錯誤（使用defvector方法）
已經(jīng)有25人回復(fù)
- vasp自旋極化計算時，初始設(shè)置自旋能否不被優(yōu)化
已經(jīng)有6人回復(fù)
- vasp5.2如何計算DOS和band structure？
已經(jīng)有10人回復(fù)
- 計算彈性常數(shù)時的問題
已經(jīng)有10人回復(fù)
- 彈性常數(shù)的計算
已經(jīng)有12人回復(fù)
- VASP5.2 直接計算彈性常數(shù)
已經(jīng)有52人回復(fù)
- vasp計算中，如何選擇正確的K點
已經(jīng)有7人回復(fù)
- 【求助】vasp 計算 INCAR的設(shè)置
已經(jīng)有9人回復(fù)
- 【求助】一定壓強下的彈性常數(shù)計算方法
已經(jīng)有25人回復(fù)
- 【求助】算VASP彈性常數(shù)時的能量-應(yīng)變出錯
已經(jīng)有18人回復(fù)
- 【討論】用vasp計算過彈性常數(shù)的請進��！更新中！！
已經(jīng)有52人回復(fù)
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ù)值？ 還請不吝指教。 另外你有即時溝通工具嗎？ 有的話，方便提供一下嗎？ |
9樓2015-06-14 16:04:56
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
感謝參與，應(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進行多少組變形，計算量很大。 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) |
2樓2015-05-17 18:49:02
3樓2015-05-17 19:08:45
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
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的運算的效率 NCORE和NPAR是差不多的，INCAR內(nèi)聲明一個即可。NCORE默認(rèn)值為1。parallelisation只會影響運算效率，一般不影響運算的結(jié)果。 但IBRION=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 |
4樓2015-05-17 19:24:54
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線與星海 2026-03-19 | 8/400 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|