请教CASTEP计算弹性常数失败的原因
请高手们帮忙解答:
下面是我在计算弹性常数时在C21.castep里出现的提示:
*Warning* max. SCF cycles performed but system has not reached the groundstate.
Current total energy, E = -10221.90318714 eV
Current free energy (E-TS) = -10221.90428269 eV
(energies not corrected for finite basis set)
NB est. 0K energy (E-0.5TS) = -10221.90373491 eV
****************************************************************************
Warning: electronic minimisation did not converge when finding ground state.
****************************************************************************
Writing checkpoint file...
Writing model to 3D_Atomistic_cij__2__1.check
Error check_elec_ground_state : electronic_minimisation of initial cell failed.
计算的结构比较复杂,正交结构,原来的SCF为100,现在只增大SCF能否使计算完成?有没有人遇到过这种情况?该如何修改参数呢?
返回小木虫查看更多
京公网安备 11010802022153号
根据这个提示*Warning* max. SCF cycles performed but system has not reached the groundstate.应该增大SCF cycles。你的体系这么大怎么也应该调到500那样吧。你可以逐步往上增大,直至收敛
还有一个问题不明白,为什么不能将SCF 值设的很大而要逐渐收敛呢?为什么不能一开始就设1000呢,是不是SCF值设的越大计算的就越慢?
试试减小density mixing参数. 打开Castep calculation-Electronic, 点击右下角more..., 出现新参数框, 在菜单栏点击SCF, 找density mixing, 把Charge后面的参数改小些, 比如0.05或更小. 100步自洽不收敛, 我估计再增加迭代次数也是无济于事的.
楼上能否解释一下为什么修改那个参数?我算的这个体系没有十天半月是算不完的,所以很想知道修改参数后的可行性,非常感谢!
自洽场(SCF)计算时, 程序利用上两步计算得到的电荷密度进行混合得到新的输入电荷密度, mix_charge就是混和参数, 如果混合参数太大, 会引起自洽场计算过程中, 电荷密度有较大的振荡而不收殓.
另外, 在计算很大的体系的时候, 最好适当降低收殓精度, 否则, 收殓可能会有困难,
非常感谢楼上的回答。我前天按照你的建议将mix_charge减小为0.05,同时将SCF增大为200,但是计算了两天在Initial geometry时就失败了,也没有给出失败的原因。现在是不是只能减小收敛的精度了(原来是10-6)?这样的话还要重新定结构,很郁闷.........我再试一下
很是郁闷,已经半个月过去了,还是没有一点进展。
精度降低了还是不收敛,错误跟原来一模一样。金属体系就这么难收敛吗,还是有别的原因?