计算化学公社

标题: 求助：xTB运行程序优化过程崩溃 [打印本页]

作者
Author: SpongeBob 时间: 2021-6-7 22:38
标题: 求助：xTB运行程序优化过程崩溃
本帖最后由 SpongeBob 于 2021-6-8 08:55 编辑

我在尝试用xTB计算时出现了一些问题。我的体系中包含五个BH4-，一个TFSI-，和一个Mg2+，应用命令xtb xtb.xyz --chrg -4 --uhf 0 --gfn 1 --omd -I md.inp，却无法运行，在优化时就会崩溃如下图，请问一下是什么原因啊？ (, 下载次数 Times of downloads: 28)
用gfn0和gfn2可以优化成功，但是结果似乎不是很合理。
(, 下载次数 Times of downloads: 28)

(, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 3)

作者
Author: hdhxx123 时间: 2021-6-7 22:49
贴的图里写了，结构不合理不收敛，建议先用gfn0优化或者其他半经验方法优化一下。

作者
Author: SpongeBob 时间: 2021-6-8 08:54

hdhxx123 发表于 2021-6-7 22:49
贴的图里写了，结构不合理不收敛，建议先用gfn0优化或者其他半经验方法优化一下。

我尝试过用gfn0和gfn2进行优化，都可以优化成功，但是优化的结果似乎不合理。我把图放在了第一层。

作者
Author: sobereva 时间: 2021-6-8 19:39
体系里那么多阴离子，整体带显著负电荷，在真空下可不容易通过静电斥力崩飞了

作者
Author: SherryLiu 时间: 2021-6-10 06:09
下图是我用gfn 2 优化的结构，但在md.inp 里加了势函数以保持离子间的相互作用，结构似乎比较合理。
另外不知道这个体系是否需要考虑溶剂效应，楼主可以试试看是否有所差异。
此外，gfn 0, 1, 2 可能适合某种体系可能不适合某种体系，需要根据具体的研究的体系确定。
$wall
potential=logfermi
sphere: auto, all
$end
$cma

作者
Author: SpongeBob 时间: 2021-6-10 08:26

sobereva 发表于 2021-6-8 19:39
体系里那么多阴离子，整体带显著负电荷，在真空下可不容易通过静电斥力崩飞了

嗯嗯，谢谢sob老师。

作者
Author: SpongeBob 时间: 2021-6-10 08:27

SherryLiu 发表于 2021-6-10 06:09
下图是我用gfn 2 优化的结构，但在md.inp 里加了势函数以保持离子间的相互作用，结构似乎比较合理。
另外 ...

产生这个现象的原因很可能是溶剂效应，我会去尝试一下，谢谢

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