 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1569 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
參數(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ù)此樓» 猜你喜歡
334分 一志愿武理-080500 材料求調(diào)劑
已經(jīng)有10人回復(fù)
-
324分 085600材料與化工
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有11人回復(fù)
-
309求調(diào)劑
已經(jīng)有6人回復(fù)
-
調(diào)劑推薦
已經(jīng)有7人回復(fù)
-
一志愿西交大080500材料學(xué)碩349
已經(jīng)有7人回復(fù)
-
一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有4人回復(fù)
-
352分-085602-一志愿985
已經(jīng)有6人回復(fù)
-
273求調(diào)劑
已經(jīng)有5人回復(fù)
-
349求調(diào)劑
已經(jīng)有4人回復(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)路徑能量圖怎么畫？
已經(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ù)
|
非常感謝你的關(guān)注，U值是參考文獻(xiàn)上的，贗勢(shì)也和文獻(xiàn)一樣，不把Mn當(dāng)成兩類原子的話 也是一樣的情況，用文獻(xiàn)上的encut是沒(méi)問(wèn)題，我就是納悶我怎么就測(cè)不出來(lái)，主要是我這個(gè)能量-截?cái)嗄艿淖兓厔?shì) 沒(méi)有趨于平滑，不知道跟我的磁構(gòu)型有沒(méi)有關(guān)系，文獻(xiàn)是反鐵磁構(gòu)型，我就這么設(shè)置了初始磁矩，或者是我結(jié)構(gòu)優(yōu)化的不好，沒(méi)有進(jìn)入基態(tài)構(gòu)型？但我是分步優(yōu)化，過(guò)程都挺正常的呀。 350 -200.21188 400 -199.69089 450 -199.45863 500 -199.39341 550 -199.41708 600 -199.46676 650 -199.50511 700 -199.54608 |
9樓2013-11-17 18:26:15
2樓2013-11-10 22:53:45
3樓2013-11-11 15:07:28
4樓2013-11-11 20:03:45
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +8 | 素年祭語(yǔ) 2026-03-31 | 11/550 |
|
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +8 | 睿08 2026-03-27 | 8/400 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +17 | Li李魚 2026-03-26 | 17/850 |
|
|
|
[考研] 求收留 +8 | 1943443204 2026-03-28 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +3 | Acesczlo 2026-03-29 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
|
[考研] 293分求調(diào)劑，外語(yǔ)為俄語(yǔ) +5 | 加一一九 2026-03-31 | 5/250 |
|
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
|
[考研] 359求調(diào)劑 +5 | 王了個(gè)楠 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +5 | 兔兔兔111223314 2026-03-29 | 7/350 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|