计算化学公社

标题: 有关ONIOM的电荷和多重度的设置问题 [打印本页]

作者
Author:
lei234    时间: 2016-3-30 07:18
标题: 有关ONIOM的电荷和多重度的设置问题
以上传的输入和输出文件为例,请教ONIOM的电荷和多重度是如何设置的?三组电荷和多重度分别对应的是什么?
log文件中提示的455个电子指的是哪里的?
整个体系gv里显示的电荷是-2,多重度是1,为什么log文件说多重度有问题?

感谢~~

作者
Author:
sobereva    时间: 2016-4-3 14:33
如果是ONIOM分了两层,则需要指定三个电荷和多重度,分别是低水平全部,高水平局部,低水平局部。注意同样区域的电荷多重度应当一样,比如高水平局部和低水平局部的电荷和多重度必须一样。

报错原因这里写了
ONIOM: generating point  2 -- high level on model system.
Standard basis: 6-31G(d) (6D, 7F)
The combination of multiplicity 1 and   455 electrons is impossible.
显然正在用DFT算局部体系

当前这种划分很糟糕。碳纳米管甭用ONIOM来分割,共轭全都破坏了。

作者
Author:
lei234    时间: 2016-4-4 01:24
sobereva 发表于 2016-4-3 14:33
如果是ONIOM分了两层,则需要指定三个电荷和多重度,分别是低水平全部,高水平局部,低水平局部。注意同样 ...

谢谢您的回复~我原来写的是-2 1 -2 1 -2 1,所以应该改成-2 1 -2 2 -2 2对吗?
这种碳纳米管上吸附小分子的体系,只能用您说的MOPAC吗?我试过了您说的全部用DFT结合混合基组,实在是太慢了。。。不知您还有什么建议?
作者
Author:
sobereva    时间: 2016-4-4 01:27
lei234 发表于 2016-4-4 01:24
谢谢您的回复~我原来写的是-2 1 -2 1 -2 1,所以应该改成-2 1 -2 2 -2 2对吗?
这种碳纳米管上吸附小分 ...


你可以试试,不过ONIOM怎么设对此体系都不理想。即便看似以能接受的耗时算完,结果也不好,折腾来折腾去反倒时间都耽误了。

想进一步节约时间,可以用ORCA用纯泛函算,参见
大体系弱相互作用计算的解决之道
http://sobereva.com/214
或者把纳米管尺寸减小点,并不需要那么长的管子。
很推荐在MOPAC里用PM7预优化一下。
作者
Author:
lei234    时间: 2016-4-4 06:06
sobereva 发表于 2016-4-4 01:27
你可以试试,不过ONIOM怎么设对此体系都不理想。即便看似以能接受的耗时算完,结果也不好,折腾来折腾 ...

好的,多谢~
作者
Author:
cwindy    时间: 2016-4-4 11:21
改成-2 1 -2 2 -2 2应该能运行

但是你确认你的高层的多重度是2么?
作者
Author:
lei234    时间: 2016-4-4 23:51
cwindy 发表于 2016-4-4 11:21
改成-2 1 -2 2 -2 2应该能运行

但是你确认你的高层的多重度是2么?

455个电子多重度不是最小是2么?
作者
Author:
cwindy    时间: 2016-4-7 15:57
lei234 发表于 2016-4-4 23:51
455个电子多重度不是最小是2么?

一般算激发态才会出现多重度大于1的情况

我一般都是这样认为的,不知道对不对




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3