计算化学公社

标题: acpype命令出现错误 [打印本页]

作者
Author:
黄舒伟    时间: 2021-4-9 10:12
标题: acpype命令出现错误
请问一下,我按照社长说的http://bbs.keinsci.com/thread-428-1-1.html 用acpype脚本去生成拓扑,用http://svn.code.sf.net/p/ccpn/co ... python/acpype/test/ 中的例子AAA.mol2去复现的时候,出现了报错如下,请问是什么情况呢?(已经chmod 777 ./acpype.py)似乎是引用Antechamber出错了,请问应该怎么解决?(Antechamber在我服务器单独是可以使用的)
  1. (base) [learning@rtx2 2text]$ ./acpype.py -i AAA.mol2
  2. ============================================================================
  3. | ACPYPE: AnteChamber PYthon Parser interfacE v. 0 0 Rev: 0 (c) 2021 AWSdS |
  4. ============================================================================
  5. WARNING: no 'babel' executable, no PDB file as input can be used!
  6. ==> ... charge set to 0
  7. WARNING: Old version of antechamber. Strongly consider upgrading to AmberTools
  8. WARNING: Setting mopac for antechamber
  9. ==> Executing Antechamber...
  10. ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

  11. Welcome to antechamber 19.0: molecular input file processor.

  12. acdoctor mode is on: check and diagnose problems in the input file.
  13. -- Check Format for mol2 File --
  14.    Status: pass
  15. Info: Finished reading file (AAA.mol2); atoms read (33), bonds read (32).
  16. -- Check Unusual Elements --
  17.    Status: pass
  18. -- Check Open Valences --
  19.    Status: pass
  20. -- Check Geometry --
  21.       for those bonded   
  22.       for those not bonded   
  23.    Status: pass
  24. -- Check Weird Bonds --
  25.    Status: pass
  26. -- Check Number of Units --
  27.    Status: pass
  28. acdoctor mode has completed checking the input file.

  29. Running: /home/alvin/apps/anaconda3/bin/bondtype -j full -i ANTECHAMBER_BOND_TYPE.AC0 -o ANTECHAMBER_BOND_TYPE.AC -f ac

  30. ---Judge bond type for Residue 1 with ID of 1 and Name of ALA ---

  31. ---Judge bond type for Residue 2 with ID of 2 and Name of ALA ---

  32. ---Judge bond type for Residue 3 with ID of 3 and Name of ALA ---

  33. Running: /home/alvin/apps/anaconda3/bin/atomtype -i ANTECHAMBER_AC.AC0 -o ANTECHAMBER_AC.AC -p gaff
  34. Info: Total number of electrons: 124; net charge: 0

  35. Running: /home/alvin/apps/anaconda3/bin/mopac.sh
  36. sh: /home/alvin/apps/anaconda3/bin/mopac.sh: No such file or directory
  37. /home/alvin/apps/anaconda3/bin/to_be_dispatched/antechamber: Fatal Error!
  38. Cannot properly run "/home/alvin/apps/anaconda3/bin/mopac.sh".

  39. ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  40. ERROR: Antechamber failed
  41. ++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  42. Cannot open file (AAA_bcc_gaff.mol2) with mode (r).
  43. No such file or directory
  44. ++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  45. ERROR: Parmchk failed
  46. ERROR: Tleap failed
  47. ==> Removing temporary files...
  48. ACPYPE FAILED: [Errno 2] No such file or directory: 'AAA_AC.prmtop'
  49. Total time of execution: less than a second
复制代码



作者
Author:
sobereva    时间: 2021-4-10 01:32
看上去是调用MOPAC出错
要求acpype不做几何优化再试 (, 下载次数 Times of downloads: 44)

作者
Author:
黄舒伟    时间: 2021-4-10 10:01
社长,永远的神!
作者
Author:
黄舒伟    时间: 2021-4-10 16:54
本帖最后由 黄舒伟 于 2021-4-11 13:19 编辑

社长,我想请问一下,我按您说的更改了文件参数后,可以正确生成AAA.mol2的拓扑,但是对于我想要生成的H.mol2的拓扑却出现了如下报错(在您说的在线的acpype网站上能正确H.mol2的拓扑):
  1. (base) [learning@rtx2 text]$ ./acpype.py -i H.mol2 -c user
  2. ============================================================================
  3. | ACPYPE: AnteChamber PYthon Parser interfacE v. 0 0 Rev: 0 (c) 2021 AWSdS |
  4. ============================================================================
  5. WARNING: no 'babel' executable, no PDB file as input can be used!
  6. ACPYPE FAILED: [Errno 2] No such file or directory: 'tmp'
  7. Total time of execution: less than a second
复制代码
想请问一下为什么会出现这种情况呢?

作者
Author:
黄舒伟    时间: 2021-4-23 11:22
黄舒伟 发表于 2021-4-10 16:54
社长,我想请问一下,我按您说的更改了文件参数后,可以正确生成AAA.mol2的拓扑,但是对于我想要生成的H.mo ...

查阅了前人的经验,应该是ambertools的版本有问题,不过也没去安装最新的去试,比较在线的就可以的话那本地是否可以也就无所谓了。




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3