24小时热门版块排行榜    

查看: 1797  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-28 15:41 by winsaint
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +12 晴天加油 2026-08-26 13/650 2026-08-28 11:32 by skyjudy
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见