Re: [AMBER] slope

From: Jio M <jiomm.yahoo.com>
Date: Mon, 16 Nov 2009 20:07:07 -0800 (PST)

Dear Sir,





>the temperature plot should look like what you requested in your input.

>check the mdout file- does sander recognize your request for temperature

>change?



here is the part of output annealing

#############################

      Begin reading energy term weight changes/NMR restraints

 WEIGHT CHANGES:

 TEMP0         0
140000  300.000000 2000.000000     
0      0

 TEMP0    140001 200000 2000.000000
2000.000000     
0      0

 TEMP0    200001 340000 2000.000000
1000.000000     
0      0

 TEMP0    340001 480000 1000.000000
2000.000000     
0      0

 TEMP0    480001 620000 2000.000000
1000.000000     
0      0

 TEMP0    620001 760000 1000.000000
2000.000000     
0      0

 TEMP0    760001 900000 2000.000000
1000.000000     
0      0

 TEMP0    9000011050000 1000.000000 
300.000000     
0      0



 RESTRAINTS:

                         
** No restraint defined **



                 
Done reading weight changes/NMR restraints



I think (please correct me if I am wrong) above indicates temp change.



 Also I plotted graph of Temp vs Time. it seems that the temp
(y)  approx. corresponds to y = mxi + c equation [ where m should
be one  and i is correction factor ,to consider different division
factor of x and y axis i.e. i = delta y/delta x]



>note that a temp0 of 2000 will almost certainly give you problems unless you

>add additional restraints for things like chirality.

>this is a very high energy state and you can easily get cis-trans

>isomerizations and chirality inversions.



Actually I have a macromolecule (> 1300 atoms) that donot have
chiral center or cis trans possibility. But in future I will definately
consider this while doing high temp annealing, this is nice point to be
noticed that I didnt know



thanks and regards;

JIomm


--- On Mon, 11/16/09, Carlos Simmerling <carlos.simmerling.gmail.com> wrote:

From: Carlos Simmerling <carlos.simmerling.gmail.com>
Subject: Re: [AMBER] slope
To: "AMBER Mailing List" <amber.ambermd.org>
Date: Monday, November 16, 2009, 11:11 AM

the temperature plot should look like what you requested in your input.
check the mdout file- does sander recognize your request for temperature
change?

note that a temp0 of 2000 will almost certainly give you problems unless you
add additional restraints for things like chirality.
this is a very high energy state and you can easily get cis-trans
isomerizations and chirality inversions.

On Sun, Nov 15, 2009 at 11:42 PM, Jio M <jiomm.yahoo.com> wrote:

> Dear Sir,
>
>
>
> I am using following input file
>
> &cntrl
>
>
>   nstlim=1050000,
>
>   ntt=3, gamma_ln=5,
>
>
>   scee=1.2,
>
>
>   ntpr=500, pencut=0.1,
>
>
>   nmropt=1,
>
>
>   vlimit=10,
>
>
>   ntb=0,dt = 0.0005,cut=999,ntwx=500
>
> /
>
>  &wt type='TEMP0', istep1=0,istep2=140000,
>
>    value1=300.0, value2=2000.0
>
>  /
>
>  &wt type='TEMP0', istep1=140001,istep2=200000,
>
>    value1=2000.0, value2=2000.0
>
>  /
>
>  &wt type='TEMP0', istep1=200001,istep2=340000,
>
>    value1=2000.0, value2=1000.0
>
>  /
>
>  &wt type='TEMP0', istep1=340001,istep2=480000,
>
>    value1=1000.0, value2=2000.0
>
>  /
>
>  &wt type='TEMP0', istep1=480001,istep2=620000,
>
>    value1=2000.0, value2=1000.0
>
>  /
>
>  &wt type='TEMP0', istep1=620001,istep2=760000,
>
>    value1=1000.0, value2=2000.0
>
> /
>
>  &wt type='TEMP0', istep1=760001,istep2=900000,
>
>    value1=2000.0, value2=1000.0
>
> /
>
>  &wt type='TEMP0', istep1=900001,istep2=1050000,
>
>    value1=1000.0, value2=300.0
>
>  /
>
>  &wt type='END'
>
>  /
>
>
>
>
>
> --- On Sun, 11/15/09, Carlos Simmerling <carlos.simmerling.gmail.com>
> wrote:
>
> From: Carlos Simmerling <carlos.simmerling.gmail.com>
> Subject: Re: [AMBER] slope {no reply}
> To: "AMBER Mailing List" <amber.ambermd.org>
> Date: Sunday, November 15, 2009, 7:07 PM
>
> for (1) it all depends on how you set the thermostat temp0, which
> thermostat and what coupling constant. no way to predict your data
> without that.
>
> On 11/15/09, Jio M <jiomm.yahoo.com> wrote:
> >
> >
> >
> >
> > Dear Amber users,
> >
> >
> >
> > 1) I am performing annealing and want to plot temp vs steps graph ..is
> > it neccessary that  temp vs  steps plot  should give
> > linear rise or fall of graph of slope 45 degree or it should be linear
> > what so ever slope may be ?
> >
> >
> >
> > 2) same query as before.. can we solubilize the system with single
> molecule
> > instead of box solvent
> >
> >
> >
> > please suggest
> >
> >
> >
> > thanks and regards;
> >
> > Jiomm
> >
> >
> >
> >
> >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
> >
> >
> >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
>
>
> --
> ===================================================================
> Carlos L. Simmerling, Ph.D.
> Professor, Department of Chemistry
> CMM Bldg, Room G80           Phone: (631) 632-1336   Fax: 632-1555
> Stony Brook University           E-mail: carlos.simmerling.gmail.com
> Stony Brook, NY 11794-5115 Web: http://www.simmerlinglab.org
> ===================================================================
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Nov 16 2009 - 20:30:04 PST
Custom Search