24小时热门版块排行榜    

查看: 1798  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +4 gouxfjh 2026-08-28 4/200 2026-08-29 04:46 by alongwaytogo
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 04:18 by gy1nBQXYQJqL
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 03:58 by gy1nBQXYQJqL
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 00:28 by winsaint
[基金申请] 国自然评审意见 +11 wangmingqi 2026-08-28 15/750 2026-08-29 00:17 by fangyl2005
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见