24小时热门版块排行榜    

查看: 2445  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 材料专硕274一志愿陕西师范大学求调剂 +8 薛云鹏 2026-03-13 8/400 2026-03-19 15:36 by haoshis
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[教师之家] 焦虑 +9 水冰月月野兔 2026-03-13 13/650 2026-03-19 09:50 by otani
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
信息提示
请填处理意见