计算化学公社

标题: 关于Oniom QM/MM输入文件的问题：1.能量一直震荡 2.基态能量和激发态能量差不多 [打印本页]

作者
Author: 赤赤媛 时间: 2020-11-11 10:04
标题: 关于Oniom QM/MM输入文件的问题：1.能量一直震荡 2.基态能量和激发态能量差不多
第一次做Oniom ，开始发现能量震荡99999，然后用了 opt=gdiis ， opt(gdiis,maxstep=4,notrust)，还是老样子：震荡。后来发现基态和激发态能量相似。请问是输入文件出现什么问题了吗？如何改？

作者
Author: 赤赤媛 时间: 2020-11-11 10:13
基态输入文件

作者
Author: 赤赤媛 时间: 2020-11-11 10:15
激发态输入文件

作者
Author: biogon 时间: 2020-11-11 10:28
加大格点精度，要不然用其他泛函，基组更没必要用这么大的

作者
Author: 赤赤媛 时间: 2020-11-11 10:43

biogon 发表于 2020-11-11 10:28
加大格点精度，要不然用其他泛函，基组更没必要用这么大的

谢谢，还想问一下请问输入文件有错吗？特别是激发态输入文件，可以写成td/M062X/defTZVP吗

作者
Author: 让你变成回忆 时间: 2020-11-11 11:13
我觉得输入文件不妥，
（1）你这里使用的力学嵌入的方式，你确定考虑力学嵌入可以，而不是考虑电学嵌入？
（2）如果是电学嵌入，输入文件既没有在元素名处指定原子电荷，也没有在UFF处告诉程序用什么方式计算原子电荷，用默认的电荷是否可靠？（这里是用默认的还是不用电荷我忘了，之前有测试过）
（3）激发态应该写成：td m062x/defTZVP，不确定你这样写是否有问题。

作者
Author: Frank 时间: 2020-11-11 17:56
要用Qeq, #p oniom(m062x/def2TZVP TD (NStates=10):UFF=QEq)=Embed

作者
Author: sobereva 时间: 2020-11-11 21:24
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充

不要随便选择“量化理论”分类，给你改了。仔细看本版的版头了解各个分类都适合什么帖子！

本来M06-2X就容易震荡，算价层激发能严重偏高，而且耗时还高，干嘛要用M06-2X。没特殊情况就改成PBE0-D3(BJ)

作者
Author: 赤赤媛 时间: 2020-11-12 22:30

sobereva 发表于 2020-11-11 21:24
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充

不要随便选择“量化 ...

好的，记住了，社长

作者
Author: 赤赤媛 时间: 2020-11-12 22:33

让你变成回忆发表于 2020-11-11 11:13
我觉得输入文件不妥，
（1）你这里使用的力学嵌入的方式，你确定考虑力学嵌入可以，而不是考虑电学嵌入？
...

我才接触这种方法，大神，能不能说一下力学嵌入和电学嵌入的区别

作者
Author: 赤赤媛 时间: 2020-11-12 22:37

Frank 发表于 2020-11-11 17:56
要用Qeq, #p oniom(m062x/def2TZVP TD (NStates=10):UFF=QEq)=Embed

这种是电学嵌入的关键词吗？谢谢

作者
Author: wzkchem5 时间: 2020-11-12 23:24

赤赤媛发表于 2020-11-12 22:33
我才接触这种方法，大神，能不能说一下力学嵌入和电学嵌入的区别

电学嵌入考虑了MM区和QM区之间的静电相互作用，力学嵌入不考虑

作者
Author: sobereva 时间: 2020-11-13 08:57

赤赤媛发表于 2020-11-12 22:33
我才接触这种方法，大神，能不能说一下力学嵌入和电学嵌入的区别

(, 下载次数 Times of downloads: 69)
(, 下载次数 Times of downloads: 72)

作者
Author: 赤赤媛 时间: 2020-11-13 10:57

wzkchem5 发表于 2020-11-12 23:24
电学嵌入考虑了MM区和QM区之间的静电相互作用，力学嵌入不考虑

谢谢老师

作者
Author: 赤赤媛 时间: 2020-11-13 10:59

sobereva 发表于 2020-11-13 08:57

谢谢社长

欢迎光临计算化学公社 (http://bbs.keinsci.com/)