24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1728  |  回复: 2

贺仪

铁杆木虫 (著名写手)

[求助] ONIOM电荷设置问题

最近我在计算分子内电荷转移的问题,想利用高斯的layer来计算,准备把一个分子分成三部分,这三部分都是QM的计算。
但是我看说明书一直没明白其中的一些概念。最主要的是不知道怎么填写电荷和自旋的那一栏。
说明书中是这么说的:
Values and defaults for three-layer ONIOM calculations follow an analogous pattern (in the subscripts below, the first character is one of: Real, Int=Intermediate system, and Mod=Model system, and the second character is one of: H, M and L for the High, Medium and Low levels):
cRealL  sRealL   [cIntM   sIntM   [cIntL  sIntL  [cModH  sModH   [cModM sModM   [cModL  sModL]]]]]

我能明白的就是这里的real system只的就是我们模拟的整个的体系,但是不明白那个Intermediate system 和Model system指的是什么,在网上搜了也没找见,看文章的话估计得话很多时间才能明白,各位同仁能不能指点下呢?
下面是我做的一个分子,分成三层,我想设置的是Low这层带一个负电荷,M层带一个正电荷,H层不带电荷,总体也是不带电荷的。大家觉得我应该怎么设置呢?

%nprocshared=4
%chk=E:\cooperation\cnn\SYS-after.chk
# oniom(b3lyp/6-31g(d):b3lyp/6-31g(d):b3lyp/6-31g(d))

sys

-1 2 1 2 1 2 0 1 0 1 0 1
C                0   -5.80978800   -0.70835900   -0.98827500 M
C                0   -4.69805200   -1.40248800   -0.56586200 M
C                0   -3.53582800   -0.71606100   -0.12291700 M
C                0   -3.53571800    0.71609900   -0.12332600 M
C                0   -4.69780600    1.40245400   -0.56674200 M
C                0   -5.80971100    0.70826300   -0.98858900 M
C                0   -2.36897800    1.39771200    0.31449900 M
C                0   -1.25423500    0.71369800    0.74485100 M H 14
C                0   -1.25430400   -0.71351500    0.74514700 M H 11
C                0   -2.36918700   -1.39760300    0.31527100 M
C                0   -0.00000400   -1.41705800    1.22539400 H
C                0    1.25418500   -0.71353900    0.74481600 L H 11
C                0    1.25427300    0.71367600    0.74486600 L H 14
C                0    0.00002100    1.41731600    1.22510900 H
C                0    2.36922400    1.39767100    0.31501900 L
C                0    3.53581900    0.71603300   -0.12313800 L
C                0    3.53560200   -0.71612400   -0.12351900 L
C                0    2.36889700   -1.39764100    0.31453400 L
C                0    4.69768300   -1.40257200   -0.56679500 L
C                0    5.80973100   -0.70846900   -0.98841600 L
C                0    5.80993800    0.70814700   -0.98806800 L
C                0    4.69811300    1.40236400   -0.56605000 L
H                0   -6.69272400   -1.24575500   -1.32386400 M
H                0   -4.69490400   -2.49020600   -0.56581300 M
H                0   -4.69456200    2.49017000   -0.56712600 M
H                0   -6.69249900    1.24560600   -1.32464900 M
H                0   -2.36627500    2.48628700    0.30286100 M
H                0   -2.36660200   -2.48618400    0.30411800 M
H                0    0.00016800   -1.42797900    2.32895900 H
H                0   -0.00008100   -2.46641500    0.91176500 H
H                0    0.00001000    2.46660900    0.91129000 H
H                0    0.00003700    1.42844500    2.32867700 H
H                0    2.36673100    2.48625200    0.30387300 L
H                0    2.36616200   -2.48621800    0.30304300 L
H                0    4.69437800   -2.49028800   -0.56710200 L
H                0    6.69264100   -1.24588400   -1.32404100 L
H                0    6.69306500    1.24547000   -1.32327000 L
H                0    4.69521300    2.49008100   -0.56558100 L
回复此楼
嘟啦啦嘟啦啦嘟
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

★ ★
gmy1990: 金币+2 2013-01-07 22:06:58
你的ONIOM分层电荷和自旋指定应该没有问题,因为GaussianView打开后可以看到的确显示的是低层电荷-1,中层电荷+1,就相信系统默认的吧。不过我不确定你用分层计算是否能说明电荷转移的问题。
2楼2013-01-07 21:50:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by kaegi at 2013-01-07 21:50:26
你的ONIOM分层电荷和自旋指定应该没有问题,因为GaussianView打开后可以看到的确显示的是低层电荷-1,中层电荷+1,就相信系统默认的吧。不过我不确定你用分层计算是否能说明电荷转移的问题。

谢谢,昨天发完贴不就,就意识到我理解错了,layer的概念与我想的完全不一样。
现在在试fragment,对高斯不熟,摸起来慢呀,
谢谢回复!
嘟啦啦嘟啦啦嘟
3楼2013-01-08 09:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 4/200 2026-05-15 17:43 by x0mp7owy2b
[教师之家] 上海大学实验技术岗位非升即走 +3 嘻嘻哈哈乐呵呵 2026-05-15 3/150 2026-05-15 15:40 by 187385hgh
[论文投稿] 有带发论文的吗 +4 山楂之术 2026-05-09 4/200 2026-05-15 15:40 by 妹子不好惹
[基金申请] 青C资助名额大幅增加! +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 12:03 by onwj4wpxp2
[基金申请] 精华III评审感受-评审感受-评审感受 +14 ferrarichen 2026-05-11 18/900 2026-05-15 11:12 by cmhchen
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +9 dd921ww 2026-05-12 10/500 2026-05-15 10:41 by muyiliuhui
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 3/150 2026-05-15 09:23 by onwj4wpxp2
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 cjf4bx70cj 2026-05-14 4/200 2026-05-15 09:03 by gagyerk94e
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[考博] 申博自荐 +4 食品的橙子 2026-05-09 6/300 2026-05-14 16:05 by great1919
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
信息提示
请填处理意见