[gmx-users] 1/viscosity in g_energy

[Alif M Latif]

Dear GROMACS users and developers,

I'm having trouble getting values for my 1/viscosity calculation which obtained 
from g_energy:

g_energy -f md1_vis.edr -o md1_vis_1perv.xvg


the output:

Statistics over 1000001 steps [ 0.0000 through 2000.0000 ps ], 1 data sets
All statistics are over 100001 points

Energy                      Average   Err.Est.       RMSD  Tot-Drift
-------------------------------------------------------------------------------
1/Viscosity                       0          0          0          0  (m s/kg)

gcq#310: "Shoot them in the back now" (The Ramones)


 Have I done something wrong with my simulation? If that's the case, then this 
is my mdp parameter:

cpp             = /lib/cpp
include          = -I../top
integrator         = md
dt             = 0.002
nsteps           = 10000
nstxout          = 100
nstvout          = 100
nstlog           = 1000
nstenergy         = 100
nstxtcout         = 100
xtc_grps         = System    
energygrps         = System  
nstlist          = 10 
ns_type          = grid

rlist             = 1.0

; PME chosen as the best option to calculate viscosity ;
coulombtype         = PME
rcoulomb         = 1.0

rvdw             = 1.4

pbc                      = xyz

tcoupl           = berendsen
tc-grps          = System   
tau_t             = 0.1

; temperature supposed to set to 20*C @ 293.15 K as reference work ; 
ref_t             = 293.15

Pcoupl           = no

gen_vel          = yes
gen_temp         = 298.15
gen_seed         = 173529

constraints         = none


; NEMD ;

acc_grps                  = system
accelerate                = 0.1 0.0 0.0
cos_acceleration          = 0.02


Any suggestion are most welcomed and thank you in advance!



MUHAMMAD ALIF MOHAMMAD LATIF
Laboratory of Theoretical and Computational Chemistry
Department of Chemistry
Faculty of Science
Universiti Putra Malaysia
43400 UPM Serdang, Selangor
MALAYSIA



      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110317/b3c6c94e/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: --static--liam_fetch_bl.gif
Type: image/gif
Size: 9852 bytes
Desc: not available
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110317/b3c6c94e/attachment.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: --static--liam_fetch_lb.gif
Type: image/gif
Size: 4545 bytes
Desc: not available
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20110317/b3c6c94e/attachment-0001.gif>