24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2497  |  回复: 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时候就没有错误提示。请问这是什么原因导致的。两种方法分别适合什么时候用呢
回复此楼

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

amber

» 猜你喜欢

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

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

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的回帖

lujunyan1118

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
jiaoyixiong: 金币+2, 鼓励交流 2012-10-30 20:26:18
倩430: 金币+5, ★★★★★最佳答案 2012-10-31 15:18:17
#B3LYP/6-31G* integral(Grid=UltraFine) pop(mk,readradii) scf=xqc iop(6/33=2)

我用的是这个,其中最主要的是iop(6/33=2), 有它才能用log文件来拟合RESP
5楼2012-10-30 19:51:42
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 倩430 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 083200 305分 求二轮调剂 不接受跨专业 +3 Claireyyyy 2026-04-09 3/150 2026-04-09 23:17 by xiaohe123
[考研] 301求调剂 +3 XYPLR 2026-04-05 5/250 2026-04-09 20:55 by lbsjt
[考研] 297求调剂 +22 ljy20040718! 2026-04-03 24/1200 2026-04-09 20:48 by yanenwang
[考研] 291分调剂 +3 上岸小莹加油 2026-04-09 4/200 2026-04-09 17:28 by 1753564080
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 085600材料与化工,求调剂 +11 won_qii 2026-04-07 11/550 2026-04-09 17:03 by luoyongfeng
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +14 Naiko 2026-04-04 14/700 2026-04-09 16:56 by luoyongfeng
[硕博家园] 新一代电子信息294求调剂 不挑学校 +5 Ytyt11 2026-04-09 6/300 2026-04-09 14:40 by Ytyt11
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 328求调剂 +17 lftmya 2026-04-07 18/900 2026-04-09 08:05 by 5268321
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 求调剂 +9 月@163.com 2026-04-07 11/550 2026-04-08 14:48 by qlm5820
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-05 5/250 2026-04-06 15:40 by lin-da
[考研] 一志愿北京化工大学,初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊!
信息提示
请填处理意见