| 查看: 2273 | 回复: 1 | |||
脉脉如烟金虫 (小有名气)
|
[求助]
如何使用高斯ONIOM模型中mechanical embedding方法 已有1人参与
|
|
【求助】ONIOM的mechanical embedding方法 Gaussian的ONIOM模型对MOF材料QM/MM计算，考虑MM对QM电荷嵌入的时候，有mechanical embedding（法一） 和 electronic embedding（法二）两种方法。我先尝试两种方法优化初始复合物比较其对结构的影响发现： 法一的路径部分输入：#p ONIOM(B3LYP/6-31G*:UFF=QEQ) opt scfcon=6 scfcyc=400 optcyc=200 结果：.log文件发现MM部分的电荷均为0.0000。没加成功 请问1：使用mechanical embedding方法的正确输入是什么样子的呢；是需要提前计算或查找MM区片段电荷进行嵌入吗，如何实现呢？ 法二的路径部分输入：#p ONIOM(B3LYP/6-31G*:UFF=QEQ)=scalecharge=5 opt scfcon=6 scfcyc=400 optcyc=200 结果：.log文件发现MM部分均有电荷。但是优化过程中能量竟在10-2000kcal/mol的范围内强烈震荡！！ 请问2：我的体系QM部分有破坏共价键，electronic embedding电子嵌入不合适导致的能量震荡吗？ 请问3：体系中有溴阴离子，用什么基组比较合适呢？6-31+G*可以吗？ 以上三个问题，恳请帮忙啊。谢谢啊！！  |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有236人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关商家推荐: (我也要在这里推广)
xytk04
铁杆木虫 (正式写手)
- 应助: 35 (小学生)
- 金币: 7908.2
- 红花: 5
- 帖子: 615
- 在线: 1302小时
- 虫号: 345152
- 注册: 2007-04-14
- 专业: 理论和计算化学
【答案】应助回帖
|
>>>>>>法一的路径部分输入：#p ONIOM(B3LYP/6-31G*:UFF=QEQ) opt scfcon=6 scfcyc=400 optcyc=200 这样写是ME，没有问题。问题是你在用QEQ电荷，QEQ电荷的产生与初始结构有关，如果用这样写的话，你每次计算使用的电荷是不一样的，最后你计算出来的能量会出现问题的。 >>>>>>结果：.log文件发现MM部分的电荷均为0.0000。没加成功 你看到的电荷是Mulliken charges吧？如果是这样的话，MM部分的电荷为0.0000是正常的。你看看哪些原子有charge，是不是都是QM区域的原子？找找这个charge是怎么来的？ |
2楼2018-05-07 09:34:57
20
投票:
