计算化学公社

标题: 用gfn1-xtb优化一个单胞600多个原子的晶体时，scf不收敛 [打印本页]

作者
Author: 2532251334 时间: 2024-4-13 15:38
标题: 用gfn1-xtb优化一个单胞600多个原子的晶体时，scf不收敛
各位老师，学生优化一个单胞600多个原子的晶体时，scf一直不收敛，能量都变正了。想问下到底哪里出问题了。

作者
Author: wzkchem5 时间: 2024-4-13 16:51
可以先用GFN-FF初步优化一下，至少优化一下氢原子，这个结构里C-H、N-H键普遍太短

作者
Author: 2532251334 时间: 2024-4-13 17:21

wzkchem5 发表于 2024-4-13 16:51
可以先用GFN-FF初步优化一下，至少优化一下氢原子，这个结构里C-H、N-H键普遍太短

好的，谢谢回复

作者
Author: sobereva 时间: 2024-4-14 03:12
GFN1-xTB在CP2K里不用OT而用对角化的话很容易SCF不收敛
改成OT。对于分子晶体，当前晶胞尺寸不考虑k点勉强可以接受

作者
Author: 2532251334 时间: 2024-4-14 16:07

sobereva 发表于 2024-4-14 03:12
GFN1-xTB在CP2K里不用OT而用对角化的话很容易SCF不收敛
改成OT。对于分子晶体，当前晶胞尺寸不考虑k点勉强 ...

谢谢sob老师，就说很奇怪，用PBE都可以收敛，我试试。

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