【文章內(nèi)容簡介】
groundstate, you can try the following setting: LMAXMIX = 4 for d elements LMAXMIX = 6 for f elements 這個(gè) FLAG對(duì)于含 d電子和 f電子的體系是非常重要的,很大一部分體系的收斂問題可以通過設(shè)置合適的 LMAXMIX值來解決。 ALGO, IALGO, LDIAG If the selfconsistency loop does not converge within 40 steps, it will probably not converge at all. In this case you should reconsider the tags IALGO, LDIAG, and the mixingparameters. 一般情況下,或使用 IALGO=48時(shí)遇到收斂問題的話,可以考慮設(shè) IALGO為 38,或設(shè)置 ALGO=Normal or Fast (in VAS and later versions)。 Default ALGO = Normal NELMDL NELMDL gives the number of nonselfconsistent steps at the beginning In some cases (for instance MD’ s, or ionic relaxation) you might set NELMIN to a larger value (4 to 8) mixingparameters 盡管 VASP說明書中給出了調(diào)節(jié) AMIX和 BMIX的一些較為明確的建議,但是實(shí)際去調(diào)節(jié)的時(shí)候,還是挺難的,但原則上說,是可以通過調(diào)節(jié)這兩個(gè) Flag來使得收斂問題得以解決的,只是得有耐心。 收斂問題還跟 kmesh及 SIGMA(當(dāng)使用 ISMEAR不等于 5 和 4時(shí) )的設(shè)置有關(guān)。要達(dá)到同樣的精度,較小的 SIGMA則需要較大的 kmesh;而且,當(dāng) SIGMA較小時(shí),若 kpoints不夠多,也會(huì)出現(xiàn)難收斂的情況。 NELM = 40 maximum of 40 electronic steps kmesh, SIGMA GGA=PS(選用的贗勢(shì)文件為 PBBEsol, =PE為 PBE的方法 ) LHFCALC = .True. PRECFOCK = Normal NKRED = 2 (設(shè)置此參數(shù)容易報(bào)錯(cuò),不知為何 ) TIME = HFSCREEN = AEXX = the exact exchange is used NPAR=16(使用核數(shù)為 4,16,64此種 ) ALGO=ALL 雜化泛函的計(jì)算 HF Functional 一般來說,現(xiàn)在都用新方法來計(jì)算極化,因?yàn)橛闷饋矸奖恪? 老方法的話還要把三個(gè)方向分來來計(jì)算,等于要跑三個(gè)任務(wù),相對(duì)復(fù)雜,但是準(zhǔn)確性可能好一些。 LBERRY=.TRUE. IGPAR=3 1,2,3對(duì)應(yīng) x, y, z方向 NPPSTR=7 DIPOL= EDIFF=1E5 Berry Phase計(jì)算極化,老方法 設(shè)置在計(jì)算離子的 dipole時(shí)的參考點(diǎn)即設(shè)置 DIPOL(注意的是,它的 設(shè)置需要使得原子移動(dòng)前后的原子都在這個(gè)參考點(diǎn)的一側(cè)。比如這