24小时热门版块排行榜    

查看: 2722  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

倩430

铁虫 (小有名气)

[求助] 高斯输出文件转化为mol2文件时出现问题 -c resp 与 -c bcc

将高斯输出文件转化为mol2文件时,使用命令如下
antechamber -i a.log -fi gout -o a.mol2 -fo mol2 -c resp
得到错误提示
unit   10 Error on open: ANTECHAMBER.ESP

unit     3 Error on open:qout

cannot open charge file to read:QOUT ,exit
而使用-c bcc时候就没有错误提示。请问这是什么原因导致的。两种方法分别适合什么时候用呢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倩430

铁虫 (小有名气)

引用回帖:
5楼: Originally posted by lujunyan1118 at 2012-10-30 19:51:42
#B3LYP/6-31G* integral(Grid=UltraFine) pop(mk,readradii) scf=xqc iop(6/33=2)

我用的是这个,其中最主要的是iop(6/33=2), 有它才能用log文件来拟合RESP

明白些了呵呵,就是说把那个加上就应该差不多了吧。我用antechamber生成的输入文件如下,貌似比我用gv,关键词多了好多呢
--Link1--
%chk=molecule
#HF/6-31G* SCF=tight Test Pop=MK iop(6/33=2) iop(6/42=6) opt
6楼2012-10-31 15:17:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

wangyan10

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
chaizhm: 金币+2, 谢谢~ 2012-10-30 08:54:55
-c后面的bcc和resp是两种不同的charge model。resp采用的方法是HF,而bcc采用AM1方法,这两者是不同的。http://ambermd.org/tutorials/adv ... resp_paper_jacs.pdf   这是一篇介绍resp charge model的原始文献,想了解的话可以看一下。至于哪种charge model适用于什么样的体系,自己做科研去领悟吧!这个不是一会半会能说明白的问题。
努力的生活着
2楼2012-10-29 20:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
chaizhm: 金币+3, 谢谢~ 2012-10-30 08:55:05
你算高斯的文件使用antechamber生成的吗?
resp拟合的时候需要用到高斯log文件里面写的esp信息。要在gjf的头上加特殊关键字。如果gjf文件使用antechamber生成的,应该没有问题

一般用amber力场的话就用resp电荷,amber立场参数是这么来的。

bcc我不怎么用,bcc要快一些,不需要高斯算。可能做对接的时候,处理大批量小分子用用
3楼2012-10-29 20:05:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

倩430

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by lujunyan1118 at 2012-10-29 20:05:05
你算高斯的文件使用antechamber生成的吗?
resp拟合的时候需要用到高斯log文件里面写的esp信息。要在gjf的头上加特殊关键字。如果gjf文件使用antechamber生成的,应该没有问题

一般用amber力场的话就用resp电荷 ...

高斯输入文件不是antechamber 生成的,用gassianview生成的.gif文件。自己修改部分后,内容如下
%chk=egcg.chk
%nprocshared=8
%mem=8000MB
# opt freq rb3lyp/6-31g

egcg.mol2

0 1
C                  0.18226000    1.41376000   -0.12743000
C                  1.37057000    2.14000000    0.07324000请问:如果用antechamber生成.gjf的话,生成之后还需要修改哪里(就是你说的关键词是哪个,难道改rb3lyp为HF吗)才能保证-resp这个命令能成功运行。
4楼2012-10-30 17:26:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +3 基诺咪客 2026-08-17 3/150 2026-08-17 15:25 by xinshengwu
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] filecode=后面第一个是大写字母 +7 wangze12014 2026-08-14 9/450 2026-08-17 14:59 by wangze12014
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +33 医学老男孩 2026-08-13 80/4000 2026-08-17 14:40 by yq199112088
[基金申请] 感觉是下周放榜了 +3 angus9576 2026-08-17 6/300 2026-08-17 14:32 by 启萌科技
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见