Re: [AMBER] trajectory of minimization using generalized Born

From: Jason Swails <jason.swails.gmail.com>
Date: Wed, 18 Mar 2015 17:03:42 -0400

On Wed, Mar 18, 2015 at 12:41 PM, Edwin Helbert Aponte Angarita <
helbert2a.gmail.com> wrote:

> Dear All,
>
> I am able to write the trajectory of a minimization of a peptide salvated
> in TIP3PBOX. I have used this .in file
>
> &cntrl
> imin=1,
> maxcyc=20000,
> ncyc=4000,
> cut=10.0,
> ntb=1,
> ntr=0,
> ntpr=100,
> ntwx=100,
> ioutfm=1,
> /
>
> Then I calculate the rmsd to the first frame using ptraj and whith:
> trajin ./min.mdcrd
> rms first out rmsd_to_first.dat .N,CA,C
>
> and I get a very nice plot showing how the peptide converges to a local
> minimum and does not change noticeable any more. However, I cannot get
> sander nor sander.MPI to write the trajectory file if I use generalized
> Born implicit solvent. In that case I use:
>
> &cntrl
> imin=1,
> maxcyc=2000,
> ncyc=1000,
> cut=999.,
> rgbmax=999.,
> igb=1,
> ntb=0,
> ntpr=50,
> ntwx=50,
> ioutfm=1,
> /
>
> and I am setting approppriatelly -i -o -p -c -r -x.
>
> Do you know how can I get the trajectory when using generalized Born
> implicit solvent? I only found a very old similar post (
> http://archive.ambermd.org/200502/0076.html) but is not really helpful.
>

​I have no idea why it's not working for you. It works just fine for me
with Amber 12 and AmberTools 13 (fully updated):

swails.batman /usr/local/amber12/test/gb7_trx $ ../../update_amber --version

Version is reported as <version>.<patches applied>


 AmberTools version 13.24
    Amber version 12.21


​if you're sure it's not working for you, then let us know what operating
system and compilers you are using to see if we can try to reproduce it
somewhere...

HTH,
Jason

-- 
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 18 2015 - 14:30:02 PDT
Custom Search