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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】80金币求助oniom的问题，需要您的解答（只要回复就给金币啦）{金币发完了}
北京石油化工学院2026年研究生招生接收调剂公告
221/3返回列表123下一页
查看: 1695  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主

[交流] 【求助】80金币求助oniom的问题，需要您的解答（只要回复就给金币啦）{金币发完了}

半个月没人回复了，决定放弃。当时悬赏的80金币，现在还剩75金币没发出去。现在，只要回复就给金币！


80金币求助oniom的问题。

我一共不到90金币，没办法，这80金币已经是我能给出的最后一个铜板了。如果某大侠知道答案，但是觉得80金币太少，我可以再用Q币继续补偿您，呵呵。

我会用的分层算，一直是Gaussian98格式的oniom，也就是那种极其简单的分层算法。现在看Gaussian03格式的分层算法，晕了。

先给一个具体的分子，酞菁锂，结构如图1（PcLi2）。我分3层了。第一层是那两个Li原子，第二层是与Li直接相连的四个N原子，第三层是分子的其余部分。
我在GaussView4.1中定义各层，如图2.
图片有些大，请稍等。

我用GaussView4.1直接生成的文件如下：
%nprocshared=1
%mem=500MB
%chk=G:\ByHand\oniom\Li2PcONIOM.chk
# opt freq=noraman oniom(mp2/6-311+g(3df,3pd):b3pw91/6-31g(d):b3lyp/3-21g)

Li2Pc Oniom

0 1 -2 1 -2 1 2 1 2 1 2 1
C                0   -1.06827940   -2.59205361    0.00000000 L H 5
C                0   -0.64817712   -4.05208072    0.00000000 L
C                0    0.64817712   -4.05208072    0.00000000 L
C                0    1.06827940   -2.59205361    0.00000000 L H 5
N                0    0.00000000   -1.69380459    0.00000000 M
C                0   -0.64817712    4.05208072    0.00000000 L
C                0    0.64817712    4.05208072    0.00000000 L
C                0   -1.06827940    2.59205361    0.00000000 L H 10
C                0    1.06827940    2.59205361    0.00000000 L H 10
N                0    0.00000000    1.69380459    0.00000000 M
C                0    2.59205361   -1.06827940    0.00000000 L H 13
C                0    4.05208072   -0.64817712    0.00000000 L
N                0    1.69380459    0.00000000    0.00000000 M
C                0    4.05208072    0.64817712    0.00000000 L
C                0    2.59205361    1.06827940    0.00000000 L H 13
C                0   -2.59205361   -1.06827940    0.00000000 L H 18
C                0   -4.05208072   -0.64817712    0.00000000 L
N                0   -1.69380459    0.00000000    0.00000000 M
C                0   -4.05208072    0.64817712    0.00000000 L
C                0   -2.59205361    1.06827940    0.00000000 L H 18
N                0    2.29712177    2.29712177    0.00000000 L
N                0   -2.29712177    2.29712177    0.00000000 L
N                0   -2.29712177   -2.29712177    0.00000000 L
N                0    2.29712177   -2.29712177    0.00000000 L
C                0    5.23533458    1.38094601    0.00000000 L
H                0    5.20203797    2.45044856    0.00000000 L
C                0    6.46048480    0.68548141    0.00000000 L
H                0    7.38597263    1.22242753    0.00000000 L
C                0    6.46048480   -0.68548141    0.00000000 L
H                0    7.38597263   -1.22242753    0.00000000 L
C                0    5.23533458   -1.38094601    0.00000000 L
H                0    5.20203797   -2.45044856    0.00000000 L
C                0    1.38094601   -5.23533458    0.00000000 L
H                0    2.45044856   -5.20203797    0.00000000 L
C                0    0.68548141   -6.46048480    0.00000000 L
H                0    1.22242753   -7.38597263    0.00000000 L
C                0   -0.68548141   -6.46048480    0.00000000 L
H                0   -1.22242753   -7.38597263    0.00000000 L
C                0   -1.38094601   -5.23533458    0.00000000 L
H                0   -2.45044856   -5.20203797    0.00000000 L
C                0   -5.23533458   -1.38094601    0.00000000 L
H                0   -5.20203797   -2.45044856    0.00000000 L
C                0   -6.46048480   -0.68548141    0.00000000 L
H                0   -7.38597263   -1.22242753    0.00000000 L
C                0   -6.46048480    0.68548141    0.00000000 L
H                0   -7.38597263    1.22242753    0.00000000 L
C                0   -5.23533458    1.38094601    0.00000000 L
H                0   -5.20203797    2.45044856    0.00000000 L
C                0   -1.38094601    5.23533458    0.00000000 L
H                0   -2.45044856    5.20203797    0.00000000 L
C                0   -0.68548141    6.46048480    0.00000000 L
H                0   -1.22242753    7.38597263    0.00000000 L
C                0    0.68548141    6.46048480    0.00000000 L
H                0    1.22242753    7.38597263    0.00000000 L
C                0    1.38094601    5.23533458    0.00000000 L
H                0    2.45044856    5.20203797    0.00000000 L
Li               0    0.00000000    0.00000000    1.15363497 H
Li               0    0.00000000    0.00000000   -1.15363497 H

