计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2519|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助:关于cp2k进行QMMM计算时,MM力场和QM区域设置的问题

[复制链接 Copy URL]

15

帖子

0

威望

485

eV
积分
500

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
各位老师:
大家好!
我最近在做一个溶剂里面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盒子 ?

多谢各位大神指点

geo-opt.inp

7.96 KB, 下载次数 Times of downloads: 36

LiKCl800-Zn10.cif.xyz

25.21 KB, 下载次数 Times of downloads: 8

geo-opt-out文件.rar

73.41 KB, 下载次数 Times of downloads: 4

11

帖子

0

威望

125

eV
积分
136

Level 2 能力者

2#
发表于 Post on 2022-10-12 15:12:25 | 只看该作者 Only view this author
QMMM教程里有说,QM算的部分,MM部分也要声明,ε设为0就行了

186

帖子

1

威望

4551

eV
积分
4757

Level 6 (一方通行)

3#
发表于 Post on 2022-10-12 16:27:59 | 只看该作者 Only view this author
2.2 : 关于 QM 盒子中心见如下 manual 的 CENTER/CENTER_TYPE
https://manual.cp2k.org/cp2k-202 ... ORCE_EVAL/QMMM.html

25

帖子

0

威望

404

eV
积分
429

Level 3 能力者

4#
发表于 Post on 2022-12-9 10:03:01 | 只看该作者 Only view this author
也是“ GEOMETRY wrong or EMAX_SPLINE too small!  ”, 请问楼主解决了吗?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 20:52 , Processed in 0.168843 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list