24小时热门版块排行榜    

查看: 763  |  回复: 2

z5476t4508

新虫 (正式写手)

明天家

[交流] 【求助】关于ONIOM输入文件 已有1人参与

初学ONIOM,摆脱大家指点一下!
我在高斯手册里发现一个例子:

# ONIOM(B3LYP/6-31G(d):Amber) Geom=Connectivity

2 layer ONIOM job

0 1 0 1 0 1   Charge/spin for entire molecule (real system), model system-high level & model-low
C-CA--0.25   0   -4.703834   -1.841116   -0.779093 L
C-CA--0.25   0   -3.331033   -1.841116   -0.779093 L H-HA-0.1  3
C-CA--0.25   0   -2.609095   -0.615995   -0.779093 H
C-CA--0.25   0   -3.326965    0.607871   -0.778723 H
C-CA--0.25   0   -4.748381    0.578498   -0.778569 H
C-CA--0.25   0   -5.419886   -0.619477   -0.778859 L H-HC-0.1  5
H-HA-0.1     0   -0.640022   -1.540960   -0.779336 L
H-HA-0.1    0   -5.264565   -2.787462   -0.779173 L
H-HA-0.1    0   -2.766244   -2.785438   -0.779321 L
C-CA--0.25  0   -1.187368   -0.586452   -0.779356 L H-HA-0.1  3
C-CA--0.25  0   -2.604215    1.832597   -0.778608 H
H-HA-0.1    0   -5.295622    1.532954   -0.778487 L H-HA-0.1  5
H-HA-0.1    0   -6.519523   -0.645844   -0.778757 L
C-CA--0.25  0   -1.231354    1.832665   -0.778881 L H-HC-0.1  11
C-CA--0.25  0   -0.515342    0.610773   -0.779340 L
H-HA-0.1    0   -3.168671    2.777138   -0.778348 L H-HA-0.1  11
H-HA-0.1    0   -0.670662    2.778996   -0.779059 L
H-HA-0.1    0    0.584286    0.637238   -0.779522 L

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

如果 01 01 01 分别指的是整个体系的电荷、自旋;量子区域的电荷和自旋的话,那么也就是说整个体系和量子区域的总电量都是0啦,可是吧输入文件里的原子电荷加起来并不是0呀,C的是-0.25, H的是0.1....
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Itrustyouarekeepingwell?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianweiwang

铁虫 (小有名气)

z5476t4508(金币+10):xiexie 2010-09-06 13:41:04
加上下面的几行试试
HrmStr1  N    CA    490.0   1.335
HrmStr1  C    H2    490.0   1.335
HrmStr1  CA   CA    490.0   1.335
HrmBnd1  CA  N   CA    70.0  116.000
HrmBnd1  CA  N   C     70.0  116.000
HrmBnd1  CA  CA  N     70.0  116.000
HrmBnd1  CA  N   H2    70.0  116.000
HrmBnd1  N   CA  H4     70.0  116.000
HrmBnd1  N   C   H2    70.0  116.000
HrmBnd1  O   C   H2    70.0  116.000
2楼2010-08-31 17:09:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dragonli

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请解释什么意思!谢谢
3楼2010-09-13 14:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 z5476t4508 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 国社科又开始会评了,不知道这次命运如何 +6 雨打竹帘 2026-08-30 10/500 2026-08-31 09:29 by huixian257
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 中青基了要发朋友圈吗? +6 349506619 2026-08-28 6/300 2026-08-31 08:43 by kingkocxr
[基金申请] 面上意见出来了 +10 黄鸟于飞Chao 2026-08-29 19/950 2026-08-31 08:01 by blueearth171
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 怎么看青基中了没有啊 +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 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合现在查不到了吗? +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
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见