 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1567 | 回復: 10 | ||||
[求助]
參數(shù)測試能量怎么都不穩(wěn)定，求救！
|
||||
|
最近在做一個體系的參數(shù)測試，搞了好久了，但是一直沒有得出結果，雖然可以用文獻上的參數(shù)，但是想自己驗證一下，掌握方法。 腳本是這樣 #!/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 結果是這樣： 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點也取過很大的 趨勢也是這樣 sigma 也試了不同的值 能量始終都到不了穩(wěn)定的值 另外我這個POSCAR是用的優(yōu)化好的結構，我這個體系是反鐵磁結構，所以加自旋加Ｕ，但我也試過不加的，能量趨勢也是這樣，然后用文獻的值520 來測K點 已經(jīng)增大到7 15 19了 能量還是一直降低 根本看不到能量變化小與1mev的希望啊…… 各位大神出來指點一二吧，我實在沒招了…… 感激不盡�。。�！  |
回復此樓» 猜你喜歡
262求調劑
已經(jīng)有5人回復
-
0856調劑
已經(jīng)有6人回復
-
張芳銘-中國農(nóng)業(yè)大學-環(huán)境工程專碩-298
已經(jīng)有6人回復
-
考研生物與醫(yī)藥調劑
已經(jīng)有3人回復
-
274求調劑
已經(jīng)有5人回復
-
化學工程085602 305分求調劑
已經(jīng)有26人回復
-
334求調劑
已經(jīng)有3人回復
-
一志愿 南京航空航天大學 ，080500材料科學與工程學碩
已經(jīng)有5人回復
-
看看我！301求調劑 081700
已經(jīng)有3人回復
-
求調劑 生物學 377分
已經(jīng)有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- NEB的計算出現(xiàn)負值，測試的其他能量都比初始的00和終止的10低
已經(jīng)有17人回復
- ENCUT收斂測試時，能量差為什么總不能達到1meV呢？
已經(jīng)有20人回復
- hartree 能量單位和KJ/mol單位換算
已經(jīng)有10人回復
- 求助：北京低碳清潔能源研究所怎么樣？碩士進去待遇如何？
已經(jīng)有19人回復
- 求助：關于肋片效率問題以及換熱系數(shù)測量誤差分析
已經(jīng)有14人回復
- 求助，如何使用Forster能量轉移理論計算藥物小分子與蛋白之間的能量轉移��？
已經(jīng)有11人回復
- 求助castep計算中：晶胞總能量E 與晶胞體積V 的關系圖如何得到，看內容
已經(jīng)有24人回復
- 【求助】用DMOL3怎么算鈦原子能量
已經(jīng)有11人回復
- 【求助】vasp 算出的0k下的能量如何擴展到其他溫度壓力
已經(jīng)有12人回復
- 【求助】反應路徑能量圖怎么畫？
已經(jīng)有7人回復
- 【求助】VASP里的幾個能量怎么區(qū)別呢？？？？？？？？
已經(jīng)有23人回復
- 【求助】請問分子最小能量結果怎么不一致？
已經(jīng)有12人回復
- 【求助】怎么求N在最穩(wěn)定單質形態(tài)下的能量
已經(jīng)有5人回復
- 【求助】能量密度是怎么計算的
已經(jīng)有7人回復
- 【求助】中科院青島能源所怎么樣
已經(jīng)有32人回復
- 【求助】單原子的總能量怎么算？
已經(jīng)有37人回復
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
|
?????????????ж?????? ?????^???? ?]???????? ???O???????? 1) ?z??Y???Л]????e lattice constant and atom position ?????Л]????? key in ?e?` 2) POTCAR ?Л]????e (ex: element ?????e???????e) 3) ?????? VASP ??????????? (?Еr?? vasp , ?o?e lib ??????????~?????????????e?? ??????_?J?F??????VASP?????????????) ???_?J???????? ??? input files ?G???? ?п?????????????? ???????F????F?? (script, POSCAR, POTCAR, KPOINTS) |
» 本帖已獲得的紅花（最新10朵）
kelly6127樓2013-11-14 00:23:03
2樓2013-11-10 22:53:45
3樓2013-11-11 15:07:28
4樓2013-11-11 20:03:45
jpchou
專家顧問 (著名寫手)
-
專家經(jīng)驗: +20 - 1ST強帖: 5
- 應助: 231 (大學生)
- 金幣: 2342
- 散金: 4367
- 紅花: 103
- 帖子: 1061
- 在線: 757.3小時
- 蟲號: 585915
- 注冊: 2008-08-08
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 第一性原理
5樓2013-11-13 17:50:09
6樓2013-11-13 20:54:05
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★
kelly612: 金幣+8 2013-11-17 18:36:29
kelly612: 金幣+8 2013-11-17 18:36:29
| LDAU中的U值應該指定4個就可以了，你指定了5個，大概是你吧自旋上和自旋下的Mn看成兩個不同原子了�？梢园炎孕虾妥孕碌腗n當成同一個原子指定U和POTCAR；然后指定U值。U值不知道你根據(jù)什么選的5.5-1的組合，可以參照文獻指定；截斷能也可以參考文獻值--也就是你找到文獻上用VASP計算的文章，采用和文獻值相同的PAW--也就是電子組態(tài)要相同，然后直接拿來文獻的截斷能用即可，無需測試。另外，能量差小于1meV這個也不是嚴格的，能量-截斷能的變化趨勢趨于平滑即可。。。另外，截斷能太大了也不好，我以前計算發(fā)現(xiàn)過，ECUT太大反而會讓體系收斂緩慢，甚至scf不收斂。 |
8樓2013-11-14 10:24:12
|
非常感謝你的關注，U值是參考文獻上的，贗勢也和文獻一樣，不把Mn當成兩類原子的話 也是一樣的情況，用文獻上的encut是沒問題，我就是納悶我怎么就測不出來，主要是我這個能量-截斷能的變化趨勢 沒有趨于平滑，不知道跟我的磁構型有沒有關系，文獻是反鐵磁構型，我就這么設置了初始磁矩，或者是我結構優(yōu)化的不好，沒有進入基態(tài)構型？但我是分步優(yōu)化，過程都挺正常的呀。 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
10樓2013-11-17 18:29:15
相關版塊跳轉
我要訂閱樓主 kelly612 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學工程085602 305分求調劑 +26 | RichLi_ 2026-03-25 | 26/1300 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿 南京航空航天大學 ，080500材料科學與工程學碩 +5 | @taotao 2026-03-30 | 5/250 |
|
|
[考研] 335求調劑 +3 | 321* 2026-03-31 | 3/150 |
|
|
|
[考研] 083000學碩274求調劑 +11 | Li李魚 2026-03-26 | 11/550 |
|
|
|
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調劑 +12 | 芋泥小鈴鐺 2026-03-28 | 12/600 |
|
|
|
[考研] 283求調劑（080500） +14 | A child 2026-03-27 | 14/700 |
|
|
|
[考研] 材料與化工272求調劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 317求調劑 +10 | 蛋黃咸肉粽 2026-03-26 | 10/500 |
|
|
|
[考研] 299求調劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
|
[考研] 329求調劑 +10 | 鈕恩雪 2026-03-25 | 10/500 |
|
|
|
[考研] 356求調劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 311（085601）求調劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 314求調劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 打過很多競賽，085406控制工程300分，求調劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
|
[考研] 網(wǎng)絡空間安全0839招調劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|