计算化学公社

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

[GROMACS] 求助:acpype生成拓扑文件时报错,提示结构间距太大

[复制链接 Copy URL]

23

帖子

0

威望

49

eV
积分
72

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
以下为运行命令:
acpype -i ligand.mol2 -a amber2 -c user -d
报错如下:
[root@192 CXCL8_5d]# acpype -i ligand.mol2 -a amber2 -c user -d
============================================================================
| ACPYPE: AnteChamber PYthon Parser interfacE v. 2023.10.27 (c) 2024 AWSdS |
============================================================================
DEBUG: CLI: -i ligand.mol2 -a amber2 -c user -d
DEBUG: Python Version 3.11.7 | packaged by conda-forge | (main, Dec 23 2023, 14:43:09) [GCC 12.3.0]
DEBUG: Max execution time tolerance is 3h
DEBUG: /home/moi/miniconda3/bin/antechamber -dr no -i ligand.mol2 -fi mol2 -o tmp -fo ac -pf y
DEBUG:
Welcome to antechamber 22.0: molecular input file processor.

Info: The atom type is set to gaff; the options available to the -at flag are
gaff, gaff2, amber, bcc, and sybyl.

Warning: Ignoring Mol2 record type (@<TRIPOS>SET).
Warning: For atom (ID: 41, Name: O6) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 42, Name: O7) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 43, Name: O8) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 44, Name: O9) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 45, Name: O10) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 46, Name: O11) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 47, Name: O12) the best APS is not zero.
Bonds involving this atom are frozen.
Warning: For atom (ID: 48, Name: O13) the best APS is not zero.
Bonds involving this atom are frozen.

ERROR: Atoms TOO scattered (> 3.0 Ang.)
ERROR: ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
ERROR: ['ATOM 45 O10 ']
['ATOM 46 O11 ']
['ATOM 47 O12 ']
['ATOM 48 O13 ']
ERROR: ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
ERROR: Use '-f' option if you want to proceed anyway. Aborting ...
ERROR: Coordinates issues with your system
ERROR: ACPYPE FAILED: Coordinates issues with your system
Traceback (most recent call last):
  File "/home/moi/miniconda3/lib/python3.11/site-packages/acpype/cli.py", line 106, in init_main
    molecule = ACTopol(
               ^^^^^^^^
  File "/home/moi/miniconda3/lib/python3.11/site-packages/acpype/topol.py", line 3325, in __init__
    self.setResNameCheckCoords()
  File "/home/moi/miniconda3/lib/python3.11/site-packages/acpype/topol.py", line 564, in setResNameCheckCoords
    raise Exception(msg)
Exception: Coordinates issues with your system
Total time of execution: less than a second
DEBUG: No tmp folder left to be removed
Log tmp location: /tmp/tmpmo9ti7fz

实在是不懂结构上出了什么问题,能否请大家看看,谢谢!!
结尾附上mol2文件

ligand.mol2

26.08 KB, 下载次数 Times of downloads: 1

mol2文件

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

2#
发表于 Post on 2024-1-26 02:14:05 | 只看该作者 Only view this author
用sobtop(http://sobereva.com/soft/Sobtop)就没这些破毛病,敲几下键盘拓扑文件就出来了,又省事又方便
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

23

帖子

0

威望

49

eV
积分
72

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-1-26 13:12:21 | 只看该作者 Only view this author
sobereva 发表于 2024-1-26 02:14
用sobtop(http://sobereva.com/soft/Sobtop)就没这些破毛病,敲几下键盘拓扑文件就出来了,又省事又方便

请问sob老师,如果我大分子选择的是amber03力场,TIP3P力场,小分子要选择哪个力场呢?

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

4#
发表于 Post on 2024-1-27 03:59:09 | 只看该作者 Only view this author
MOI_001 发表于 2024-1-26 13:12
请问sob老师,如果我大分子选择的是amber03力场,TIP3P力场,小分子要选择哪个力场呢?

AMBER03完全过时了
如今用AMBER力场描述蛋白质,至少也该用AMBER14SB
有机小分子用GAFF
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 20:42 , Processed in 0.160407 second(s), 23 queries , Gzip On.

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