计算化学公社

标题: 求助QM/MM边界连接原子报错问题 [打印本页]

作者
Author: romiishere 时间: 2024-9-12 22:53
标题: 求助QM/MM边界连接原子报错问题
本帖最后由 romiishere 于 2024-9-12 23:36 编辑

请教一下各位，我在对用cp2k做qmmm的时候报错了CPASSERT faild，我的体系sapphire0001表面，XY方向周期性的，并且将在之后对吸附进行研究，想先试着跑一下没有吸附的表面

(, 下载次数 Times of downloads: 7) (, 下载次数 Times of downloads: 9)

(, 下载次数 Times of downloads: 8)

经检查输入文件和多次测试，我发现是在QMMM边界的连接原子出现了问题，例如下面这种情况

如果把选中的部分删除则程序能够正常运行，但是一但有不同的link模块指向了同一个MM原子就会报错，指向同一个QM原子不会出现问题。想请问一下

1是不是对连接原子的处理就只能对一个MM原子连接一个QM原子？

2如果是这样的话，对于这种有多个共价键的边界应该怎么处理呢？

不胜感激！

输入文件如下，为了排查错误只留下了两个link模块其它的link都删掉了

TYPE设成GHO也一样

(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1)

作者
Author: sobereva 时间: 2024-9-13 01:59
当前这样需要截断半离子键的情况实际上并不怎么适合QM/MM。建议尽量直接DFT算整体，更精确还省得折腾

作者
Author: romiishere 时间: 2024-9-13 10:28

sobereva 发表于 2024-9-13 01:59
当前这样需要截断半离子键的情况实际上并不怎么适合QM/MM。建议尽量直接DFT算整体，更精确还省得折腾

好的好的，十分感谢！

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