多谢指点！实际上我做个测试，将来扩胞用到大的体系上，然后上面吸附团簇，进而吸附小分子催化。
这样的话，DFT计算就十分吃力了，因此采用QMMM，基底用MM，上面活性部分用QM。
这时候，运行速度会大大加快，可以搜索到更多构型。

这么小的体系，鼓捣分子力场设置的功夫，DFT都跑完了，结果还更精确

CP2K做MM计算要掌握很多知识和逻辑才能正确使用，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里面给了详细的介绍和诸多例子，其中和当前问题相关的一页：

不是，难道计划研究的体系尺寸大到DFT算不动、研究目标又不需要对电子状态描述准确的结论吗？（而且DFT这种量子化学方法的计算不能简称为“量子计算”，看http://sobereva.com/680）就算DFT算不动也还有GFN-xTB呢。

这种结构在CP2K里恰当设置DFT不可能处理不了，真有什么问题不如直接拿出来（比如因为SCF收敛困难而算得慢，就按http://sobereva.com/665仔细排查），而不是绕弯路折腾什么经典分子力场；即便模拟很快，把找参数、设置输入的准备流程也算上可能根本就省不下多少时间，而且结果能否被认可而发表都十分存疑。单说相互作用势函数及单位的设置就有GENPOT字段的一堆讲究，真心不建议缺乏经验就贸然尝试。