AMBER: (no subject)

From: Syed Tarique Moin <tarisyed.yahoo.com>
Date: Thu, 21 Feb 2008 22:10:36 -0800 (PST)

Hello,
   
  I am trying to minimize a complex, but i was stuck by the linmin error. I am attaching output file of minimization.
   
  Thanks in advance
   
  Regards

       
---------------------------------
Be a better friend, newshound, and know-it-all with Yahoo! Mobile. Try it now.


          -------------------------------------------------------
          Amber 9 SANDER 2006
          -------------------------------------------------------

| Run on 02/22/2008 at 00:53:04
  [-O]verwriting output

File Assignments:
| MDIN: min500res.in
| MDOUT: min500res.out
|INPCRD: a_wat.inpcrd
| PARM: a_wat.prmtop
|RESTRT: min500res.restrt
| REFC: a_wat.inpcrd
| MDVEL: mdvel
| MDEN: mden
| MDCRD: mdcrd
|MDINFO: mdinfo
|INPDIP: inpdip
|RSTDIP: rstdip

 
 Here is the input file:
 
Minimization Amber9
 &cntrl
 imin=1, maxcyc=500,
 ntpr=5, ntr=1,
 &end
Restraints
 25.0
RES 1 633
END
END

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

| Flags: MPI
 getting new box info from bottom of inpcrd
| INFO: Old style inpcrd file read

| peek_ewald_inpcrd: Box info found
|Largest sphere to fit in unit cell has radius = 41.112
| New format PARM file being parsed.
| Version = 1.000 Date = 02/21/08 Time = 16:43:42
 NATOM = 59490 NTYPES = 18 NBONH = 54486 MBONA = 5133
 NTHETH = 10821 MTHETA = 6970 NPHIH = 20811 MPHIA = 15935
 NHPARM = 0 NPARM = 0 NNB = 120143 NRES = 17214
 NBONA = 5133 NTHETA = 6970 NPHIA = 15935 NUMBND = 51
 NUMANG = 109 NPTRA = 50 NATYP = 34 NPHB = 1
 IFBOX = 1 NMXRS = 24 IFCAP = 0 NEXTRA = 0
 NCOPY = 0


| Memory Use Allocated
| Real 4071891
| Hollerith 374156
| Integer 2024999
| Max Pairs 2478750
| nblistReal 713880
| nblist Int 2496746
| Total 66196 kbytes
| Duplicated 0 dihedrals
| Duplicated 0 dihedrals

     BOX TYPE: RECTILINEAR

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

                                                                                

General flags:
     imin = 1, nmropt = 0

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

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

Potential function:
     ntf = 1, ntb = 1, 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 = 1

Energy minimization:
     maxcyc = 500, ncyc = 10, ntmin = 1
     dx0 = 0.01000, drms = 0.00010

Ewald parameters:
     verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1
     vdwmeth = 1, eedmeth = 1, netfrc = 0
     Box X = 82.224 Box Y = 96.666 Box Z = 91.967
     Alpha = 90.000 Beta = 90.000 Gamma = 90.000
     NFFT1 = 90 NFFT2 = 96 NFFT3 = 96
     Cutoff= 8.000 Tol =0.100E-04
     Ewald Coefficient = 0.34864
     Interpolation order = 4

    LOADING THE CONSTRAINED ATOMS AS GROUPS


   5. REFERENCE ATOM COORDINATES

      
    ----- READING GROUP 1; TITLE:
 Restraints

     GROUP 1 HAS HARMONIC CONSTRAINTS 25.00000
 GRP 1 RES 1 TO 633
      Number of atoms in this group = 9747
    ----- END OF GROUP READ -----

--------------------------------------------------------------------------------
   3. ATOMIC COORDINATES AND VELOCITIES
--------------------------------------------------------------------------------

                                                                                
 begin time read from input coords = 0.000 ps

 Number of triangulated 3-point waters found: 16581
| Atom division among processors:
| 0 7436 14874 22308 29745 37182 44616 52053
| 59490

     Sum of charges from parm topology file = -0.32360874
     Assuming uniform neutralizing plasma
| Running AMBER/MPI version on 8 nodes


--------------------------------------------------------------------------------
   4. RESULTS
--------------------------------------------------------------------------------

 ---------------------------------------------------
 APPROXIMATING switch and d/dx switch using CUBIC SPLINE INTERPOLATION
 using 5000.0 points per unit in tabled values
 TESTING RELATIVE ERROR over r ranging from 0.0 to cutoff
| CHECK switch(x): max rel err = 0.2738E-14 at 2.422500
| CHECK d/dx switch(x): max rel err = 0.8314E-11 at 2.736960
 ---------------------------------------------------
