[AMBER] Questions about nmode Command

From: 张明焜 <zhangmk69.gmail.com>
Date: Sun, 4 May 2014 15:59:01 +0800

Questions about nmode Command

Hellow specialists:

I want to caculate the normal mode frequences of a protein with nmode
command.first ,I do Energy minimization in sander and got the local energy
minimum conformation Coordinate “min.rst” .then I Edit nmode script
“nmd.in”like
below:

normal mode analysis

 &data

 ntrun=1,

 maxcyc=0,

 nsave=250, nprint=250,

 cut=10.0,

 scnb=2.0, scee=1.2,

 nvect=30, ilevel=0,

 t=310.0,

and Enter the command like :

$AMBERHOME/bin/nmode -O -i nmd.in -o nmd.out -c min.rst -p Glu.prmtop -r
Glu.rst -v Glu.vecs

Then I get an error:

***************** ***************** *****************

   step = 0

     F = -0.526752E+02 GRDMAX = 0.525215E+01 GNORM = 0.155207E+01

       E-NONB E-ELE E-HBOND E-BOND

     0.11450E+02 -0.84061E+02 0.00000E+00 0.38022E+00

       E-ANGLE E-DIHED E-NB14 E-EEL14

     0.24024E+01 0.40266E+01 0.15051E+01 0.11621E+02

       E-POL E-3BOD

     0.00000E+00 0.00000E+00

 Root-mean-square gradient of input coords is 1.55207279918410

 This is greater than the requested maximum: 1.000000000000000E-005



Then check Each energy component of sander result and nmode result,I found
the E-EEL14 item is different, which is 38.2626 in sander. What's this
about? So this is my first queston.

     Later ,I turn up the drms to 5.0 to to avoid this error,then I
got the normal
mode eigenvalues in “Glu.out”,like

   freq. E Cv S

          cm**-1 kcal/mol cal/mol-kelvin cal/mol-kelvin

--------------------------------------------------------------------------------

 Total 75.455 37.957 102.989

 translational 0.923 2.979 41.114

 rotational 0.923 2.979 29.847

 vibrational 73.608 32.000 32.027

     1 -0.014

     2 0.000

     3 0.000

     4 0.000

     5 15.956

     6 17.041

     7 51.622 0.619 1.976 4.827

     8 90.323 0.625 1.957 3.725

     9 124.291 0.633 1.932 3.104

    10 131.321 0.635 1.925 2.998

    11 148.172 0.640 1.909 2.767

    12 194.258 0.657 1.856 2.256

    13 260.820 0.689 1.760 1.722

    14 267.339 0.693 1.749 1.679

    15 286.182 0.704 1.718 1.561

    16 379.039 0.767 1.544 1.101

    17 399.342 0.783 1.503 1.021

18 428.541 0.807 1.442 0.917

But I can’t get the eigenvectors in Glu.vecs, This file was empty! Where
can I got normal mode eigenvectors? So this is my second question.


thank you for your early reply!


-- 
---------------------------------------------------------------------
ZhangMK
+8623-63063582
Chongqing Institute of Green and Intelligent Technology of Chinese Academy
Sciences.
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun May 04 2014 - 01:30:03 PDT
Custom Search