我的问题如下：
（1）关于电荷和自旋多重度这一行：0 1 -2 1 -2 1 2 1 2 1 2 1
请问0 1是不是整个分子的电荷和自旋多重度？
后面的-2 1 -2 1 2 1 2 1 2 1，我知道是五对“电荷和自旋多重度”，请问它们分别对应的是什么？
（2）关于分子说明，如第一行分子说明：
C                0   -1.06827940   -2.59205361    0.00000000 L H 5
请问“H 5”是什么意思？是不是表示它和H层之间仅仅隔着一个5号原子？如果H层和L层之间隔着不止一个原子，是不是就不会有“H 5”这样的声明？还有，当H层和L层之间仅有一个原子的时候，请问为什么要有这样的声明？
（3）请问H层和L层能不能直接连接？H层和L层之间必须被M层隔开吗？
（4）在Gaussian03的说明书中有这么一段：
Atom layer assignment is done as part of the molecule specification, via additional parameters on each line according to the following syntax:
atom coordinate-spec layer [link-atom [bonded-to [scale-fac1 [scale-fac2 [scale-fac3]]]]]
where atom and coordinate-spec represent the normal molecule specification input for the atom. Layer is a keyword indicating the layer assignment for the atom, one of High, Medium and Low. The other optional parameters specify how the atoms located at a layer boundary are to be treated. You use link-atom to specify the atom with which to replace the current atom (it can include atom type and partial charge and other parameters). Link atoms are necessary when covalent bonding exists between atoms in different layers in order to saturate the (otherwise) dangling bonds.
请问：为什么要存在一种Link atom（链接原子，我不知道我这样翻译对不对）？它们的存在是为了补全电荷还是什么别的意图？在规定各个连接原子完毕之后，这些链接原子是否也要算进各个层的电荷和自旋多重度中去？
（5）[link-atom [bonded-to [scale-fac1 [scale-fac2 [scale-fac3]]]]]
请问换算因子scale-fac1、scale-fac2、scale-fac3分别是0.0、0.0、0.0时候，oniom在各层之间使用常规换算因子。请问如果不用常规换算因子，那么我是不是直接写上这个因子的数值？

我一共不到90金币，没办法，这80金币已经是我能给出的最后一个铜板了。如果某大侠知道答案，但是觉得80金币太少，我可以再用Q币继续补偿您，先谢谢了。

[ Last edited by yjcmwgk on 2009-7-2 at 20:56 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-06-06 17:04:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
无奈了。没有人回答吗？唉……等待达人回答ing

[ Last edited by yjcmwgk on 2009-6-6 at 21:43 ]
赞一下 回复此楼
2楼2009-06-06 20:53:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)
★ ★ ★ ★
gwdavid(金币+3,VIP+0):3q 6-8 18:45
yjcmwgk(金币+1,VIP+0):只要回复就给金币啦！ 7-2 20:17
建议把chemdraw那个图换成gjf格式,并标出lablel

楼主别急,可能还没有人看到,我帮你加成红题,醒目一下.

另外有个疑问,这样选三层,怎么保证每层就是你想设置的呢?
赞一下(25人) 回复此楼
3楼2009-06-06 21:13:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
引用回帖:
Originally posted by erylingjet at 2009-6-6 21:13:
建议把chemdraw那个图换成gjf格式,并标出lablel

楼主别急,可能还没有人看到,我帮你加成红题,醒目一下.

另外有个疑问,这样选三层,怎么保证每层就是你想设置的呢?

呵呵，谢谢斑竹！
建议把chemdraw那个图换成gjf格式,并标出lablel----我已经把gjf文件完整的贴出来了，就在正文里，只要复制进gjf文件，就能打开。
这样选三层,怎么保证每层就是你想设置的呢?----我用GV4.1的Layer工具，直接点击各个原子，用这种方式保证我想设置的HLM层分别就是我要的那些原子
赞一下 回复此楼
4楼2009-06-06 21:41:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wuhanhgf2002

