计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3377|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[Amber] 求助,使用AmberTools+Gaussian创建gromacs的拓扑文件出错

[复制链接 Copy URL]

12

帖子

0

威望

67

eV
积分
79

Level 2 能力者

使用Ambertools + Gaussian 09 创建拓扑文件时先用gaussian09计算出了 .out文件,然后使用 antechamber 做了 .prepin文件 并用 parmchk2补充了 frcmod.ala文件在tleap中对alo阴离子的pdb 文件检测中并没有找到键的缺失(图1),但是对于packmol扩包的mix1.pdb文件在tleap中检测却出现了大量的阴离子键角缺失(图2)(阳离子部分没有报错),请问要怎么解决(相关的各类文件如附件所示,其中.lib文件是由.prepin文件在xleap中修改后保存的)

QQ图片20211210112031.png (161.42 KB, 下载次数 Times of downloads: 13)

QQ图片20211210112031.png

QQ图片20211210112044.png (101.76 KB, 下载次数 Times of downloads: 13)

QQ图片20211210112044.png

alo.lib

2.4 KB, 下载次数 Times of downloads: 3

alo.mol2

1.39 KB, 下载次数 Times of downloads: 2

alo.pdb

818 Bytes, 下载次数 Times of downloads: 3

frcmod.ala

3.18 KB, 下载次数 Times of downloads: 2

parmchk2生成的文件

input.mix

227 Bytes, 下载次数 Times of downloads: 2

packmol输入文件

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

2#
发表于 Post on 2021-12-10 13:00:51 | 只看该作者 Only view this author
本帖最后由 ggdh 于 2021-12-10 13:02 编辑

你做frcmod文件的时候产生全部参数试试(加上选项 -a Y)
parmchk2 -a Y  ....

12

帖子

0

威望

67

eV
积分
79

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2021-12-10 19:10:34 | 只看该作者 Only view this author
ggdh 发表于 2021-12-10 13:00
你做frcmod文件的时候产生全部参数试试(加上选项 -a Y)
parmchk2 -a Y  ....

生成的新的frcmod文件和之前的一样...还是有同样的问题

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

4#
发表于 Post on 2021-12-10 19:32:15 | 只看该作者 Only view this author
thetn 发表于 2021-12-10 19:10
生成的新的frcmod文件和之前的一样...还是有同样的问题

搞不定就用ztop。。。或者是acpype把。然后mix的top直接用gmx 的include itp方式,方便多了。

12

帖子

0

威望

67

eV
积分
79

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2021-12-10 19:56:12 | 只看该作者 Only view this author
ggdh 发表于 2021-12-10 19:32
搞不定就用ztop。。。或者是acpype把。然后mix的top直接用gmx 的include itp方式,方便多了。

主要是之前做另外一种的时候也是用的这个方法,当时也没遇到这问题,

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

6#
发表于 Post on 2021-12-11 03:03:23 | 只看该作者 Only view this author
用acpype就好了嘛,或者用在线的生成就好了,不优化结构也不算原子电荷。

完了直接用IQMol的GAFF优化一下存个结构高斯跑个单点Multiwfn算下RESP填进itp里去就好了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 18:32 , Processed in 0.201898 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list