 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2409 | 回復(fù): 12 | ||
[求助]
ENCUT測(cè)試
|
|
體系為280個(gè)原子，選用PAW_PBE勢(shì)。2種原子A,B的 ENMIN ，ENMAX分別為 A B 184:245 300:400 為了測(cè)試ENCUT，之前我首先在ENCUT=700時(shí)測(cè)試了KPOINTS，INCAR如下 SYSTERM ISTART=0 #ICHARG= LWAVE=.False. PREC=high ALGO=Fast LREAL=Auto LELF=.TRUE. #LCHARG=.FALSE. ################electronic relaxation############ NELM=80 ENCUT=700 EDIFF=1E-05 #ISPIN=2 ###################ion relaxtion################ #NSW=400 #IBRION=1 #ISIF=2 #EDIFFG=-0.001 ISMEAR=-5 分別測(cè)試了1x1x1，2x2x2，3x3x3，4x4x4的K點(diǎn)，KPOINTS格式如下： Automatic 0 Gamma x x x 0 0 0 其能量分別為-2290.2356 -2290.7178 -2290.7170 -2290.7176，所以測(cè)試ENCUT時(shí)我采用了2x2x2的KPOINTS。 分別測(cè)試ENCUT=400 500 600 700 800 900 1000，能量如下： 400 -2289.7561 500 -2290.1378 600 -2290.1329 700 -2290.7178 800 -2290.7593 900 -2290.8174 1000 -2290.8880 我的問(wèn)題是為什么都達(dá)到1000了，能量還沒(méi)有收斂。哪里出現(xiàn)了問(wèn)題嗎？ |
回復(fù)此樓» 猜你喜歡
11408軟件工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿211 初試270分 求調(diào)劑
已經(jīng)有8人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有6人回復(fù)
-
學(xué)碩274求調(diào)劑
已經(jīng)有9人回復(fù)
-
本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以
已經(jīng)有7人回復(fù)
-
各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩
已經(jīng)有3人回復(fù)
-
317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑
已經(jīng)有7人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
316求調(diào)劑
已經(jīng)有7人回復(fù)
-
322求調(diào)劑
已經(jīng)有7人回復(fù)
emilyoyang
木蟲(chóng) (正式寫手)
- 應(yīng)助: 232 (大學(xué)生)
- 金幣: 4340.4
- 紅花: 42
- 帖子: 667
- 在線: 189.9小時(shí)
- 蟲(chóng)號(hào): 1630291
- 注冊(cè): 2012-02-20
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2013-01-09 13:43:16
3樓2013-01-09 14:54:23
emilyoyang
木蟲(chóng) (正式寫手)
- 應(yīng)助: 232 (大學(xué)生)
- 金幣: 4340.4
- 紅花: 42
- 帖子: 667
- 在線: 189.9小時(shí)
- 蟲(chóng)號(hào): 1630291
- 注冊(cè): 2012-02-20
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2013-01-09 15:54:10
colinchen2
金蟲(chóng) (小有名氣)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 1652
- 帖子: 59
- 在線: 33.7小時(shí)
- 蟲(chóng)號(hào): 1536228
- 注冊(cè): 2011-12-13
- 性別: GG
- 專業(yè): 半導(dǎo)體物理
5樓2013-01-09 16:00:43
6樓2013-01-09 16:03:03
7樓2013-01-09 16:40:14
emilyoyang
木蟲(chóng) (正式寫手)
- 應(yīng)助: 232 (大學(xué)生)
- 金幣: 4340.4
- 紅花: 42
- 帖子: 667
- 在線: 189.9小時(shí)
- 蟲(chóng)號(hào): 1630291
- 注冊(cè): 2012-02-20
- 性別: MM
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★
xiaoluoj: 金幣+2, ★有幫助 2013-01-09 22:27:46
xiaoluoj: 金幣+2, ★有幫助 2013-01-09 22:27:46
|
可能是你自洽計(jì)算不收斂 雖然你自洽計(jì)算結(jié)束了 但計(jì)算的能量可能不對(duì) 你計(jì)算的體系有沒(méi)有磁性？ 自洽計(jì)算磁矩是否收斂到了優(yōu)化時(shí)的磁矩？（我猜可能沒(méi)有） 所以造成你的能量計(jì)算本身是不準(zhǔn)確的 進(jìn)而造成不同cutoff自洽能不收斂 另外 那么大的體系 文獻(xiàn)上的400是合理的 你用500肯定沒(méi)問(wèn)題 但計(jì)算量可能大了點(diǎn) |
8樓2013-01-09 17:30:34
|
running on 11 nodes distr: one band on 1 nodes, 11 groups vasp.5.2.2 15Apr09 complex POSCAR found type information on POSCAR N Si POSCAR found : 2 types and 280 ions scaLAPACK will be used LDA part: xc-table for Pade appr. of Perdew POSCAR, INCAR and KPOINTS ok, starting setup FFT: planning ...( 1 ) reading WAVECAR entering main loop N E dE d eps ncg rms rms(c) DAV: 1 0.648341706010E+04 0.64834E+04 -0.62749E+05 6292 0.148E+03 DAV: 2 -0.198633704846E+04 -0.84698E+04 -0.82035E+04 7590 0.316E+02 DAV: 3 -0.260881746868E+04 -0.62248E+03 -0.61848E+03 7282 0.898E+01 DAV: 4 -0.262567178229E+04 -0.16854E+02 -0.16800E+02 7876 0.148E+01 DAV: 5 -0.262604758145E+04 -0.37580E+00 -0.37522E+00 7381 0.226E+00 0.103E+02 RMM: 6 -0.229243796515E+04 0.33361E+03 -0.13557E+03 6700 0.332E+01 0.288E+01 RMM: 7 -0.229130484895E+04 0.11331E+01 -0.13174E+02 7342 0.996E+00 0.187E+01 RMM: 8 -0.229118137113E+04 0.12348E+00 -0.16175E+01 7178 0.440E+00 0.578E+00 RMM: 9 -0.229007908022E+04 0.11023E+01 -0.38694E+00 7377 0.183E+00 0.131E+00 RMM: 10 -0.229033387579E+04 -0.25480E+00 -0.28771E+00 7375 0.147E+00 0.247E+00 RMM: 11 -0.229013411261E+04 0.19976E+00 -0.59105E-01 7374 0.741E-01 0.318E-01 RMM: 12 -0.229013715473E+04 -0.30421E-02 -0.43355E-02 7378 0.204E-01 0.451E-02 RMM: 13 -0.229013781928E+04 -0.66455E-03 -0.17970E-03 7418 0.521E-02 0.856E-02 RMM: 14 -0.229013780946E+04 0.98257E-05 -0.24151E-03 7377 0.437E-02 0.276E-02 RMM: 15 -0.229013780307E+04 0.63870E-05 -0.86009E-05 6487 0.860E-03 1 F= -.22901378E+04 E0= -.22901378E+04 d E =0.000000E+00 ----------------------------------------------- 這是ENCUT=500時(shí)，自洽計(jì)算的輸出，INCAR，KPOINTS（2x2x2）如1樓所示。體系沒(méi)有磁性，“雖然你自洽計(jì)算結(jié)束了 但計(jì)算的能量可能不對(duì)”什么意思，為什么自洽計(jì)算結(jié)束了，能量會(huì)不正確？ |
9樓2013-01-09 20:52:11
10樓2013-01-09 20:56:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiaoluoj 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +6 | leaflight 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +7 | 芋泥小鈴鐺 2026-03-28 | 7/350 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +5 | 崔wj 2026-03-26 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +8 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 9/450 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 329求調(diào)劑 +6 | 星野? 2026-03-26 | 6/300 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +4 | 前行必有光 2026-03-28 | 5/250 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛(ài)吃番茄的旭 2026-03-24 | 5/250 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 307求調(diào)劑 +8 | 超級(jí)伊昂大王 2026-03-24 | 9/450 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研]
材料學(xué)碩，求調(diào)劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|