计算化学公社

标题: Dmol3频率计算问题 [打印本页]

作者
Author: 蔓莉青萝 时间: 2019-5-18 19:44
标题: Dmol3频率计算问题
请教一下各位，我用Dmol3进行结构优化成功，一些性质计算也可以，但是将优化好的结构进行频率计算的时候，一直失败说scf不收敛，或者如下图所示 (, 下载次数 Times of downloads: 30) ，无论我是降低精度，增加smearing值，或者降低density mixing Charge值都无法解决问题，不知道各位有没有什么好的解决办法？

作者
Author: kkwang 时间: 2019-5-30 08:53
增大迭代次数

作者
Author: 蔓莉青萝 时间: 2019-5-31 20:58

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

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

作者
Author: 卡开发发 时间: 2019-6-1 00:40

蔓莉青萝发表于 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，实在不行还可以适当调整下有限差分的步长，其余没什么特别好的办法。

作者
Author: 蔓莉青萝 时间: 2019-6-1 15:45
本帖最后由蔓莉青萝于 2019-6-1 16:01 编辑

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

好的，谢谢，我试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)