 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1042 | 回復(fù): 0 | ||
jk射手座新蟲 (著名寫手)
|
[求助]
求教優(yōu)化大晶格結(jié)構(gòu)，優(yōu)化出的能量總是非常大怎么辦
|
|
如題，最近優(yōu)化PKT 1:5~1:7的結(jié)構(gòu)，晶胞比較大，調(diào)試過好幾次參數(shù)經(jīng)常計(jì)算出錯(cuò)，而且算出來的能量經(jīng)常是這樣的形式： free energy TOTEN = -16745.72548078 eV free energy TOTEN = -354385.714304 eV free energy TOTEN = -3343.64830070 eV free energy TOTEN = -271897.42515522 eV free energy TOTEN = -72412.16945720 eV free energy TOTEN = -105104.91766061 eV free energy TOTEN = -223893.13581697 eV 按理說，稍大一些的結(jié)構(gòu)也就偶爾會出現(xiàn)這樣的情況，之前做過1:1~1:4的結(jié)構(gòu)，都沒這樣過，但是現(xiàn)在一下都是這樣，只有個(gè)別會優(yōu)化出正常的結(jié)構(gòu)，比如下邊： free energy TOTEN = -267.661869 eV free energy TOTEN = -112091.39839325 eV free energy TOTEN = -583495.94751343 eV free energy TOTEN = -316312.26362194 eV 所以正常的能量范圍應(yīng)該在-267電子伏附近。。。 我想請教各位前輩這是什么原因造成的（目前通過在網(wǎng)上查找各種原因，我個(gè)人猜測可能服務(wù)器的vasp5.0安裝的時(shí)候就有問題，但是那是兩年前安裝的，當(dāng)時(shí)我根本不在這邊，所以都不知道如果真是那樣要怎么解決這個(gè)問題） 我把我的文件貼出來一下： INCAR SYSTEM = PKT ENCUT=400 PREC=N ISMEAR =0 SIGMA=0.05 EDIFF=1E-3 ISTART=0 ICHARG=2 IBRION=2 ISIF=2 NSW=100 LREAL=Auto（這個(gè)是為了解決計(jì)算中的一個(gè)錯(cuò)誤提示：BRMIX: very serious problems :the old and the new charge density differ） KPOINTS Monkhorst Pack 0 Monkhorst Pack 8 8 4 0 0 0 poscar我寫的是腳本，變換晶格常數(shù)c： #!/bin/bash #for ps in 0.5 10 20 30 50 100 150 200 250 300 350 400 450 500 550 600; do for c in 23.80 23.81 23.82 23.83 23.84 23.85 23.86 23.87 23.88 23.89 23.90 ; do cat > POSCAR <<! data_PKT 1.0 3.90 0.00 0.00 0.00 3.90 0.00 0.00 0.00 $c 1 1 5 5 18 Selective dynamics Direct 0.0 0.0 0.000 F F T 0.5 0.5 0.083 F F T 0.0 0.0 0.167 F F T 0.0 0.0 0.333 F F T 0.0 0.0 0.500 F F T 0.0 0.0 0.667 F F T 0.0 0.0 0.833 F F T 0.5 0.5 0.250 F F T 0.5 0.5 0.417 F F T 0.5 0.5 0.583 F F T 0.5 0.5 0.750 F F T 0.5 0.5 0.917 F F T 0.5 0.0 0.083 F F T 0.0 0.5 0.083 F F T 0.5 0.5 0.000 F F T 0.5 0.0 0.250 F F T 0.0 0.5 0.250 F F T 0.5 0.5 0.167 F F T 0.5 0.0 0.417 F F T 0.0 0.5 0.417 F F T 0.5 0.5 0.333 F F T 0.5 0.0 0.583 F F T 0.0 0.5 0.583 F F T 0.5 0.5 0.500 F F T 0.5 0.0 0.750 F F T 0.0 0.5 0.750 F F T 0.5 0.5 0.677 F F T 0.5 0.0 0.917 F F T 0.0 0.5 0.917 F F T 0.5 0.5 0.833 F F T ! /public/jzh/vasp5 cp CONTCAR CONTCAR$c cp OUTCAR OUTCAR$c grep "TOTEN" OUTCAR | tail -1 >>E.dat done POTCAR用的是GGA型的，之前PAW 和PBE 的（都是vasp自帶的贗勢包）比較容易出錯(cuò)，GGA算的東西比較穩(wěn)定。 跪求大神們提提意見，怎么才能得到結(jié)果差不多可以接受的優(yōu)化結(jié)構(gòu)參數(shù)，不然瞎猜一個(gè)a0 c0 也不是那回事兒呀�。�！小弟先謝謝了！ |
回復(fù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有14人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
加自旋幾何優(yōu)化晶格常數(shù)偏差較大
已經(jīng)有6人回復(fù)
- 優(yōu)化后晶格常數(shù)變化很大原因
已經(jīng)有10人回復(fù)
- VASP做結(jié)構(gòu)優(yōu)化時(shí)用ISIF=3只想對晶格常數(shù)中的一個(gè)維度放開經(jīng)行優(yōu)化
已經(jīng)有7人回復(fù)
- VASP優(yōu)化幾何時(shí)，ISIF=3時(shí)，如何只改變晶格，而不改變位矢方向.
已經(jīng)有11人回復(fù)
- 無法優(yōu)化出文獻(xiàn)中的晶格常數(shù)
已經(jīng)有7人回復(fù)
- 求教各位大俠，優(yōu)化石墨烯的晶格常數(shù)時(shí)，IBRION與ISIF是如何設(shè)置的呢？
已經(jīng)有4人回復(fù)
- 優(yōu)化完成，發(fā)現(xiàn)晶格常數(shù)跟實(shí)驗(yàn)的數(shù)據(jù)相差太大
已經(jīng)有8人回復(fù)
- Castep優(yōu)化結(jié)果的晶格常數(shù)最終結(jié)果怎么看？
已經(jīng)有13人回復(fù)
- castep晶格優(yōu)化問題
已經(jīng)有15人回復(fù)
- 求教Mn的晶格結(jié)構(gòu)，MS的metals文件夾里沒有
已經(jīng)有11人回復(fù)
- 一次型NiAl結(jié)構(gòu)優(yōu)化，晶格矢量變化非常大，求助~
已經(jīng)有5人回復(fù)
- 用LDA贗勢優(yōu)化出來的晶格常數(shù)一定比實(shí)驗(yàn)值小嗎？謝謝！
已經(jīng)有10人回復(fù)
- VASP一次性優(yōu)化銳鈦礦TiO2晶格常數(shù)時(shí)出錯(cuò)！！
已經(jīng)有6人回復(fù)
- 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化
已經(jīng)有17人回復(fù)
- 【求助】vasp優(yōu)化晶格參數(shù)a
已經(jīng)有13人回復(fù)
- 【求助成功】晶格優(yōu)化的結(jié)果很詭異
已經(jīng)有13人回復(fù)
- 【求助成功】晶格優(yōu)化中，murn.f的使用方法
已經(jīng)有7人回復(fù)
- 【求助】侯老師的vasp入門的晶格優(yōu)化部分沒看懂，問一下
已經(jīng)有25人回復(fù)
- 【求助】晶體結(jié)構(gòu)優(yōu)化，計(jì)算出問題，請大家?guī)兔Γ?/a>
已經(jīng)有7人回復(fù)
- 【原創(chuàng)】晶格優(yōu)化中ISIF取值的探究（六方MgO為例）
已經(jīng)有9人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 jk射手座 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +9 | 冬十三 2026-03-15 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷 +8 | 曦熙兮 2026-03-15 | 8/400 |
|
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|