 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1644 | 回復: 12 | ||||
| 當前主題已經(jīng)存檔。 | ||||
nkleof銀蟲 (正式寫手)
|
[交流]
【求助】晶格優(yōu)化問題 已有7人參與
|
|||
|
我用vasp優(yōu)化一個體系，所得到的結果8.36A遠小于實驗值8.51A。 所采用的贗勢是PAWPBE的，過渡金屬元素Fe和Mn分別試過標準勢和pv勢。采用的優(yōu)化方法分別試過先優(yōu)化原胞體積（ISIF=7），后馳豫離子位置（ISIF=2）；以及全優(yōu)化（ISIF=3）。馳豫方式分別試過CG法（IBRION=2）和準牛頓法（IBRION=1）。為了消除pulay stress，每種方法都以得到的CONTCAR作為POSCAR反復優(yōu)化了3次。 但是以上所有方法所得到的結果都是8.36附近，遠小于實驗值，也小于其它文獻所給出的GGA計算值8.49。 請經(jīng)驗豐富的蟲友幫忙想想可能是在什么地方出問題了？ |
回復此樓» 收錄本帖的淘帖專輯推薦
 VASP |
» 猜你喜歡
求調(diào)劑院校信息
已經(jīng)有4人回復
-
085600材料與化工306
已經(jīng)有4人回復
-
286求調(diào)劑
已經(jīng)有10人回復
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有9人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有4人回復
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有8人回復
-
資源與環(huán)境 調(diào)劑申請(333分)
已經(jīng)有5人回復
-
280求調(diào)劑
已經(jīng)有12人回復
-
269專碩求調(diào)劑
已經(jīng)有5人回復
-
材料學碩301分求調(diào)劑
已經(jīng)有7人回復
gfei2008
木蟲 (小有名氣)
- 應助: 6 (幼兒園)
- 金幣: 1716.7
- 紅花: 1
- 帖子: 176
- 在線: 116.2小時
- 蟲號: 544222
- 注冊: 2008-04-12
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
★ ★ ★
qli0218(金幣+1):謝謝 2010-04-19 21:14
nkleof(金幣+2):感謝建議！ 2010-04-19 21:14
hometownlove(金幣+2):謝謝交流 2010-04-20 15:24
qli0218(金幣+1):謝謝 2010-04-19 21:14
nkleof(金幣+2):感謝建議！ 2010-04-19 21:14
hometownlove(金幣+2):謝謝交流 2010-04-20 15:24
| 你在優(yōu)化之前要進行K-point和ENCUT的測試。得到一個可靠的K點和ENCUT，然后晶格常數(shù)的求法也有相關方面的介紹，就是先求出體系能量與晶格之間的關系然后通過一個方程可以擬合出一條曲線。這個能量最點對應的晶格長度就是晶格常數(shù)。還有，我也用過PAWPBE計算過晶格常數(shù)，是要比試驗值稍微小一點。希望我的回答對你有幫助 |
2樓2010-04-19 21:11:59
nkleof
銀蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時
- 蟲號: 152581
- 注冊: 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
|
我用的k網(wǎng)格疏密程度是k點之間距離0.09 /A，PREC=Accurate因此ENCUT是贗勢內(nèi)默認值的1.3倍，所以這兩個優(yōu)化條件我覺得精度還可以，如果有問題請經(jīng)驗豐富的高人再給點建議。 擬合的優(yōu)化正打算做，但是我覺得vasp的自動優(yōu)化功能誤差應該不會太大，因此用ISIF=2或3的設置算出來的值和靜態(tài)計算擬合優(yōu)化出來的值應該不會有太大差異，是不是有可能我的什么設置出錯？ 如下是我做優(yōu)化用的INCAR PREC = Accurate LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 NELMIN = 4 NSW = 200 LORBIT = 12 EDIFF = 1E-5 IBRION = 2 ISIF = 3 ISPIN = 2 MAGMOM = 8*-6 16*6 32*0 [ Last edited by nkleof on 2010-4-19 at 21:27 ] |
3樓2010-04-19 21:23:24
4樓2010-04-20 10:09:34
calos818
木蟲 (著名寫手)
新人菜鳥
- 應助: 11 (小學生)
- 金幣: 3125.7
- 散金: 77
- 紅花: 3
- 帖子: 1201
- 在線: 730.8小時
- 蟲號: 698507
- 注冊: 2009-02-08
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
5樓2010-04-20 10:26:22
nkleof
銀蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時
- 蟲號: 152581
- 注冊: 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
6樓2010-04-20 11:47:51
7樓2010-04-20 12:25:08
8樓2010-04-20 14:37:49
nkleof
銀蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時
- 蟲號: 152581
- 注冊: 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
9樓2010-04-20 16:05:53
fanchen021
鐵桿木蟲 (正式寫手)
- 應助: 7 (幼兒園)
- 金幣: 7994.9
- 散金: 1844
- 紅花: 4
- 帖子: 776
- 在線: 569小時
- 蟲號: 773188
- 注冊: 2009-05-18
- 專業(yè): 理論和計算化學
10樓2010-04-20 16:31:39
相關版塊跳轉(zhuǎn)
我要訂閱樓主 nkleof 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
|
[考研] 085600材料與化工 +8 | 安全上岸！ 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|