计算化学公社

标题: 求助：基于GFN1-xtb进行预平衡SCF总是不收敛 [打印本页]

作者
Author: Akaiconan 时间: 2022-6-11 17:00
标题: 求助：基于GFN1-xtb进行预平衡SCF总是不收敛
本帖最后由 Akaiconan 于 2022-6-11 17:00 编辑

社长和各位同仁好：
   我对含有100个原子的Na-Al-F体系进行基于GFN1-xtb理论的预平衡，灵感来源于另外一个社员的问题和社长的回答（http://bbs.keinsci.com/thread-24001-1-1.html）。尝试了OT和对角化算法，但是一直不收敛，震荡的很厉害。
   inp文件是Multiwfn生成的，除了改改eps和steps这种无关痛痒的参数外，未进行其他的改动。
   体系为packmol生成，30个NaF和10个AlF3混一起。之前对只含30个NaF和只含10个AlF3的盒子进行过基于GFN1-xtb的平衡，很顺利，但是混在一起就出问题了。

   在建立盒子时，我使盒子尽可能的致密（tolerance=2.0），我不太清楚这是不是造成不收敛的一个因素，如果不是，想请哪位大神不吝赐教一下不收敛的原因。

下方为输入文件和结构文件。
(, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 9)

作者
Author: sobereva 时间: 2022-6-11 17:34
这种DFT算起来不费劲的体系，而且和GFN-xTB原本参数化过程中考虑的体系差异太大的情况，就甭用GFN1-xTB了。否则就算能跑，可靠度也不高。来回折腾的功夫DFT都算完了。

作者
Author: Akaiconan 时间: 2022-6-11 18:42

sobereva 发表于 2022-6-11 17:34
这种DFT算起来不费劲的体系，而且和GFN-xTB原本参数化过程中考虑的体系差异太大的情况，就甭用GFN1-xTB了。 ...

明白了社长，谢谢指导，我就不折腾了。

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