计算化学公社

标题: 求助：关于cp2k进行QMMM计算时，MM力场和QM区域设置的问题 [打印本页]

作者
Author: haisenburg 时间: 2022-10-6 10:30
标题: 求助：关于cp2k进行QMMM计算时，MM力场和QM区域设置的问题
各位老师：
大家好！
我最近在做一个溶剂里面Zn团簇的QMMM计算！溶剂采用MM，Zn(10个原子)团簇采用QM(gpw,DFT)，计算过程有如下两个疑问：
1. 关于MM力场的设置问题
我的溶剂是LiCl-KCl熔盐，规定了Li-Li、Li-K、K-K、Li-Cl、K-Cl、Cl-Cl的力场参数（BMHFT势函数），计算时候报错，说是没有规定 Zn-K、Zn-Li、Zn-Cl、Zn-Zn的关键力场参数。关于 Zn-K、Zn-Li、Zn-Cl三种相互力场没有规定，我能理解，因为涉及溶剂和团簇之间的相互作用。但是Zn-Zn是在QM区域，采用DFT计算，为何还需要规定Zn-Zn的力场，如果这里规定了Zn-Zn，那么QM区域计算什么？其与QM区域的DFT计算不是冲突吗？

2. 关于QM区域设置的问题
2.1 我的模型由于是团簇溶解在溶剂里面，如果将QM区域设置成与MM盒子一样大小，计算太耗时，我就想将QM区域设置成MM盒子的一半大小，但是老是报错“ GEOMETRY wrong or EMAX_SPLINE too small! ”，不管我是否将QM区域的周期性打开，还是加上WALLS，均报错！查了一下结构，确实是两个原子考的太近了导致的几何结构错误！inp和out文件见附件

2.2 另一个疑问是：QM区域的盒子是以指定进行QM计算的原子为中心的盒子，还是以MM盒子为参考，建立的相对于MM盒子原点建立的QM盒子？

多谢各位大神指点

作者
Author: 名字长的厉害 时间: 2022-10-12 15:12
QMMM教程里有说，QM算的部分，MM部分也要声明，ε设为0就行了

作者
Author: highlight 时间: 2022-10-12 16:27
2.2 ：关于 QM 盒子中心见如下 manual 的 CENTER/CENTER_TYPE
https://manual.cp2k.org/cp2k-202 ... ORCE_EVAL/QMMM.html

作者
Author: tbw@zyc 时间: 2022-12-9 10:03
也是“ GEOMETRY wrong or EMAX_SPLINE too small! ”，请问楼主解决了吗？

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