计算化学公社

标题: MM 调用GROMACS力场报错 [打印本页]

作者
Author: zifei 时间: 2021-8-19 21:04
标题: MM 调用GROMACS力场报错
我用acpype生成的GROMACS和AMBER的力场，用CP2K调用。AMBER的可以正常运行，GROMACS会报“ Caught signal 11 (Segmentation fault: address not mapped to object at address 0x1)”的错误，请问是什么原因呢？体系很简单，就是一个水分子测试的。CP2K的输入文件和GROMACS拓扑文件在附件中。
&GLOBAL
PROJECT test
PRINT_LEVEL LOW
RUN_TYPE ENERGY
&END GLOBAL

&FORCE_EVAL
  METHOD FIST
  &MM
&FORCEFIELD
      PARMTYPE G96
      PARM_FILE_NAME h2o_GMX.top
      &SPLINE
         EMAX_SPLINE 100000
         RCUT_NB 10.0
      &END
&END FORCEFIELD

&POISSON             ! solver for non periodic calculations
   &EWALD
      EWALD_TYPE NONE
   &END EWALD
&END POISSON
  &END MM

  &SUBSYS                ! system description
&CELL
   ABC  20. 20. 20.
&END CELL
&TOPOLOGY
      CONN_FILE_FORMAT G96
      CONN_FILE_NAME h2o_GMX.top
      COORD_FILE_FORMAT xyz
      COORD_FILE_NAME h2o.xyz
&END TOPOLOGY
&END SUBSYS
&END FORCE_EVAL

作者
Author: zifei 时间: 2021-8-20 18:26
找到原因了。。G96是gromos的格式，不是gromacs。。。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)