Hi AMBER folks,
I used NAMD to run a box of 901 molecules of water using TIP3P model. I used
tleap to generate a tip.top for 901 molecules
and using TIP3P model. I tried using ptraj to read the namd file using the
tip.top as parm file to compute the diffusion coefficient of
water at 298 K I get 3.01*10-5 cm2/sec. which is incorrect.
AT the moment I dont understand why I am getting the wrong number. Can
anyone please help me.
thanks,
ganesh
*********details*****************
[ganesh.willis water]$ ~/ptraj10 tip.top
\-/
-/- PTRAJ: a utility for processing trajectory files
/-\
\-/ Version: "AMBER 10.0 integrated" (2/15/2008)
-/- Executable is: "/home/ganesh/ptraj10"
/-\
\-/ Residue labels:
WAT WAT WAT WAT WAT WAT WAT WAT WAT WAT
...
WAT
PRMTOP does not contain %FLAG RADIUS_SET
PTRAJ: Processing input from "STDIN" ...
trajin equil3.dcd
PTRAJ: trajin equil3.dcd
Checking coordinates: equil3.dcd
diffusion :1-901 1.0 average diffusion
PTRAJ: diffusion :1-901 1.0 average diffusion
Mask [:1-901] represents 2703 atoms
go
[No output trajectory specified (trajout)]
PTRAJ: Successfully read the input file.
Coordinate processing will occur on 100 frames.
Summary of I/O and actions follows:
INPUT COORDINATE FILES
File (equil3.dcd) is a CHARMM trajectory in big endian binary format
with box information representing 100 sets
NFILE = 100 ISTEP = 1000 NINTV = 1000 NSAVC = 100000 NSAVV
= 0
NDEGF = 0 NFREA = 0 DELTA = 0.0205 QCRYS = 1 QDIM4
= 0
VERNU = 24 NATREC= 2703 NFREAT = 2703
Dumping the title:
REMARKS FILENAME=equil3.dcd CREATED BY
NAMD
REMARKS DATE: 01/21/11 CREATED BY USER:
jpotoff
NO OUTPUT COORDINATE FILE WAS SPECIFIED
ACTIONS
1> DIFFUSION
Only the average results will be dumped to diffusion_?.xmgr
The time between frames in psec is 1.000.
To calculated diffusion constants, calculate the slope of the lines(s)
and multiply by 10.0/6.0; this will give units of 1x10**-5 cm**2/s
The atoms in the calculation follow: * (All atoms are selected)
Processing CHARMM trajectory file equil3.dcd
Set 1 .................................................
Set 50 .................................................
Set 100
PTRAJ: Successfully read in 101 sets and processed 100 sets.
Dumping accumulated results (if any)
[ganesh.willis water]$ xmgrace diffusion_r.xmgr &
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
- application/octet-stream attachment: tip.top
Received on Fri Jan 21 2011 - 13:00:03 PST