计算化学公社

标题: 求助 gromacs grompp步骤出错 [打印本页]

作者
Author:
lqhenwunai    时间: 2020-8-28 17:06
标题: 求助 gromacs grompp步骤出错
我用packmol生成300个水分子,box尺寸为30*30*30


生成h2o-out.pdb 文件,部分内容如下:
HEADER
TITLE     Built with Packmol
REMARK   Packmol generated pdb file
REMARK   Home-Page: http://m3g.iqm.unicamp.br/packmol
REMARK
HETATM    1  O   SOL A   1      -3.286  -9.267  -0.230  1.00  0.00           O
HETATM    2  H   SOL A9999      -2.832 -10.000   0.215  1.00  0.00           H
HETATM    3  H   SOL A9999      -4.189  -9.187   0.119  1.00  0.00           H



然后我在top里写道:
#include "oplsaa.ff/forcefield.itp"
#include "oplsaa.ff/tip3p.itp"

[ System ]
TIP3P in water

[ Molecules ]
SOL               300

mdp文件如下:
integrator  = steep         ; Algorithm (steep = steepest descent minimization)
emtol       = 1000.0        ; Stop minimization when the maximum force < 1000.0 kJ/mol/nm
emstep      = 0.01          ; Minimization step size
nsteps      = 200000         ; Maximum number of (minimization) steps to perform

; Parameters describing how to find the neighbors of each atom and how to calculate the interactions
nstlist         = 1         ; Frequency to update the neighbor list and long range forces
cutoff-scheme   = Verlet    ; Buffered neighbor searching
ns_type         = grid      ; Method to determine neighbor list (simple, grid)
coulombtype     = PME       ; Treatment of long range electrostatic interactions
rcoulomb        = 1.0       ; Short-range electrostatic cut-off
rvdw            = 1.0       ; Short-range Van der Waals cut-off
pbc             = xyz       ; Periodic Boundary Conditions in all 3 dimensions



但是当我gmx grompp -f minim.mdp -o energy-minimization.tpr -c h2o-out.pdb时,我得到以下报错信息:
ERROR 1 [file topol.top, line 12]:
  ERROR: The cut-off length is longer than half the shortest box vector or
  longer than the smallest box diagonal element. Increase the box size or
  decrease rlist.

请问这应该怎么修改呢?我之前是用gmx solvate  -box 2.3 2.3 2.3 -o input.gro -p topol.top来生成水盒子,用的是同样的mdp文件,没有问题。但是现在我想从packmol生成的pdb文件来做初始构型,于是就有以上报错信息。请问这是怎么回事呢?

谢谢大家的帮助。



作者
Author:
sobereva    时间: 2020-8-29 00:34
盒子太小,提示已经说明了原因

这年头用300个水来模拟水的性质也太次了,最起码500个
用packmol构建的水盒子的初始结构明显不如用gmx solvate

作者
Author:
lqhenwunai    时间: 2020-8-29 13:41
sobereva 发表于 2020-8-29 00:34
盒子太小,提示已经说明了原因

这年头用300个水来模拟水的性质也太次了,最起码500个

谢谢,但是我之前用gmx solvate  -box 2.3 2.3 2.3 -o input.gro -p topol.top来生成水盒子,用的是同样的mdp文件,没有问题。

但是现在增大盒子尺寸反而有问题了。
作者
Author:
sobereva    时间: 2020-8-30 14:01
lqhenwunai 发表于 2020-8-29 13:41
谢谢,但是我之前用gmx solvate  -box 2.3 2.3 2.3 -o input.gro -p topol.top来生成水盒子,用的是同样 ...

要比就用solvate的情况比
如果你不设add_box_sides,packmol根本不产生盒子信息





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