AMBER: TI calculation using sander

From: Latifa Douali <latifa.douali.pnl.gov>
Date: Tue, 07 Jun 2005 17:10:54 -0700

Hi all,


I am trying to redo the TI example given in the tutorial of amber8. I
run the job as described in the tutorial. Sudenly, the calculation
crashed with this error message:

> unreasonably small cut for non-periodic run: 8.00
>
> *** input error(s)


My output file looks like:

Here is the input file:
                                                                                                                      

 neutralize toluene
&cntrl
  ntr=0,
  nstlim = 100000, nscm = 2000, ntave = 10000,
  ntx =5, irest=1, ntb=0, ntpr=200,
  dt=0.001, nrespa=2,
  ntt=1, temp0=300., tautp=2.0,
  ntc=2, ntf=2, tol=0.000001,
  ntwr = 10000, ntwx=0,
  icfe=1, clambda=0.5,
/
                                                                                                                      

--------------------------------------------------------------------------------
   1. RESOURCE USE:
--------------------------------------------------------------------------------
                                                                                                                      

| Flags:
| NONPERIODIC ntb=0 and igb=0: Setting up nonperiodic simulation
| *** cutoff > system size, list only builds once
|Largest sphere to fit in unit cell has radius = 31.913
| New format PARM file being parsed.
| Version = 1.000 Date = 06/07/05 Time = 10:37:12
 NATOM = 15 NTYPES = 4 NBONH = 8 MBONA = 7
 NTHETH = 16 MTHETA = 8 NPHIH = 27 MPHIA = 9
 NHPARM = 0 NPARM = 0 NNB = 68 NRES = 1
 NBONA = 7 NTHETA = 8 NPHIA = 9 NUMBND = 4
 NUMANG = 5 NPTRA = 3 NATYP = 4 NPHB = 0
 IFBOX = 0 NMXRS = 15 IFCAP = 0 NEXTRA = 0
 NCOPY = 0
                                                                                                                      

                                                                                                                      

| Memory Use Allocated
| Real 2216
| Hollerith 93
| Integer 388287
| Max Pairs 105
| Max Rstack 50195
| Max Istack 75
| Total 1927 kbytes
                                                                                                                      

Running a free energy calculation with lambda = 0.500
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals
                                                                                                                      

--------------------------------------------------------------------------------
   2. CONTROL DATA FOR THE RUN
--------------------------------------------------------------------------------
                                                                                                                      

PHE
                                                                                                                      

General flags:
     imin = 0, nmropt = 0
                                                                                                                      

Nature and format of input:
     ntx = 5, irest = 1, ntrx = 1
                                                                                                                      

Nature and format of output:
     ntxo = 1, ntpr = 200, ntrx = 1, ntwr
= 10000
     iwrap = 0, ntwx = 0, ntwv = 0, ntwe
= 0
     ioutfm = 0, ntwprt = 0, idecomp = 0,
rbornstat= 0
                                                                                                                      

Potential function:
     ntf = 2, ntb = 0, igb = 0, nsnb
= 25
     ipol = 0, gbsa = 0, iesp = 0
     dielc = 1.00000, cut = 8.00000, intdiel = 1.00000
     scnb = 2.00000, scee = 1.20000
                                                                                                                      

Frozen or restrained atoms:
     ibelly = 0, ntr = 0
                                                                                                                      

Molecular dynamics:
     nstlim = 100000, nscm = 2000, nrespa = 2
     t = 0.00000, dt = 0.00100, vlimit = 20.00000
                                                                                                                      

Berendsen (weak-coupling) temperature regulation:
     temp0 = 300.00000, tempi = 0.00000, tautp = 2.00000
                                                                                                                      

SHAKE:
     ntc = 2, jfastw = 0
     tol = 0.00000
                                                                                                                      

Free energy options:
     klambda = 1
     clambda = 0.50000
                                                                                                                      

 unreasonably small cut for non-periodic run: 8.00
                                                                                                                      

 *** input error(s)
--------------------------------------------------------


Thank you for your help

Latifa






-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jun 08 2005 - 01:53:01 PDT
Custom Search