有关ONIOM的电荷和多重度的设置问题

lei234 · 发表于 Post on 2016-3-30 07:18:22

以上传的输入和输出文件为例，请教ONIOM的电荷和多重度是如何设置的？三组电荷和多重度分别对应的是什么？
log文件中提示的455个电子指的是哪里的？
整个体系gv里显示的电荷是-2，多重度是1，为什么log文件说多重度有问题？

感谢~~

sobereva · 发表于 Post on 2016-4-3 14:33:55

如果是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来分割，共轭全都破坏了。

lei234 · 发表于 Post on 2016-4-4 01:24:41

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

谢谢您的回复~我原来写的是-2 1 -2 1 -2 1，所以应该改成-2 1 -2 2 -2 2对吗？
这种碳纳米管上吸附小分子的体系，只能用您说的MOPAC吗？我试过了您说的全部用DFT结合混合基组，实在是太慢了。。。不知您还有什么建议？

sobereva · 发表于 Post on 2016-4-4 01:27:23

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预优化一下。

lei234 · 发表于 Post on 2016-4-4 06:06:20

sobereva 发表于 2016-4-4 01:27
你可以试试，不过ONIOM怎么设对此体系都不理想。即便看似以能接受的耗时算完，结果也不好，折腾来折腾 ...

好的，多谢~

cwindy · 发表于 Post on 2016-4-4 11:21:11

改成-2 1 -2 2 -2 2应该能运行

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

lei234 · 发表于 Post on 2016-4-4 23:51:45

cwindy 发表于 2016-4-4 11:21
改成-2 1 -2 2 -2 2应该能运行

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

455个电子多重度不是最小是2么？

cwindy · 发表于 Post on 2016-4-7 15:57:17

lei234 发表于 2016-4-4 23:51
455个电子多重度不是最小是2么？

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 有关ONIOM的电荷和多重度的设置问题