本帖最后由 蔓莉青萝 于 2019-6-1 16:01 编辑 卡开发发 发表于 2019-6-1 00:40 好的,谢谢,我试试 |
蔓莉青萝 发表于 2019-5-31 20:58 计算频率难收敛是因为DMol3采用的是有限差分方法算频率,即通过计算6*原子数个力差分来得到Hessian。这6*原子数的力是通过原子分别沿着x y z三个方向正负分别移动产生结构并计算每个原子的力得到的,这些结构里面有一部分是SCF比较难收敛。我个人经验上,尤其是一些有磁性的窄带体系或是本来结构就比较怪异的结构(基本上都是那种多参考效应比较强的)SCF本来就比较难收敛,在进行移动移动原子时产生的某些结构确实更加难收敛。 只是增加迭代次数真的没用(参考http://sobereva.com/61的第2节第17点)。对于DMol3,smearing过大既不可靠也不见得解决问题,在精度可以保证的基础上调整density mixing和DIIS size,实在不行还可以适当调整下有限差分的步长,其余没什么特别好的办法。 |
| 参与人数Participants 1 | eV +2 | 收起 理由Reason |
|---|---|---|
|
| + 2 |
kkwang 发表于 2019-5-30 08:53 我scf已经设置到500了,感觉不是这个问题 |
| 增大迭代次数 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-25 03:31 , Processed in 0.167968 second(s), 26 queries , Gzip On.