计算化学公社

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

[GROMACS] 使用网站acpype生成拓扑文件出现以下错误

[复制链接 Copy URL]

4

帖子

0

威望

13

eV
积分
17

Level 1 能力者

跳转到指定楼层 Go to specific reply
楼主
网站:https://www.bio2byte.be/acpype/

生成拓扑文件时出现错误不知道如何解决,希望有大佬指教!

Logs for bQ3vRfx56lYRyGcIskHq | n20.pdb | Sat May 04 14:15 2024Details
Hash ID: bQ3vRfx56lYRyGcIskHq
Submission date: May 4, 2024, 2:15 p.m.
Submission filename: n20.pdb
Job status: Failed
Input type: INPUT_FILE
Smiles?: None
Charge type: bcc
Net charge: auto
Multiplicity: 1
Atom type: gaff

Log file===========================================================================| ACPYPE: AnteChamber PYthon Parser interfacE v. 2022.7.21 (c) 2024 AWSdS |===========================================================================DEBUG: Python Version 3.7.7 (default, Mar 10 2020, 15:16:38) [GCC 7.5.0]DEBUG: Max execution time tolerance is 1hDEBUG: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/antechamber -dr no -i n20.pdb.mol2 -fi mol2 -o tmp -fo ac -pf yDEBUG: Welcome to antechamber 21.0: molecular input file processor.The atom type is set to gaff; the options available to the -at flag are gaff, gaff2, amber, bcc, and sybyl.Info: Bond types are assigned for valence state (2) with penalty (1).DEBUG: setResNameCheckCoords doneWARNING: no charge value given, trying to guess one...DEBUG: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/antechamber -dr no -i n20.pdb.mol2 -fi mol2 -o tmp -fo mol2 -c gas -pf nDEBUG: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/antechamber -dr no -i tmp -fi mol2 -o tmp -fo mol2 -c wc -cf tmp.crg -pf nDEBUG: ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++DEBUG: Welcome to antechamber 21.0: molecular input file processor.The atom type is set to gaff; the options available to the -at flag are gaff, gaff2, amber, bcc, and sybyl.Info: Bond types are assigned for valence state (2) with penalty (1).DEBUG: ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++DEBUG: readMol2TotalCharge: -1.0DEBUG: Net charge drift '0.000000'==> ... charge set to -1==> Executing Antechamber...DEBUG: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/antechamber -dr no -i n20.pdb.mol2 -fi mol2 -o n20.pdb.mol2_bcc_gaff.mol2 -fo mol2 -c bcc -nc -1 -m 1 -s 2 -df 2 -at gaff -pf n DEBUG: No file left to be removedERROR: ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ERROR: Welcome to antechamber 21.0: molecular input file processor.The atom type is set to gaff; the options available to the -at flag are gaff, gaff2, amber, bcc, and sybyl.Info: Finished reading file (n20.pdb.mol2); atoms read (31), bonds read (32).Info: Determining atomic numbers from atomic symbols which are case sensitive.Running: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/bondtype -j full -i ANTECHAMBER_BOND_TYPE.AC0 -o ANTECHAMBER_BOND_TYPE.AC -f acInfo: Bond types are assigned for valence state (2) with penalty (1).Running: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/atomtype -i ANTECHAMBER_AC.AC0 -o ANTECHAMBER_AC.AC -p gaffInfo: Total number of electrons: 171; net charge: -1Info: The number of electrons is odd (171). Please check the total charge (-nc flag) and spin multiplicity (-m flag).Running: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/sqm -O -i sqm.in -o sqm.out/home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/wrapped_progs/antechamber: Fatal Error!Cannot properly run "/home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/sqm -O -i sqm.in -o sqm.out".ERROR: ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ERROR: Antechamber failedDEBUG: /home/b2b-user/web/lib/python3.7/site-packages/acpype/amber21-11_linux/bin/parmchk2 -i n20.pdb.mol2_bcc_gaff.mol2 -f mol2 -o n20.pdb.mol2_AC.frcmodERROR: ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ERROR: Cannot open file (n20.pdb.mol2_bcc_gaff.mol2) with mode (r).No such file or directoryERROR: ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ERROR: Parmchk failedERROR: Tleap failedACPYPE FAILED: [Errno 2] No such file or directory: 'n20.pdb.mol2_AC.prmtop'  File "/home/b2b-user/web/lib/python3.7/site-packages/acpype/acs_api.py", line 129, in acpype_api    molecule.createMolTopol()  File "/home/b2b-user/web/lib/python3.7/site-packages/acpype/topol.py", line 1116, in createMolTopol    self.topFileData = open(self.acTopFileName).readlines()Total time of execution: less than a secondServer messages
  • Waiting for initializing the ACPYPE process. [Attempt 1/2] Error ACPYPE: ERROR: [Errno 2] No such file or directory: 'n20.pdb.mol2_AC.prmtop'. [Attempt 2/2] Error ACPYPE: ERROR: [Errno 2] No such file or directory: 'n20.pdb.mol2_AC.prmtop'.

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

