 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1568 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||||
[求助]
參數(shù)測(cè)試能量怎么都不穩(wěn)定，求救！
|
||||
|
最近在做一個(gè)體系的參數(shù)測(cè)試，搞了好久了，但是一直沒(méi)有得出結(jié)果，雖然可以用文獻(xiàn)上的參數(shù)，但是想自己驗(yàn)證一下，掌握方法。 腳本是這樣 #!/bin/sh for i in 350 400 450 500 550 600 650 700 do cat > INCAR << ! SYSTEM = LMP ENCUT = $i ISTART = 0 ; ICHARG = 2 LDAU = .TRUE LDAUTYPE = 1 LDAUL = -1 2 2 -1 -1 LDAUU = 0 5.5 5.5 0 0 LDAUJ = 0 1 1 0 0 ISPIN = 2 MAGMOM = 4*0 1*7 2*-7 1*7 4*0 16*0 ISMEAR = 0; SIGMA = 0.05 NELM = 30 LMIXMIX =4 LREAL = Auto PREC = High ! echo "ENCUT = $i eV" ; mpirun -np 16 vasp5 E=$(grep "energy without entropy" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}') echo $i $E >>comment done KPOINTS 用的 3 6 9 結(jié)果是這樣： 350 -200.25397 400 -199.72309 450 -199.49069 500 -199.42421 550 -199.44941 600 -199.49795 650 -199.53697 700 -199.57734 K點(diǎn)也取過(guò)很大的 趨勢(shì)也是這樣 sigma 也試了不同的值 能量始終都到不了穩(wěn)定的值 另外我這個(gè)POSCAR是用的優(yōu)化好的結(jié)構(gòu)，我這個(gè)體系是反鐵磁結(jié)構(gòu)，所以加自旋加Ｕ，但我也試過(guò)不加的，能量趨勢(shì)也是這樣，然后用文獻(xiàn)的值520 來(lái)測(cè)K點(diǎn) 已經(jīng)增大到7 15 19了 能量還是一直降低 根本看不到能量變化小與1mev的希望啊…… 各位大神出來(lái)指點(diǎn)一二吧，我實(shí)在沒(méi)招了…… 感激不盡�。。�！  |
回復(fù)此樓» 猜你喜歡
070300化學(xué)求調(diào)劑
已經(jīng)有9人回復(fù)
-
085404 22408 315分
已經(jīng)有6人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有22人回復(fù)
-
347求調(diào)劑
已經(jīng)有10人回復(fù)
-
354求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿a區(qū)211，085601-307分求調(diào)劑
已經(jīng)有11人回復(fù)
-
284求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有18人回復(fù)
-
362求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- NEB的計(jì)算出現(xiàn)負(fù)值，測(cè)試的其他能量都比初始的00和終止的10低
已經(jīng)有17人回復(fù)
- ENCUT收斂測(cè)試時(shí)，能量差為什么總不能達(dá)到1meV呢？
已經(jīng)有20人回復(fù)
- hartree 能量單位和KJ/mol單位換算
已經(jīng)有10人回復(fù)
- 求助：北京低碳清潔能源研究所怎么樣？碩士進(jìn)去待遇如何？
已經(jīng)有19人回復(fù)
- 求助：關(guān)于肋片效率問(wèn)題以及換熱系數(shù)測(cè)量誤差分析
已經(jīng)有14人回復(fù)
- 求助，如何使用Forster能量轉(zhuǎn)移理論計(jì)算藥物小分子與蛋白之間的能量轉(zhuǎn)移�。�
已經(jīng)有11人回復(fù)
- 求助castep計(jì)算中：晶胞總能量E 與晶胞體積V 的關(guān)系圖如何得到，看內(nèi)容
已經(jīng)有24人回復(fù)
- 【求助】用DMOL3怎么算鈦原子能量
已經(jīng)有11人回復(fù)
- 【求助】vasp 算出的0k下的能量如何擴(kuò)展到其他溫度壓力
已經(jīng)有12人回復(fù)
- 【求助】反應(yīng)路徑能量圖怎么畫(huà)？
已經(jīng)有7人回復(fù)
- 【求助】VASP里的幾個(gè)能量怎么區(qū)別呢？？？？？？？？
已經(jīng)有23人回復(fù)
- 【求助】請(qǐng)問(wèn)分子最小能量結(jié)果怎么不一致？
已經(jīng)有12人回復(fù)
- 【求助】怎么求N在最穩(wěn)定單質(zhì)形態(tài)下的能量
已經(jīng)有5人回復(fù)
- 【求助】能量密度是怎么計(jì)算的
已經(jīng)有7人回復(fù)
- 【求助】中科院青島能源所怎么樣
已經(jīng)有32人回復(fù)
- 【求助】單原子的總能量怎么算？
已經(jīng)有37人回復(fù)
3樓2013-11-11 15:07:28
2樓2013-11-10 22:53:45
4樓2013-11-11 20:03:45
jpchou
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +20 - 1ST強(qiáng)帖: 5
- 應(yīng)助: 231 (大學(xué)生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線(xiàn): 757.3小時(shí)
- 蟲(chóng)號(hào): 585915
- 注冊(cè): 2008-08-08
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 第一性原理
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
請(qǐng)問(wèn)系統(tǒng)是 bulk or surface or molecule (cluster?) unit cell 多大? 計(jì)算是否收斂? (因?yàn)槟?NELM =30 比預(yù)設(shè)的 40 還要低， 這樣是有危險(xiǎn)的 大多數(shù) VASP 的預(yù)設(shè)值都是有測(cè)過(guò) NELM=40 是個(gè)好數(shù)字 你將它調(diào)低，那就要確保計(jì)算有收斂 請(qǐng)檢查 OSZICAR 如果 electronic loops 是在 30 步以?xún)?nèi) 那就是有收斂 vice versa) |
5樓2013-11-13 17:50:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑：085600材料與化工，考材科基，總分319 +15 | 678lucky 2026-03-31 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 362求調(diào)劑 +7 | 西南交材料專(zhuān)碩3 2026-03-31 | 7/350 |
|
|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +4 | GENJIOW 2026-03-31 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +5 | akie... 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿中國(guó)科學(xué)院大學(xué)265求調(diào)劑 +6 | 恬淡ye 2026-03-31 | 7/350 |
|
|
|
[考研] 材料專(zhuān)碩調(diào)劑 +12 | 椰椰。 2026-03-29 | 12/600 |
|
|
|
[考研] 看看我！301求調(diào)劑 081700 +3 | 一個(gè)小明同學(xué) 2026-03-24 | 3/150 |
|
|
|
[考研] 0710生物學(xué)求調(diào)劑！ +3 | 敘述文 2026-03-31 | 3/150 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +7 | YYYYYNNNNN 2026-03-26 | 9/450 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +7 | JKSOIID 2026-03-26 | 7/350 |
|
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚(yú) 2026-03-27 | 12/600 |
|
|
|
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
|
[考研] 289求調(diào)劑 +16 | 新時(shí)代材料 2026-03-27 | 16/800 |
|
|
|
[考研] 071010 323 分求調(diào)劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|