 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 495 | 回復(fù): 0 | |||
qphll金蟲 (正式寫手)
|
[求助]
收斂測試, LaFeO3, 能量變化異常
|
|
嘗試用vasp計(jì)算LaFeO3, 在收斂測試的時(shí)候發(fā)現(xiàn)隨著encut的增加, 總能反而升高. (1) INCAR腳本 m WAVECAR for i in 200 250 300 350 400 450 500 550 600 do cat > INCAR << ! SYSTEM = LF orthorhombic ENCUT = $i ISMEAR = 0 ! Gaussian smearing SIGMA = 0.1 LREAL= Auto PREC= Accurate ! echo "ENCUT= $i eV" ; time; ./vasp5211 ; E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` echo $i $E >> compare done (2) KPOINTS 文件 Automatic generation 0 Monhkorst-Pack 4 4 4 0.0 0.0 0.0 (3) POSCAR 文件 rthorhombic 1.0 5.5349998474 0.0000000000 0.0000000000 0.0000000000 7.8880000114 0.0000000000 0.0000000000 0.0000000000 5.5989999771 12 4 4 ! O Fe La Direct 0.270900011 0.034699999 0.229399994 0.729099989 0.965300024 0.770600021 0.229099989 0.965300024 0.729399979 0.770900011 0.034699999 0.270600021 0.729099989 0.534699976 0.770600021 0.270900011 0.465299994 0.229399994 0.770900011 0.465299994 0.270600021 0.229099989 0.534699976 0.729399979 0.493999988 0.250000000 0.564400017 0.506000042 0.750000000 0.435599983 0.006000012 0.750000000 0.064400017 0.993999958 0.250000000 0.935599983 0.000000000 0.000000000 0.000000000 0.500000000 0.000000000 0.500000000 0.000000000 0.500000000 0.000000000 0.500000000 0.500000000 0.500000000 0.518400013 0.250000000 0.008900000 0.481599987 0.750000000 0.991100013 0.981599987 0.750000000 0.508899987 0.018400013 0.250000000 0.491100013 最后得到的圖是這樣的: http://edu.www.gaoyang168.com/attachment/16/b3/89654_1333900508_628.jpg 同樣, 如果修改腳本, 對(duì)sigma測試, 最后得到的圖中能量是線性的: SIGMA 腳本 rm WAVECAR for i in 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 do cat > INCAR << ! SYSTEM = LF orthorhombic ENCUT = 500 ISMEAR = 0 ! Gaussian smearing SIGMA = $i LREAL= Auto PREC= Accurate ISTART = 0 ICHARGE = 2 ! echo "SIGMA= $i eV" ; time; ./vasp5211 ; E=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` echo $i $E >> sigma.compare done sigma測試能量圖: http://edu.www.gaoyang168.com/attachment/87/87/89654_1333900718_282.jpg 這是怎么回事呢? 可能的問題: a. INCAR中的參數(shù)設(shè)置?? b.贗勢的選擇?? PAW-PBE的O, La, Fe有幾個(gè), 我倒沒有進(jìn)行選擇 c. 先優(yōu)化結(jié)構(gòu), 再單點(diǎn)能計(jì)算?? POSCAR直接從文獻(xiàn)報(bào)導(dǎo)的CIF得到. 我的打算是做完測試以后, 再做晶體的結(jié)構(gòu)優(yōu)化. 大家有什么建議? 多謝. ENCUT.test SIGMA.test [ Last edited by qphll on 2012-4-9 at 01:03 ] |
回復(fù)此樓» 猜你喜歡
材料292調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑，求導(dǎo)師收
已經(jīng)有3人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有6人回復(fù)
-
305分求調(diào)劑（食品工程）
已經(jīng)有7人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有4人回復(fù)
-
收08調(diào)劑生
已經(jīng)有3人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
求助關(guān)于熒光共振能量轉(zhuǎn)移的問題？
已經(jīng)有8人回復(fù)
- 用Dmol優(yōu)化收斂的結(jié)構(gòu)，用vasp不收斂？
已經(jīng)有9人回復(fù)
- 3.22文獻(xiàn)解讀：三元材料在高能量密度鋰離子電池中的應(yīng)用
已經(jīng)有236人回復(fù)
- 干涉光能量哪里去了？
已經(jīng)有5人回復(fù)
- 求助：castep計(jì)算參數(shù)設(shè)置沒有收斂
已經(jīng)有9人回復(fù)
- 用CHEM3D軟件PM3理論水平計(jì)算分子軌道卻無法顯示軌道能量
已經(jīng)有5人回復(fù)
- 關(guān)于dmol3的收斂性測試
已經(jīng)有12人回復(fù)
- 天然氣燃燒的能量利用效率應(yīng)該怎么計(jì)算？
已經(jīng)有4人回復(fù)
- 請(qǐng)教扣電測試異常
已經(jīng)有4人回復(fù)
- 兩個(gè)原子靠近成鍵的過程中，系統(tǒng)總能量會(huì)不會(huì)發(fā)生變化？
已經(jīng)有3人回復(fù)
- XPS中能量如何校準(zhǔn)
已經(jīng)有7人回復(fù)
- 計(jì)算表面時(shí)：Vacuum真空層和slab層厚度收斂性測試，保持超晶胞體積不變嗎？？
已經(jīng)有10人回復(fù)
- 【其他】關(guān)于收斂測試的問題
已經(jīng)有22人回復(fù)
- 【求助】收斂測試怎么做？？？
已經(jīng)有17人回復(fù)
- 【求助】收斂性測試的問題
已經(jīng)有40人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qphll 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +5 | niko- 2026-03-22 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研]
 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
|
[考研] 22 350 本科985求調(diào)劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研]
材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研]
295復(fù)試調(diào)劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|