 |
 |
 |
| 24小時熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4398 | 回復: 18 | ||||
babaleo木蟲 (正式寫手)
|
[交流]
【求助】算VASP彈性常數(shù)時的能量-應變出錯 已有12人參與
|
|||
|
用侯老師的腳本練習算彈性常數(shù)時出錯, 以C33為例, (1) 用侯老師文中INCAR, 開始時設ISIF=4 優(yōu)化, 優(yōu)化后的POSCAR拷成OLDPOS (2)加e3=delta, 改defvector.f 中strain(i) 并編譯,以此來算對應的基矢 (3)用侯老師的腳本算應變區(qū)間-0.012 到0.012(步長0.003)的弛豫和能量, 然后擬合能量-應變曲線. 可是得到的能量-應變曲線怎么不是原點對稱的呢?哪一步出錯了?懇請大蝦指點. 圖中紅色是計算結(jié)果,藍色是擬合deltaE/V=3C33^2/2  |
回復此樓» 收錄本帖的淘帖專輯推薦
 彈性常數(shù)【自用】 |
» 猜你喜歡
275求調(diào)劑
已經(jīng)有10人回復
-
286求調(diào)劑
已經(jīng)有6人回復
-
083000學碩274求調(diào)劑
已經(jīng)有5人回復
-
【求調(diào)劑】085601材料工程專碩 | 總分272 |
已經(jīng)有5人回復
-
275求調(diào)劑
已經(jīng)有13人回復
-
0856，材料與化工321分求調(diào)劑
已經(jīng)有8人回復
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有12人回復
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有7人回復
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有8人回復
-
295材料工程專碩求調(diào)劑
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 彈性常數(shù)的計算
已經(jīng)有6人回復
- VASP計算TiO2表面能問題 求助
已經(jīng)有8人回復
- 彈性常數(shù)原子單位轉(zhuǎn)換
已經(jīng)有3人回復
- VASP直接計算彈性常數(shù)出錯
已經(jīng)有6人回復
- vasp計算應力--應變曲線時應力值發(fā)生突變的問題
已經(jīng)有8人回復
- 求助：VASP單點能計算沒有找到能量最低點
已經(jīng)有6人回復
- 計算彈性常數(shù)時的問題
已經(jīng)有10人回復
- 彈性常數(shù)的計算
已經(jīng)有12人回復
- 在VASP4.6中如何利用應力應變關系求彈性常數(shù)
已經(jīng)有4人回復
- VASP 計算E-V曲線求助。。
已經(jīng)有10人回復
- VASP直接計算彈性常數(shù)時報錯：Very bad news! internal error in subroutine IBZKPT:
已經(jīng)有5人回復
- 求助VASP計算彈性常數(shù)的問題
已經(jīng)有5人回復
- VASP5.2 直接計算彈性常數(shù)
已經(jīng)有52人回復
- 求助 利用彈性常數(shù)計算德拜溫度求解
已經(jīng)有55人回復
- 彈性常數(shù)計算出錯
已經(jīng)有5人回復
- 計算不同晶系的彈性常數(shù)
已經(jīng)有260人回復
- 【求助】castep計算彈性常數(shù)問題
已經(jīng)有14人回復
- 【求助】一定壓強下的彈性常數(shù)計算方法
已經(jīng)有25人回復
- 【求助】怎樣用VASP計算反應速率常數(shù)
已經(jīng)有8人回復
- 【討論】用vasp計算過彈性常數(shù)的請進��！更新中��！
已經(jīng)有52人回復
thinkllh
鐵桿木蟲 (正式寫手)
- 應助: 34 (小學生)
- 金幣: 6569.1
- 散金: 50
- 紅花: 1
- 帖子: 538
- 在線: 157小時
- 蟲號: 715800
- 注冊: 2009-03-05
- 性別: GG
- 專業(yè): 極端條件下使用的金屬材料
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
這個是有步驟的： 1. 計算體系E-V關系，每個體積V下先采用ISIF=4，保持體積不變，優(yōu)化晶格形狀和離子位置，ISMEAR取0或1（對金屬），優(yōu)化完之后要增加一步靜態(tài)能量計算，一般ISMEAR=-5. 2. 擬合E-V關系，找到能量最小點的體積（推薦采用Gibbs2程序計算，如果不知道，請自行搜索，本版有我多個相關回復）。 3. 在最小體積下設置ISIF=4，優(yōu)化晶格形狀和離子位置，一遍不放心就兩遍吧，最終的CONTCAR即為最穩(wěn)定構型。 4. 基于上述穩(wěn)定構型，采用應力-應變方法或應力-應變能方法，后者可根據(jù)侯老師教程或直接在VASP中設置IBRION=6計算，采用侯老師的方法是變形之后先優(yōu)化再計算靜態(tài)能，精度更高。 Note： 1. 直接采用ISIF=3得到的CONTCAR是不能用的，尤其是較復雜的體系，不信可以grep volume vasprun.xml，一般最后幾步的體積都在個位數(shù)上震蕩；采用上述方法得到的最小能量對應體積精度一般在小數(shù)點之后，  。2. 凡是想要得到比較精確的穩(wěn)定結(jié)構的，建議都采用上述1-3步逐步施行。 3. 以上是本人失敗了N次之后總結(jié)出來的，這樣得到的應力-應變能曲線對稱性很高，特此分享，以幫助各位少走彎路。 [ Last edited by thinkllh on 2012-6-14 at 09:22 ] |
7樓2012-06-14 09:18:56
sunyang1988
至尊木蟲 (著名寫手)
- 1ST強帖: 1
- 應助: 24 (小學生)
- 貴賓: 0.984
- 金幣: 9397.7
- 散金: 1593
- 紅花: 23
- 沙發(fā): 1
- 帖子: 1744
- 在線: 1255.8小時
- 蟲號: 965316
- 注冊: 2010-03-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構、力學和
3樓2010-08-25 23:38:47
2樓2010-08-25 22:27:05
babaleo
木蟲 (正式寫手)
- 應助: 18 (小學生)
- 金幣: 2391.1
- 散金: 40
- 紅花: 1
- 帖子: 375
- 在線: 477.1小時
- 蟲號: 981394
- 注冊: 2010-03-25
- 專業(yè): 人工晶體
4樓2010-08-26 09:41:20
ufo991
鐵蟲 (正式寫手)
- 1ST強帖: 1
- 應助: 0 (幼兒園)
- 金幣: 1148.3
- 散金: 14
- 紅花: 1
- 帖子: 358
- 在線: 104.6小時
- 蟲號: 671678
- 注冊: 2008-12-10
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構、力學和
5樓2012-06-11 09:11:01
babaleo
木蟲 (正式寫手)
- 應助: 18 (小學生)
- 金幣: 2391.1
- 散金: 40
- 紅花: 1
- 帖子: 375
- 在線: 477.1小時
- 蟲號: 981394
- 注冊: 2010-03-25
- 專業(yè): 人工晶體
6樓2012-06-14 08:14:58
mousekingadv
金蟲 (小有名氣)
- 應助: 51 (初中生)
- 金幣: 1141.8
- 紅花: 1
- 帖子: 140
- 在線: 28.7小時
- 蟲號: 1721002
- 注冊: 2012-03-27
- 專業(yè): 理論和計算化學
8樓2012-06-15 20:22:09
9樓2012-11-14 11:41:44
tonglijia
木蟲 (正式寫手)
- 應助: 12 (小學生)
- 金幣: 1150.6
- 散金: 4
- 紅花: 13
- 帖子: 519
- 在線: 305.9小時
- 蟲號: 1776353
- 注冊: 2012-04-24
- 性別: GG
- 專業(yè): 金屬材料的制備科學與跨學
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
| http://www.gaoyang168.com/bbs/viewthread.php?tid=6328686麻煩您給解答下�。。�！十分感激�。。� |
10樓2013-09-07 22:08:23
相關版塊跳轉(zhuǎn)
我要訂閱樓主 babaleo 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300求調(diào)劑306分 +3 | 26要上岸 2026-03-27 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +4 | jhx777 2026-03-27 | 4/200 |
|
|
|
[考研] 328求調(diào)劑 +6 | 嗯滴的基本都 2026-03-27 | 6/300 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
|
[考研] 321求調(diào)劑 +4 | 璞玉~~ 2026-03-25 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +6 | 紅衣隱官 2026-03-21 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 一志愿北京化工大學材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +6 | 研研，接電話 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿吉大化學322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿河北工業(yè)大學0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|