| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 2410 | 回復(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)題嗎? |
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
金蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫(xiě)手)
|
可能是你自洽計(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) |
|
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ì)不正確? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +5 | 楒桉 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 320分,材料與化工專(zhuān)業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +7 | 麻辣魷魚(yú) 2026-03-27 | 8/400 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過(guò) 2026-03-27 | 6/300 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|
|
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|