| Local SIZE OF NONBOND LIST = 1405134
| TOTAL SIZE OF NONBOND LIST = 10336257


   NSTEP ENERGY RMS GMAX NAME NUMBER
      1 5.7364E+17 1.5567E+17 3.1549E+19 CA 3075

 BOND = 243646.1205 ANGLE = 4100.3728 DIHED = 5846.6820
 VDWAALS = ************* EEL = -204516.9558 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 31395.6529 RESTRAINT = 0.0000


   NSTEP ENERGY RMS GMAX NAME NUMBER
      5 3.7073E+16 8.0072E+15 1.6163E+18 CA 3075

 BOND = 243654.1386 ANGLE = 4101.4930 DIHED = 5846.8190
 VDWAALS = ************* EEL = -204516.5915 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 31125.5931 RESTRAINT = 0.0180
 EAMBER = *************


   NSTEP ENERGY RMS GMAX NAME NUMBER
     10 1.4566E+14 1.9825E+13 4.0017E+15 CA 3075

 BOND = 243680.8699 ANGLE = 4105.7876 DIHED = 5847.3354
 VDWAALS = ************* EEL = -204515.5826 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30736.0821 RESTRAINT = 0.2704
 EAMBER = *************

     .... RESTARTED DUE TO LINMIN FAILURE ...


   NSTEP ENERGY RMS GMAX NAME NUMBER
     15 3.5900E+13 4.3453E+12 8.7713E+14 CA 3075

 BOND = 243691.0310 ANGLE = 4107.5708 DIHED = 5847.5505
 VDWAALS = ************* EEL = -204515.2549 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30663.0602 RESTRAINT = 0.4214
 EAMBER = *************

     .... RESTARTED DUE TO LINMIN FAILURE ...


   NSTEP ENERGY RMS GMAX NAME NUMBER
     20 3.5782E+13 4.3299E+12 8.7401E+14 CA 3075

 BOND = 243691.0571 ANGLE = 4107.5755 DIHED = 5847.5511
 VDWAALS = ************* EEL = -204515.2541 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30662.8983 RESTRAINT = 0.4219
 EAMBER = *************

     .... RESTARTED DUE TO LINMIN FAILURE ...


   NSTEP ENERGY RMS GMAX NAME NUMBER
     25 3.5664E+13 4.3144E+12 8.7089E+14 CA 3075

 BOND = 243691.0832 ANGLE = 4107.5801 DIHED = 5847.5517
 VDWAALS = ************* EEL = -204515.2533 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30662.7366 RESTRAINT = 0.4223
 EAMBER = *************

     .... RESTARTED DUE TO LINMIN FAILURE ...


   NSTEP ENERGY RMS GMAX NAME NUMBER
     30 3.5547E+13 4.2991E+12 8.6779E+14 CA 3075

 BOND = 243691.1093 ANGLE = 4107.5848 DIHED = 5847.5522
 VDWAALS = ************* EEL = -204515.2524 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30662.5750 RESTRAINT = 0.4227
 EAMBER = *************

     .... RESTARTED DUE TO LINMIN FAILURE ...

     ***** REPEATED LINMIN FAILURE *****

     ***** SEE http://amber.scripps.edu/Questions/linmin.html FOR MORE INFO *****


                    FINAL RESULTS



   NSTEP ENERGY RMS GMAX NAME NUMBER
     31 3.5509E+13 4.2941E+12 8.6680E+14 CA 3075

 BOND = 243691.1177 ANGLE = 4107.5863 DIHED = 5847.5524
 VDWAALS = ************* EEL = -204515.2522 HBOND = 0.0000
 1-4 VDW = ************* 1-4 EEL = 30662.5229 RESTRAINT = 0.4228
 EAMBER = *************

--------------------------------------------------------------------------------
   5. TIMINGS
--------------------------------------------------------------------------------

| Read coords time 0.09 ( 0.19% of Total)
| Fast Water setup 0.01 ( 0.01% of Total)
| Build the list 0.74 (53.93% of List )
| Other 0.63 (46.07% of List )
| List time 1.36 ( 3.89% of Nonbo)
| Short_ene time 14.01 (90.69% of Direc)
| Other 1.44 ( 9.31% of Direc)
| Direct Ewald time 15.45 (45.85% of Ewald)
| Adjust Ewald time 0.07 ( 0.21% of Ewald)
| Fill Bspline coeffs 0.70 ( 4.34% of Recip)
| Fill charge grid 0.73 ( 4.56% of Recip)
| Scalar sum 0.43 ( 2.67% of Recip)
| Grad sum 0.44 ( 2.71% of Recip)
| FFT back comm time 9.35 (74.85% of FFT t)
| Other 3.14 (25.15% of FFT t)
| FFT time 12.50 (77.90% of Recip)
| Other 1.25 ( 7.82% of Recip)
| Recip Ewald time 16.04 (47.59% of Ewald)
| Force Adjust 1.60 ( 4.76% of Ewald)
| Virial junk 0.53 ( 1.58% of Ewald)
| Ewald time 33.70 (96.11% of Nonbo)
| Nonbond force 35.07 (80.14% of Force)
| Bond/Angle/Dihedral 0.22 ( 0.51% of Force)
| FRC Collect time 4.22 ( 9.64% of Force)
| Other 4.25 ( 9.71% of Force)
| Force time 43.76 (100.0% of Runmd)
| Runmd Time 43.76 (92.34% of Total)
| Other 3.54 ( 7.46% of Total)
| Total time 47.39 (100.0% of ALL )

| Highest rstack allocated: 588891
| Highest istack allocated: 12393
| Job began at 00:53:04.366 on 02/22/2008
| Setup done at 00:53:07.437 on 02/22/2008
| Run done at 00:53:51.754 on 02/22/2008
| wallclock() was called 1374 times

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Feb 24 2008 - 06:07:35 PST
Custom Search