Re: AMBER: simulated anealing ...vlimit exceed

From: Carlos Simmerling <carlos.ilion.bio.sunysb.edu>
Date: Thu, 6 May 2004 07:47:48 -0400

you have VERY high initial vdw energy-
was the structure minimized first? can you run
any normal MD? it's good to try that before a
complicated input file.
carlos

----- Original Message -----
From: <anshul.imtech.res.in>
To: "amber mail reflector" <amber.scripps.edu>
Sent: Thursday, May 06, 2004 5:50 PM
Subject: AMBER: simulated anealing ...vlimit exceed


> hi all,
> I am trying to run simulated anealing on a protein but everytime i am
> getting the error of vlimit exceed. I am running it in vacuume conditions
> and according to the archive of mails i also reduced teh force constant
> for restraints to 2 ..but still the same thing. can anyone suggest
> something.
>
>
> input file:
>
> &ewald
> eedmeth=5,
> &end
> &cntrl
> imin=0, cut=30, nstlim=90000, ntt=1, ntx=1, irest=0,
> ntpr=500, ntwx=500, ntwe=500, tempi=10,
> ntb=0, ntr=1, dielc=4, vlimit=30,
> &end
> Input for restrained residues (A1-A166)
> 2
> RES 1 167
> END
> Input for restrained residues (A182-A276)
> 2
> RES 182 276
> END
> END
> &wt type='TEMP0', istep1=0,istep2=50000,value1=10.0,
> value2=500.0, &end
> &wt type='TEMP0', istep1=50001,istep2=60000,value1=500.0,
> value2=500.0, &end
> &wt type='TEMP0', istep1=60001,istep2=70000,value1=500.0,
> value2=300.0, &end
> &wt type='TEMP0', istep1=70001,istep2=90000,value1=300.0,
> value2=300.0, &end
> &wt type='TAUTP', istep1=0,istep2=60000,value1=0.5,
> value2=0.5 &end
> &wt type='TAUTP', istep1=60001,istep2=65000,value1=2.0,
> value2=2.0 &end
> &wt type='TAUTP', istep1=65001,istep2=70000,value1=0.2,
> value2=0.2 &end
> &wt type='TAUTP', istep1=70001,istep2=90000,value1=0.5,
> value2=0.5 &end
> &wt type='END' &end
>
> output::
>
> RESULTS
> --------------------------------------------------------------------------
------
>
> ---------------------------------------------------
>
> eedmeth=4: Setting switch to one everywhere
>
> ---------------------------------------------------
> | Local SIZE OF NONBOND LIST = 4988756
> | TOTAL SIZE OF NONBOND LIST = 4988756
>
> NSTEP = 0 TIME(PS) = 0.000 TEMP(K) = 9.88 PRESS = 0.0
> Etot = ************ EKtot = 129.0390 EPtot = ************
> BOND = 772.3704 ANGLE = 1001.8422 DIHED = 2404.3438
> 1-4 NB = 1253.7810 1-4 EEL = 2187.4793 VDWAALS = ************
> EELEC = -4664.7256 EHBOND = 0.0000 RESTRAINT = 0.0000
>
> --------------------------------------------------------------------------
----
>
>  vlimit exceeded for step  0; vmax =   2.16726961E+12
>  vlimit exceeded for step  1; vmax =   1.25533206E+12
>  vlimit exceeded for step  2; vmax =   53538856.1
>  vlimit exceeded for step  3; vmax =   445647.94
>  vlimit exceeded for step  4; vmax =   1058268.05
>  vlimit exceeded for step  5; vmax =   2.00528168E+11
>  vlimit exceeded for step  6; vmax =   767387.095
>  vlimit exceeded for step  7; vmax =   322935852.
>  vlimit exceeded for step  8; vmax =   17379772.
>  vlimit exceeded for step  9; vmax =   4566640.04
>
>
>            With regards
> Anshul
>
>
>
>
>
>
>
>
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
>
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu May 06 2004 - 13:53:00 PDT
Custom Search