|
|
各位老师:
大家好!
我最近在做一个溶剂里面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盒子 ?
多谢各位大神指点
|
|