 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4400 | 回復(fù): 18 | ||||
babaleo木蟲 (正式寫手)
|
[交流]
【求助】算VASP彈性常數(shù)時的能量-應(yīng)變出錯 已有12人參與
|
|||
|
用侯老師的腳本練習(xí)算彈性常數(shù)時出錯, 以C33為例, (1) 用侯老師文中INCAR, 開始時設(shè)ISIF=4 優(yōu)化, 優(yōu)化后的POSCAR拷成OLDPOS (2)加e3=delta, 改defvector.f 中strain(i) 并編譯,以此來算對應(yīng)的基矢 (3)用侯老師的腳本算應(yīng)變區(qū)間-0.012 到0.012(步長0.003)的弛豫和能量, 然后擬合能量-應(yīng)變曲線. 可是得到的能量-應(yīng)變曲線怎么不是原點對稱的呢?哪一步出錯了?懇請大蝦指點. 圖中紅色是計算結(jié)果,藍色是擬合deltaE/V=3C33^2/2  |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 彈性常數(shù)【自用】 |
» 猜你喜歡
11408軟件工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿211 初試270分 求調(diào)劑
已經(jīng)有8人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有6人回復(fù)
-
學(xué)碩274求調(diào)劑
已經(jīng)有9人回復(fù)
-
本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有3人回復(fù)
-
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有7人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
316求調(diào)劑
已經(jīng)有7人回復(fù)
-
322求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 彈性常數(shù)的計算
已經(jīng)有6人回復(fù)
- VASP計算TiO2表面能問題 求助
已經(jīng)有8人回復(fù)
- 彈性常數(shù)原子單位轉(zhuǎn)換
已經(jīng)有3人回復(fù)
- VASP直接計算彈性常數(shù)出錯
已經(jīng)有6人回復(fù)
- vasp計算應(yīng)力--應(yīng)變曲線時應(yīng)力值發(fā)生突變的問題
已經(jīng)有8人回復(fù)
- 求助：VASP單點能計算沒有找到能量最低點
已經(jīng)有6人回復(fù)
- 計算彈性常數(shù)時的問題
已經(jīng)有10人回復(fù)
- 彈性常數(shù)的計算
已經(jīng)有12人回復(fù)
- 在VASP4.6中如何利用應(yīng)力應(yīng)變關(guān)系求彈性常數(shù)
已經(jīng)有4人回復(fù)
- VASP 計算E-V曲線求助。。
已經(jīng)有10人回復(fù)
- VASP直接計算彈性常數(shù)時報錯：Very bad news! internal error in subroutine IBZKPT:
已經(jīng)有5人回復(fù)
- 求助VASP計算彈性常數(shù)的問題
已經(jīng)有5人回復(fù)
- VASP5.2 直接計算彈性常數(shù)
已經(jīng)有52人回復(fù)
- 求助 利用彈性常數(shù)計算德拜溫度求解
已經(jīng)有55人回復(fù)
- 彈性常數(shù)計算出錯
已經(jīng)有5人回復(fù)
- 計算不同晶系的彈性常數(shù)
已經(jīng)有260人回復(fù)
- 【求助】castep計算彈性常數(shù)問題
已經(jīng)有14人回復(fù)
- 【求助】一定壓強下的彈性常數(shù)計算方法
已經(jīng)有25人回復(fù)
- 【求助】怎樣用VASP計算反應(yīng)速率常數(shù)
已經(jīng)有8人回復(fù)
- 【討論】用vasp計算過彈性常數(shù)的請進��！更新中�。�
已經(jīng)有52人回復(fù)
thinkllh
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 34 (小學(xué)生)
- 金幣: 6569.1
- 散金: 50
- 紅花: 1
- 帖子: 538
- 在線: 157小時
- 蟲號: 715800
- 注冊: 2009-03-05
- 性別: GG
- 專業(yè): 極端條件下使用的金屬材料
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
這個是有步驟的： 1. 計算體系E-V關(guān)系，每個體積V下先采用ISIF=4，保持體積不變，優(yōu)化晶格形狀和離子位置，ISMEAR取0或1（對金屬），優(yōu)化完之后要增加一步靜態(tài)能量計算，一般ISMEAR=-5. 2. 擬合E-V關(guān)系，找到能量最小點的體積（推薦采用Gibbs2程序計算，如果不知道，請自行搜索，本版有我多個相關(guān)回復(fù)）。 3. 在最小體積下設(shè)置ISIF=4，優(yōu)化晶格形狀和離子位置，一遍不放心就兩遍吧，最終的CONTCAR即為最穩(wěn)定構(gòu)型。 4. 基于上述穩(wěn)定構(gòu)型，采用應(yīng)力-應(yīng)變方法或應(yīng)力-應(yīng)變能方法，后者可根據(jù)侯老師教程或直接在VASP中設(shè)置IBRION=6計算，采用侯老師的方法是變形之后先優(yōu)化再計算靜態(tài)能，精度更高。 Note： 1. 直接采用ISIF=3得到的CONTCAR是不能用的，尤其是較復(fù)雜的體系，不信可以grep volume vasprun.xml，一般最后幾步的體積都在個位數(shù)上震蕩；采用上述方法得到的最小能量對應(yīng)體積精度一般在小數(shù)點之后，  。2. 凡是想要得到比較精確的穩(wěn)定結(jié)構(gòu)的，建議都采用上述1-3步逐步施行。 3. 以上是本人失敗了N次之后總結(jié)出來的，這樣得到的應(yīng)力-應(yīng)變能曲線對稱性很高，特此分享，以幫助各位少走彎路。 [ Last edited by thinkllh on 2012-6-14 at 09:22 ] |
7樓2012-06-14 09:18:56
sunyang1988
至尊木蟲 (著名寫手)
- 1ST強帖: 1
- 應(yīng)助: 24 (小學(xué)生)
- 貴賓: 0.984
- 金幣: 9397.7
- 散金: 1593
- 紅花: 23
- 沙發(fā): 1
- 帖子: 1744
- 在線: 1255.8小時
- 蟲號: 965316
- 注冊: 2010-03-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
3樓2010-08-25 23:38:47
2樓2010-08-25 22:27:05
babaleo
木蟲 (正式寫手)
- 應(yīng)助: 18 (小學(xué)生)
- 金幣: 2391.1
- 散金: 40
- 紅花: 1
- 帖子: 375
- 在線: 477.1小時
- 蟲號: 981394
- 注冊: 2010-03-25
- 專業(yè): 人工晶體
4樓2010-08-26 09:41:20
ufo991
鐵蟲 (正式寫手)
- 1ST強帖: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1148.3
- 散金: 14
- 紅花: 1
- 帖子: 358
- 在線: 104.6小時
- 蟲號: 671678
- 注冊: 2008-12-10
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2012-06-11 09:11:01
babaleo
木蟲 (正式寫手)
- 應(yīng)助: 18 (小學(xué)生)
- 金幣: 2391.1
- 散金: 40
- 紅花: 1
- 帖子: 375
- 在線: 477.1小時
- 蟲號: 981394
- 注冊: 2010-03-25
- 專業(yè): 人工晶體
6樓2012-06-14 08:14:58
mousekingadv
金蟲 (小有名氣)
- 應(yīng)助: 51 (初中生)
- 金幣: 1141.8
- 紅花: 1
- 帖子: 140
- 在線: 28.7小時
- 蟲號: 1721002
- 注冊: 2012-03-27
- 專業(yè): 理論和計算化學(xué)
8樓2012-06-15 20:22:09
9樓2012-11-14 11:41:44
tonglijia
木蟲 (正式寫手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 1150.6
- 散金: 4
- 紅花: 13
- 帖子: 519
- 在線: 305.9小時
- 蟲號: 1776353
- 注冊: 2012-04-24
- 性別: GG
- 專業(yè): 金屬材料的制備科學(xué)與跨學(xué)
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
| http://www.gaoyang168.com/bbs/viewthread.php?tid=6328686麻煩您給解答下！�。�！十分感激�。�！ |
10樓2013-09-07 22:08:23
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 babaleo 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +3 | 楒桉 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +8 | 辰沐5211314 2026-03-26 | 8/400 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 機械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[有機交流]
 有機合成求助
 20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎？ +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 070300，一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|