 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 5375 | 回復: 12 | |||
[求助]
關于用VASP計算彈性常數(shù)的兩個問題 已有2人參與
|
|||
|
最近在學習用VASP 5.3.5計算彈性常數(shù)，設置IBRION = 6, ISIF = 4計算。發(fā)現(xiàn)有兩個問題想請教大家： 1. 當我用16核計算時，設置NCORE = 4，走完第一個離子步后報錯，停止計算： 這是否意味著要直接取消掉NCORE一類的相關參數(shù)？取消之后是warning，提示我設置NCORE。 2. 由于計算時晶胞有一定的變形，導致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 不知道是否對結果有什么影響？ 謝謝大家了。 |
回復此樓» 收錄本帖的淘帖專輯推薦
 彈性計算 |
» 猜你喜歡
材料專碩英一數(shù)二306
已經(jīng)有6人回復
-
265求調(diào)劑
已經(jīng)有4人回復
-
265求調(diào)劑
已經(jīng)有7人回復
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線
已經(jīng)有8人回復
-
281求調(diào)劑（0805）
已經(jīng)有21人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有20人回復
-
279分求調(diào)劑 一志愿211
已經(jīng)有6人回復
-
085600材料與化工 求調(diào)劑
已經(jīng)有17人回復
-
能源材料化學課題組招收碩士研究生8-10名
已經(jīng)有13人回復
-
317求調(diào)劑
已經(jīng)有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 關于BCC-Fe的彈性常數(shù)的能量-應變計算
已經(jīng)有47人回復
- 采用VASP計算得到的禁帶寬度為什么比實驗值低
已經(jīng)有11人回復
- 救命啊！VASP計算磁性出錯
已經(jīng)有14人回復
- 感覺VASP計算缺陷存在很大問題
已經(jīng)有17人回復
- VASP算DFT+U求助
已經(jīng)有7人回復
- VASP做結構優(yōu)化時用ISIF=3只想對晶格常數(shù)中的一個維度放開經(jīng)行優(yōu)化
已經(jīng)有7人回復
- 彈性常數(shù)的計算
已經(jīng)有6人回復
- VASP計算TiO2表面能問題 求助
已經(jīng)有8人回復
- vasp計算彈性時的錯誤（使用defvector方法）
已經(jīng)有25人回復
- vasp自旋極化計算時，初始設置自旋能否不被優(yōu)化
已經(jīng)有6人回復
- vasp5.2如何計算DOS和band structure？
已經(jīng)有10人回復
- 計算彈性常數(shù)時的問題
已經(jīng)有10人回復
- 彈性常數(shù)的計算
已經(jīng)有12人回復
- VASP5.2 直接計算彈性常數(shù)
已經(jīng)有52人回復
- vasp計算中，如何選擇正確的K點
已經(jīng)有7人回復
- 【求助】vasp 計算 INCAR的設置
已經(jīng)有9人回復
- 【求助】一定壓強下的彈性常數(shù)計算方法
已經(jīng)有25人回復
- 【求助】算VASP彈性常數(shù)時的能量-應變出錯
已經(jīng)有18人回復
- 【討論】用vasp計算過彈性常數(shù)的請進��！更新中！！
已經(jīng)有52人回復
11樓2015-06-14 20:20:17
awdawdawdd
金蟲 (小有名氣)
- 應助: 24 (小學生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時
- 蟲號: 2237233
- 注冊: 2013-01-12
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
Vaucanson: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-17 19:02:39
感謝參與，應助指數(shù) +1
Vaucanson: 金幣+20, ★★★很有幫助, 非常感謝 2015-05-17 19:02:39
|
IBRION=6的時候 有關parallel的開關只能選擇默認值（NPAR NCORE之類的） 我覺得彈性常數(shù)最好還是手動算，尤其是對于對稱性不高的情況。通過施加strain，每一個類型的strain多弄幾個點，在ISIF=2下固定cell，relax所有原子，從OUTCAR中獲得stress，然后線性擬合一下就可以獲得彈性常數(shù)。這樣一般都可以獲得比較準確的值。利用stress-strain法的話，你的相關參數(shù)都要調(diào)高一點（ENCUT KPOINTS EDIFFG之類的），否則OUTCAR里面的stress不準確。 或者，也可以利用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的第二行分母我覺得應該是減號，否則不volume conserving) |
2樓2015-05-17 18:49:02
3樓2015-05-17 19:08:45
awdawdawdd
金蟲 (小有名氣)
- 應助: 24 (小學生)
- 金幣: 679
- 散金: 200
- 紅花: 6
- 帖子: 111
- 在線: 67小時
- 蟲號: 2237233
- 注冊: 2013-01-12
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Vaucanson: 金幣+20, ★★★很有幫助, 多謝回答 2015-05-18 21:07:14
franch: 金幣+4, 應助指數(shù)+1, 謝謝回帖交流。。。。 2015-05-20 11:16:54
Vaucanson: 金幣+20, ★★★很有幫助, 多謝回答 2015-05-18 21:07:14
franch: 金幣+4, 應助指數(shù)+1, 謝謝回帖交流。。。。 2015-05-20 11:16:54
|
有關parallelisation，主要是關系到multi-core multi-node的運算的效率 NCORE和NPAR是差不多的，INCAR內(nèi)聲明一個即可。NCORE默認值為1。parallelisation只會影響運算效率，一般不影響運算的結果。 但IBRION=5，6貌似都只能用默認值算(也就是你不在INCAR聲明NPAR或者NCORE，或者把這些注釋掉)，然后你用了默認值，輸出文件里面肯定會提示你NPAR或者NCORE使用默認值計算的效率不高，此時忽略警告即可（因為你只能用默認值）。 還有，IBRION=6默認應該只取兩個點（NFREE=2），加上strain=0的情況總共3個點，個人覺得不是很準確，因為擬合直線的話一般多取一些比較好（例如NFREE=4），然后調(diào)整好合適的strain的大�。≒OTIM默認為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小喵
新蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 91
- 帖子: 39
- 在線: 29.3小時
- 蟲號: 3480917
- 注冊: 2014-10-17
- 性別: MM
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
5樓2015-05-19 15:15:49
huangsunchao
木蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 1437
- 散金: 144
- 帖子: 70
- 在線: 75.1小時
- 蟲號: 1306932
- 注冊: 2011-05-26
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
6樓2015-06-12 09:27:32
7樓2015-06-13 09:46:51
vasp001
金蟲 (正式寫手)
- 應助: 69 (初中生)
- 金幣: 1513.7
- 散金: 134
- 紅花: 24
- 帖子: 702
- 在線: 326小時
- 蟲號: 727265
- 注冊: 2009-03-20
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
你好，我最近使用VASP計算D3h晶體的彈性常數(shù)。然后輸出的結果不是很懂。 不知道如何從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ù)值？ 還請不吝指教。 另外你有即時溝通工具嗎？ 有的話，方便提供一下嗎？ |
8樓2015-06-14 16:04:43
vasp001
金蟲 (正式寫手)
- 應助: 69 (初中生)
- 金幣: 1513.7
- 散金: 134
- 紅花: 24
- 帖子: 702
- 在線: 326小時
- 蟲號: 727265
- 注冊: 2009-03-20
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
你好，我最近使用VASP計算D3h晶體的彈性常數(shù)。然后輸出的結果不是很懂。 不知道如何從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
vasp001
金蟲 (正式寫手)
- 應助: 69 (初中生)
- 金幣: 1513.7
- 散金: 134
- 紅花: 24
- 帖子: 702
- 在線: 326小時
- 蟲號: 727265
- 注冊: 2009-03-20
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
你好，我最近使用VASP計算D3h晶體的彈性常數(shù)。然后輸出的結果不是很懂。 不知道如何從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ù)值？ 還請不吝指教。 另外你有即時溝通工具嗎？ 有的話，方便提供一下嗎？ |
10樓2015-06-14 16:05:03
相關版塊跳轉
我要訂閱樓主 Vaucanson 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研]
 271材料工程求調(diào)劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 一志愿985，本科211，0817化學工程與技術319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
|
[考研]
0854，計算機類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研]
11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 080500，材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|