[AMBER] Simulated annealing

From: waleed zalloum <waleed_zalloum.yahoo.com>
Date: Thu, 15 Oct 2009 10:59:29 +0000 (GMT)

Dear Amber users,
 
I am calculating a 3D structure of a oligopeptide conjugated to single stranded DNA. I run NMR restrained simulated annealing. I specified certain heat range for certain number of steps ( as in the input file bellow). Upon running,  the heat does not stick to the  specified  number of steps and there is a huge fluctuation in the temperature (eg. from 600 going to 800 and then go back to 650). What is the problem? is there any problem in the input file?
 
I specified the LISTOUT option, but I don't get any output list coming after the end of simulation? what is the problem in that?
 
ps. Bellow is the annealing input file.
 
Thank you

  
 
Waleed
                                                                             
 &cntrl                                                                       
    nstlim=190000, pencut=-0.001, nmropt=1,ipnlty=1,scee=1.2,
    ntpr=200, ntt=3, ntwx=200,imin=0,gamma_ln=1.0,irest=0,ntx=1,
    cut=9.0, ntb=0, vlimit=10, dt=0.001, igb=1, ntc=2, ntf=2,temp0=500
 /                                                                       
 &ewald
    eedmeth=5,
 /
#                                                                              
                                                                              
 &wt type='TEMP0', istep1=0,istep2=2000,value1=500.0,
            value2=500.0,    /
 &wt type='TEMP0', istep1=2001, istep2=5000, value1=500.0,
            value2=1000.0,     /
 &wt type='TEMP0', istep1=5001, istep2=10000, value1=1000.0,
            value2=1000.0,     /
 &wt type='TEMP0', istep1=10001, istep2=20000.0, value1=1000.0,
            value2=950.0,     /
 &wt type='TEMP0', istep1=20001, istep2=30000.0, value1=950.0,
            value2=900.0,     /
 &wt type='TEMP0', istep1=30001, istep2=40000.0, value1=900.0,
            value2=850.0,     /
 &wt type='TEMP0', istep1=40001, istep2=50000, value1=850.0,
            value2=800.0,     /
 &wt type='TEMP0', istep1=50001, istep2=60000, value1=800.0,
            value2=750.0,     /
 &wt type='TEMP0', istep1=60001, istep2=70000, value1=750.0,
            value2=700.0,     /
 &wt type='TEMP0', istep1=70001, istep2=80000, value1=700.0,
            value2=650.0,     /
 &wt type='TEMP0', istep1=80001, istep2=90000, value1=650.0,
            value2=600.0,     /
 &wt type='TEMP0', istep1=90001, istep2=100000, value1=600.0,
            value2=550.0,     /
 &wt type='TEMP0', istep1=100001, istep2=110000, value1=550.0,
            value2=500.0,     /
 &wt type='TEMP0', istep1=110001, istep2=120000, value1=500.0,
            value2=450.0,     /
 &wt type='TEMP0', istep1=120001, istep2=130000, value1=450.0,
            value2=400.0,     /
 &wt type='TEMP0', istep1=130001, istep2=140000, value1=400.0,
            value2=350.0,     /
 &wt type='TEMP0', istep1=140001, istep2=150000, value1=350.0,
            value2=300.0,     /
 &wt type='TEMP0', istep1=150001, istep2=160000, value1=300.0,
            value2=250.0,     /
 &wt type='TEMP0', istep1=160001, istep2=170000, value1=250.0,
            value2=200.0,     /
&wt type='TEMP0', istep1=170001, istep2=180000, value1=200.0,
            value2=150.0,     /
&wt type='TEMP0', istep1=180001, istep2=190000, value1=150.0,
            value2=100.0,     /
 

 &wt type='TAUTP', istep1=0,istep2=10000,value1=0.4,
            value2=0.4,     /
 &wt type='TAUTP', istep1=10001,istep2=170000,value1=4.0,
            value2=4.0,     /
 &wt type='TAUTP', istep1=170001,istep2=180000,value1=1.0,
            value2=1.0,     /
 &wt type='TAUTP', istep1=180001,istep2=190000,value1=0.1,
            value2=0.1,     /                   
                           
 &wt type='REST', istep1=0,istep2=5000,value1=0.0,                            
            value2=0.0,  /                                               
 &wt type='REST', istep1=5001,istep2=10000,value1=0.0,                       
            value2=0.01,   /
 &wt type='REST', istep1=10001,istep2=20000,value1=0.01,
            value2=0.1,  /
 &wt type='REST', istep1=20001,istep2=30001,value1=0.1,
            value2=0.15,  /
 &wt type='REST', istep1=30001,istep2=40000,value1=0.15,
            value2=0.2,  /
 &wt type='REST', istep1=40001,istep2=50000,value1=0.2,
            value2=0.25,  /
 &wt type='REST', istep1=50001,istep2=60000,value1=0.25,
            value2=0.3,  /
 &wt type='REST', istep1=60001,istep2=100000,value1=0.3,
            value2=0.6,  /
 &wt type='REST', istep1=100001,istep2=150000,value1=0.6,
            value2=0.9,  /
 &wt type='REST', istep1=150001,istep2=190000,value1=0.9,
            value2=1.0,  /                                             
                                                                              
 &wt type='END'  /
LISTOUT=POUT                                                                  
DISANG=RST       

 

 
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Oct 15 2009 - 04:30:03 PDT
Custom Search