[AMBER] Error when using PB implicit solvent in XMIN in NAB(PB Bomb in pb_init(): Allocation aborted)

From: Guqin Shi <shi.293.osu.edu>
Date: Fri, 9 May 2014 14:44:03 -0400

Hi Amber users,

I was using the XMIN method in NAB package to minimize my
complex(protein+small molecule). The reason I am doing that is I want to
further calculate entropy by using nmode() in NAB.

I was trying to use PB solvent in minimization and below is my input:

======================================================================================
#include "xmin_opt.h"

struct xmin_opt xo;

molecule m;
int natm;
float xyz[dynamic], fret, grad[dynamic];
float energy, grms;

xmin_opt_init(xo);

xo.maxiter = 40000;
xo.grms_tol = 0.0001;
xo.print_level = 1;
xo.method = 2;
xo.ls_maxatmov = 0.15;

m = getpdb("d1_t23_snap1.pdb");
natm = m.natoms;
allocate xyz[3*natm]; allocate grad[3*natm];
readparm(m,"d1_t23_complex_mbondi2.prmtop");

mm_options("cut=999.0, ntpr=1, nsnb=99999, ipb=2, inp=1, radiopt=0,
cavity_surften=0.005420, cavity_offset=0.920000");
mme_init(m, NULL, "::Z", xyz, NULL);
setxyz_from_mol(m,NULL,xyz);

mme(xyz,grad,0);
xmin(mme,natm,xyz,grad,energy,grms,xo);
putxyz("d1_t23_snap1.crd",natm,xyz);

=======================================================================================

And then, below is the .out file I got after successfully compiling the nab
file:

=======================================================================================
Reading parm file (d1_t23_complex_mbondi2.prmtop)
title:
default_name
        mm_options: cut=999.0
        mm_options: ntpr=1
        mm_options: nsnb=99999
        mm_options: ipb=2
        mm_options: inp=1
        mm_options: radiopt=0
        mm_options: cavity_surften=0.005420
        mm_options: cavity_offset=0.920000
      iter Total bad vdW elect nonpolar EPB
 frms
PB Warning in pb_read(): the cavity energy should be correlated to the
SASA when inp=1, so use_sav is reset to 0.
PB Warning in pb_read(): sprob=0.557 is optimized for inp=2 and should not
be used with inp=1. It has been reset to 1.4.
ff: 0 -1052.93 3041.05 -335.43 -2712.21 36.00 -1082.35
 1.76e+01
Starting L-BFGS minimisation
PB Warning in pb_read(): the cavity energy should be correlated to the
SASA when inp=1, so use_sav is reset to 0.
PB Warning in pb_read(): sprob=0.557 is optimized for inp=2 and should not
be used with inp=1. It has been reset to 1.4.
 PB Bomb in pb_init(): Allocation aborted 151 151
151
         151 151 151 151 151 151
         151 151 151 151 151 151
         151 151 151 151 151 151
         151 151 151 151 0 151
         151 151 0 151

=======================================================================================

The two warning messages, I don't think they are big problems. Seems like
they can still continue. Actually I can specify parameters for sprob in
mm_options() but use_sav is not recognized in mm_options(). Not sure how
can I eliminate this warning.

But anyway, the biggest problem is "PB Bomb in pb_init(): Allocation
aborted". I browsed around it might related to memory issue(?) but I'm not
so sure about that. I also tried to change spacing parameters for grid but
seems that doesn't help either. Can anybody help me with that? Any
suggestion or enlightenment is appreciated. I just want to know the reason
behind it...

Thank you so much!

Best,
Guqin

-- 
Guqin SHI
The Ohio State University
College of Pharmacy
500 W. 12th Ave.
Columbus, OH, 43210
(614)688-3531
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri May 09 2014 - 12:00:02 PDT
Custom Search