计算化学公社

标题: ONIOM中的赝势机组的设置咨询，谢谢指点 [打印本页]

作者
Author: 小范范1989 时间: 2016-9-12 08:49
标题: ONIOM中的赝势机组的设置咨询，谢谢指点
1：我用ONIOM来计算，分了两层，效果图如下：
2：我想对cu采用lanl2dz基组，然后我的问题是，怎么样设置这个基组，同时，这个赝势机组是不是也应该对low层的分子设置？这个怎么弄？
3：我的输入文件如下，不知道对不对。谢谢指点
# ONIOM(pbe1pbe/genecp:UFF=QEq) geom=connectivity opt

Title Card Required

0 1
C             0 1.25490195 0.99999999 0.00000000 H
H             0 1.60997671 1.50383319 0.87462793 H
H             0 1.61317108 1.50496253 -0.87267179 H
H             0 0.18490374 1.00001380 -0.00195613 H
Cu             0 1.90154642 -0.82905740 0.00000000 H
H             0 0.55811712 -1.30401305 0.36126840 H
H             0 2.24215499 -0.70863874 1.42491584 H
H             0 2.39152960 -2.21499264 0.00000000 H
H             0 3.24497572 -0.35410176 -0.36126840 H
H             0 1.56093785 -0.94947607 -1.42491584 H
H             0 -0.04261461 1.46908569 1.19543184 L
Cu             0 -0.04260197 2.49441122 2.24881095 L
C             0 -1.98260197 2.49443512 2.24881095 L
H             0 -0.04261495 1.44103211 3.27413648 L
H             0 1.42739803 2.49439311 2.24881095 L
H             0 -0.04258934 3.51973675 3.30219007 L
H             0 -0.04258899 3.54779033 1.22348542 L
H             0 -2.33927456 1.98834340 3.12148262 L
H             0 -2.33925589 3.50324341 2.25076696 L
H             0 -2.33927502 1.99173173 1.37418325 L

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 1.0 8 1.0 9 1.0 10 1.0
6
7
8
9
10
11 12 1.0
12 13 1.0 14 1.0 15 1.0 16 1.0 17 1.0
13 18 1.0 19 1.0 20 1.0
14
15
16
17
18
19
20

cu 0
lanl2dz
****
C H 0
6-31G*
****

cu 0
lanl2dz

这个设置，对high的cu起作用，对low层也起作用吗？

还是说，这个low层，根本没必要管这些赝势机组什么的，因为用力场来衡量？
对这一块不是很了解，谢谢指点。

作者
Author: sobereva 时间: 2016-9-12 14:00
low层只牵扯到力场计算，跟基组无关，你也不用考虑这个。high层才同时牵扯到QM和MM计算

作者
Author: 小范范1989 时间: 2016-9-12 15:15

sobereva 发表于 2016-9-12 14:00
low层只牵扯到力场计算，跟基组无关，你也不用考虑这个。high层才同时牵扯到QM和MM计算

恩，明白了，谢谢sob老师的指点。
如果我这么计算的话，输入文件应该没有错误吧？就是写genecp，文件最后写基组部分和赝势部分，谢谢老师

作者
Author: sobereva 时间: 2016-9-12 19:24

小范范1989 发表于 2016-9-12 15:15
恩，明白了，谢谢sob老师的指点。
如果我这么计算的话，输入文件应该没有错误吧？就是写genecp，文件最 ...

没错

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