CP2K分子动力学计算 力场参数单位如何确定？

lywbanner

大家好：
  量子计算有时候很费时间，在要求不高的情况下，用分子动力学(FIST)模拟也是个不错的选择。
不过力场设置不好，总会出现
             “GEOMETRY wrong or EMAX_SPLINE too small!”
这样错误。本人总是遇到。多方查找，也不知道问题出在哪。
   首先，模型是很简单的（初学测试）Al2O3-(100) slab,固定下面两层，其余弛豫。

其次，力场是公开发表的文献中的 ：




采用了GENPOT力场（详情见inp文件）。参数我都进行了仔细的核对，应该问题不大。
如果实在有问题，那就是单位谁设（详情见inp文件）。其他方面，实在想不出会出什么错
（inp文件是Multiwfn生成，几乎没怎么改动）。


所以特请老师、学友们看一看，是不是单位没设的原因？GENPOT力场单位究竟该怎么设？
谢谢大家!!!
MD Geometry inp.inp (4.87 KB, 下载次数 Times of downloads: 3)

2025-2-21 23:28 上传 Uploaded

点击下载
Click to download
inp文件

   

lywbanner

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

sobereva

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

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

Uus/pMeC6H4-/キ

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

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