沾沾自喜的近义词PWSCF⾃洽计算(Self Consistency)中不收敛问题pwscf⾃洽计算过程中电⼦部分设置为:
水引&electrons
diagonalization = 'david'
萝卜鸡蛋饼conv_thr = 1.0d-9
mixing_beta = 0.7
其中diagonalization⽤来设置在求KS⽅程的本征⽮量和本征值时,采⽤具体的什么算法,赋值为'david'表⽰采⽤Davidson iterative diagonalization with overlap matrix⽅法;(Davidson iterative diagonalization with overlap matrix (default). Fast, may in some rare cas fail.)
conv_thr⽤来设置⾃洽收敛标准,赋值为⾃洽循环过程总能的变化⼩于1.0e-9时,⾃洽计算停⽌;
mixing_beta⽤来设置⾃洽计算过程中前后两次电荷密度混合的参数。(default,0.7)。
但是在输出⽂件中出现了不收敛的情况,即
iteration # 6 ecut= 50.00 Ry beta=0.70
圣经主祷文
Davidson diagonalization with overlap
c_bands: 1 eigenvalues not converged
一座塔
c_bands: 1 eigenvalues not converged
大学生自我评价100字c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
c_bands: 1 eigenvalues not converged
解决办法:1.pwscf帮助⽂件指出:
小学生优秀日记大全
Self-consistency is slow or does not converge at all In most cas: your input data is bad, or el your system is metallic and you are treating it as an insulator. If this is not the ca: reduce mixing beta to 0.3 ÷ 0.1 or smaller, try the mixing mode value that is more appropriate for your problem.(如果不是input结构不好或⾦属性当成了绝缘体这种错误,可以通过降低mixing beta这个值来调整)
2.⽹上意见:a)这个warning只要不是最后⼀个电⼦步是可以忽略的。
b)混合因⼦mixing_beta的问题,⼀般来说越⼩越趋向于收敛,尤其是对收敛难的体系,应该设置更⼩的值。我都⽤0.1。
c)这是说你的hamiltonian的本征值计算不收敛。这不影响你的结果,只要你最终得到收敛结果就可以。
d)个⼈意见可以先做⼀步优化relax,再做scf计算。
⽬前正在调试中,待续。。。。
全球色