金虫 (正式写手)
★ ★
erylingjet(金币+1,VIP+0):感谢交流 6-8 10:18
yjcmwgk(金币+1,VIP+0):只要回复就给金币啦！ 7-2 20:17
你好，你的输入文件可能有错误，你手上有没有运行成功的oniom例子？如果没有的话，我可以给你一个。
赞一下(25人) 回复此楼
5楼2009-06-07 13:42:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
引用回帖:
Originally posted by wuhanhgf2002 at 2009-6-7 13:42:
你好，你的输入文件可能有错误，你手上有没有运行成功的oniom例子？如果没有的话，我可以给你一个。

十分感谢，楼上大侠能给我一个运行成功的oniom例子吗？
qdd-china@foxmail.com
赞一下 回复此楼
6楼2009-06-08 09:27:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

优秀版主
★
gwdavid(金币+1,VIP+0):3q 6-8 18:45
引用回帖:
Originally posted by wuhanhgf2002 at 2009-6-7 13:42:
你好，你的输入文件可能有错误，你手上有没有运行成功的oniom例子？如果没有的话，我可以给你一个。

我有运行成功的oniom，但是仅限于如下样式的简单的oniom：

# oniom(方法/基组:方法/基组:方法/基组) ……
自旋多重度只定义整个分子的，没有定义各层的自旋多重度
赞一下(10人) 回复此楼
7楼2009-06-08 09:33:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

magic6

金虫 (正式写手)
★
yjcmwgk(金币+1,VIP+0):只要回复就给金币啦！ 7-2 20:17
顶一下，期待高手解答～
回复此楼
8楼2009-06-12 17:19:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

qxh8419

金虫 (正式写手)
★
yjcmwgk(金币+1,VIP+0):只要回复就给金币啦！ 7-2 20:17
希望大家帮忙，虫友们都是好样的
赞一下(14人) 回复此楼
博采众长
9楼2009-06-13 21:09:38
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

niliu

铁杆木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金币+3,VIP+0):谢谢 6-14 09:48
yjcmwgk(金币+5,VIP+0):谢谢。其他问题请问能回答吗？ 6-15 19:22
yjcmwgk(金币+1,VIP+0):只要回复就给金币啦！ 7-2 20:17
引用回帖:
Originally posted by yjcmwgk at 2009-6-6 17:04:
（2）关于分子说明，如第一行分子说明：
C                0   -1.06827940   -2.59205361    0.00000000 L H 5
请问“H 5”是什么意思？是不是表示它和H层之间仅仅隔着一个5号原子？如果H层和L层之间隔着不止一个原子，是不是就不会有“H 5”这样的声明？还有，当H层和L层之间仅有一个原子的时候，请问为什么要有这样的声明？

H  表示高一层计算时这个C原子当着H原子来饱和悬挂键（dangling bond）
5  表示高一层计算时这个C原子与编号为5的原子相连
赞一下(40人) 回复此楼
10楼2009-06-13 23:09:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yjcmwgk 的主题更新
221/3返回列表123下一页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[考研] 0856调剂 +8 曲听筠 2026-03-30 8/400 2026-04-04 08:46 by tianyyysss
[考研] 265求调剂 +20 梁梁校校 2026-04-01 21/1050 2026-04-04 00:38 by userper
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 08工科，295，接受跨专业调剂 +8 lmnlzy 2026-03-30 8/400 2026-04-03 13:08 by nalakaiqi
[考研] 298求调剂 +5 种圣赐 2026-03-28 5/250 2026-04-03 12:26 by xingguangj
[考研] 320求调剂 +3 农业工程与信息� 2026-04-03 3/150 2026-04-03 11:40 by 土木硕士招生
[考研] 306分材料与化工求调剂 +4 黎吧啦啦你很有� 2026-04-03 4/200 2026-04-03 11:24 by wangjy2002
[考研] 求调剂！生物与医药专硕 +4 逆转陆先生 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[考研] 372分材料与化工（085600）一志愿湖南大学求调剂 +5 蓝笺片 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[考研] 311求调剂 +14 蓝月亮亮 2026-03-30 14/700 2026-04-02 12:18 by 1753564080
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 考研调剂0702 +3 梅子味晚霞 2026-03-28 3/150 2026-04-02 10:10 by jp9609
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 江苏苏北高校诚邀调剂同学 +3 zzll406 2026-03-31 3/150 2026-03-31 16:54 by 及时行乐fan
[考研] 生物考研337分求调剂 +4 cgxin 2026-03-30 6/300 2026-03-31 14:18 by 记事本2026
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 0703化学321分求调剂 +10 三dd. 2026-03-30 11/550 2026-03-30 19:24 by markhwc
信息提示
关闭
请填处理意见
关闭