Dmol3频率计算问题

蔓莉青萝 · 发表于 Post on 2019-5-18 19:44:52

请教一下各位，我用Dmol3进行结构优化成功，一些性质计算也可以，但是将优化好的结构进行频率计算的时候，一直失败说scf不收敛，或者如下图所示

，无论我是降低精度，增加smearing值，或者降低density mixing Charge值都无法解决问题，不知道各位有没有什么好的解决办法？

蔓莉青萝 · 发表于 Post on 2019-6-1 15:45:47

本帖最后由蔓莉青萝于 2019-6-1 16:01 编辑

卡开发发发表于 2019-6-1 00:40
计算频率难收敛是因为DMol3采用的是有限差分方法算频率，即通过计算6*原子数个力差分来得到Hessian。这6* ...

好的，谢谢，我试试

卡开发发 · 发表于 Post on 2019-6-1 00:40:00

蔓莉青萝发表于 2019-5-31 20:58
我scf已经设置到500了，感觉不是这个问题

计算频率难收敛是因为DMol3采用的是有限差分方法算频率，即通过计算6*原子数个力差分来得到Hessian。这6*原子数的力是通过原子分别沿着x y z三个方向正负分别移动产生结构并计算每个原子的力得到的，这些结构里面有一部分是SCF比较难收敛。我个人经验上，尤其是一些有磁性的窄带体系或是本来结构就比较怪异的结构（基本上都是那种多参考效应比较强的）SCF本来就比较难收敛，在进行移动移动原子时产生的某些结构确实更加难收敛。

只是增加迭代次数真的没用(参考http://sobereva.com/61的第2节第17点）。对于DMol3，smearing过大既不可靠也不见得解决问题，在精度可以保证的基础上调整density mixing和DIIS size，实在不行还可以适当调整下有限差分的步长，其余没什么特别好的办法。

蔓莉青萝 · 发表于 Post on 2019-5-31 20:58:41

kkwang 发表于 2019-5-30 08:53
增大迭代次数

我scf已经设置到500了，感觉不是这个问题

kkwang · 发表于 Post on 2019-5-30 08:53:18

增大迭代次数

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[其它量化程序] Dmol3频率计算问题

评分 Rate