 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 491 | 回復(fù): 0 | ||
qphll金蟲 (正式寫手)
|
[求助]
收斂測(cè)試, LaFeO3, 能量變化異常
|
|
嘗試用vasp計(jì)算LaFeO3, 在收斂測(cè)試的時(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測(cè)試, 最后得到的圖中能量是線性的: 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測(cè)試能量圖: http://edu.www.gaoyang168.com/attachment/87/87/89654_1333900718_282.jpg 這是怎么回事呢? 可能的問題: a. INCAR中的參數(shù)設(shè)置?? b.贗勢(shì)的選擇?? 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得到. 我的打算是做完測(cè)試以后, 再做晶體的結(jié)構(gòu)優(yōu)化. 大家有什么建議? 多謝. ENCUT.test SIGMA.test [ Last edited by qphll on 2012-4-9 at 01:03 ] |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(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ì)算分子軌道卻無(wú)法顯示軌道能量
已經(jīng)有5人回復(fù)
- 關(guān)于dmol3的收斂性測(cè)試
已經(jīng)有12人回復(fù)
- 天然氣燃燒的能量利用效率應(yīng)該怎么計(jì)算？
已經(jīng)有4人回復(fù)
- 請(qǐng)教扣電測(cè)試異常
已經(jīng)有4人回復(fù)
- 兩個(gè)原子靠近成鍵的過(guò)程中，系統(tǒng)總能量會(huì)不會(huì)發(fā)生變化？
已經(jīng)有3人回復(fù)
- XPS中能量如何校準(zhǔn)
已經(jīng)有7人回復(fù)
- 計(jì)算表面時(shí)：Vacuum真空層和slab層厚度收斂性測(cè)試，保持超晶胞體積不變嗎？？
已經(jīng)有10人回復(fù)
- 【其他】關(guān)于收斂測(cè)試的問題
已經(jīng)有22人回復(fù)
- 【求助】收斂測(cè)試怎么做？？？
已經(jīng)有17人回復(fù)
- 【求助】收斂性測(cè)試的問題
已經(jīng)有40人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qphll 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 梁梁校校 2026-03-19 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 329求調(diào)劑 +5 | 想上學(xué)吖吖 2026-03-19 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|