[AMBER] MMPBSA.py problem

From: Asmita Gupta <asmita4des.gmail.com>
Date: Tue, 25 Jun 2013 12:37:05 +0530

Dear users,


I want to do some stability calculations on an RNA aptamer structure using
MMPB(GB)/SA calculations. there are no ligands,so it's just a single
species in the system.

I submiited the following script for free energy and nmode analysis:

Input file for running GB/nmode
&general
  debug_printlevel=3,interval=10,netcdf=1
# entropy=1,
/
&gb
  igb=5, saltcon=0.100
/
&decomp
  idecomp=1
/
&nmode
  maxcyc=1000,nminterval=10,nmode_igb=1,nmode_istrng=0.1
/

and subsequently,

MMPBSA.py -O -i mmpbsa.in -o results_gb_decomp_nmode.dat -sp pc.prmtop -cp
pc_stripped.prmtop -y ../../production/run_10.crd -do decomp.out -eo
energy.out -deo decomp_energy.out

so i have basically given solvated topologies and because it's a single
molecule, a dry topology file for -cp option....is this method run correct??

a total of 1000 frames will be analyzed in this process...

my problem is after submitting it, the process just keeps on going, and now
i am getting this message:

Line minimizer aborted: step at upper bound 0.06053411
Line minimizer aborted: step at upper bound 0.10156325
Line minimizer aborted: step at upper bound 0.080443891
Line minimizer aborted: step at upper bound 0.069342559
Line minimizer aborted: step at upper bound 0.083281256
...
..

.....is it a normal behaviour?? To me it seems that the minimization
process in nmode is going wrong somewhere...Pl suggest..

Thank you

Asmi
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jun 25 2013 - 00:30:02 PDT
Custom Search