应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » ONIOM电荷设置问题
31/1返回列表
查看: 1770  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

贺仪

铁杆木虫 (著名写手)

[求助] 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
回复此楼

» 收录本帖的淘帖专辑推荐

过渡态计算

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
嘟啦啦嘟啦啦嘟
1楼 2013-01-07 13:54:45
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 贺仪 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于代码变化问题，想知道的进来 +5 且听虎啸 2026-08-07 8/400 2026-08-07 17:17 by gltch
[基金申请] 关于filecode +3 布布和一二 2026-08-07 4/200 2026-08-07 16:34 by 医学老男孩
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了，论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +23 低垂的野花 2026-07-31 31/1550 2026-08-07 12:22 by speculiar
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