2#
发表于 Post on 2024-5-23 15:06:27 | 只看该作者 Only view this author
用sobtop就完了,省事简单快速,还没那么多破毛病。acpype完全可以弃了
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

4

帖子

0

威望

13

eV
积分
17

Level 1 能力者

3#
 楼主 Author| 发表于 Post on 2024-5-23 16:51:06 | 只看该作者 Only view this author
sobereva 发表于 2024-5-23 15:06
用sobtop就完了,省事简单快速,还没那么多破毛病。acpype完全可以弃了
http://sobereva.com/soft/Sobtop

感谢你的回答!我下载了sobtop,并依据官网步骤进行操作,想要生成GAFF力场的小分子拓扑文件。在运行中显示操作完成,但拓扑文件我却不知道在哪里获得,我是在压缩包的examples文件夹内新建一个文件夹放置我的小分子,但在相应文件夹并没有发现应生成的拓扑文件,请问我是在操作上有什么错误吗?求指正!

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

4#
发表于 Post on 2024-5-24 04:33:27 | 只看该作者 Only view this author
mnmn 发表于 2024-5-23 16:51
感谢你的回答!我下载了sobtop,并依据官网步骤进行操作,想要生成GAFF力场的小分子拓扑文件。在运行中显 ...

就在当前目录下产生
搞清楚什么叫当前目录,这是计算机最基本常识。linux下运行pwd命令显示的目录,以及windows下双击运行sobtop.exe时sobtop.exe所在的目录
北京科音自然科学研究中心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

4

帖子

0

威望

13

eV
积分
17

Level 1 能力者

5#
 楼主 Author| 发表于 Post on 2024-5-24 18:28:33 | 只看该作者 Only view this author
本帖最后由 mnmn 于 2024-5-24 18:31 编辑
sobereva 发表于 2024-5-24 04:33
就在当前目录下产生
搞清楚什么叫当前目录,这是计算机最基本常识。linux下运行pwd命令显示的目录,以及 ...

找到了,感谢

4

帖子

0

威望

13

eV
积分
17

Level 1 能力者

6#
 楼主 Author| 发表于 Post on 2024-5-24 18:33:17 | 只看该作者 Only view this author
sobereva 发表于 2024-5-24 04:33
就在当前目录下产生
搞清楚什么叫当前目录,这是计算机最基本常识。linux下运行pwd命令显示的目录,以及 ...

另外,我话有一个疑问。
3 Assign GAFF atom types as much as possible
当我选择以上GAFF力场,后出现的以下选项:
Well, LJ parameters of all atoms can be found in LJ_param.dat. Now we go to next stage

How to generate bonded parameters for the topology file?
Note: "prebuilt" parameters in this context refer to those in bonded_param.dat
0 Return
1 Without bonded parameters
2 All bonded parameters are determined by mSeminario method
3 Prebuilt bonded parameters if possible, missing ones are supplemented by mSeminario method
4 Same as option 3, but missing ones are arbitrarily guessed
5 Bonded parameters related to selected atoms (usually rigid region) are determined by mSeminario method, others use prebuilt bonded parameters, then missing ones are supplemented by mSeminario method
6 Same as option 5, but missing ones are arbitrarily guessed
7 All bond and angle parameters are determined by mSeminario method, others use prebuilt bonded parameters, then missing ones are supplemented by mSeminario method

我应该如何选择?

6万

帖子

99

威望

6万

eV
积分
125128

管理员

公社社长

7#
发表于 Post on 2024-5-25 04:09:48 | 只看该作者 Only view this author
mnmn 发表于 2024-5-24 18:33
另外,我话有一个疑问。
3 Assign GAFF atom types as much as possible
当我选择以上GAFF力场,后出现 ...

把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

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 03:05 , Processed in 0.195435 second(s), 26 queries , Gzip On.

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