Mark Abraham wrote:
Vitaly V. Chaban wrote:
Mark,
In "g_energy -h" I see nothing abous viscosity except the option
"-vis". The command line is as I wrote earlier:
"g_energy -vis -acflen 1001"
The system is liquid water. The simulation time is 2 000 ps. The file
I look into is enecorr.xvg.
What the error of some sort do you mean?
-acflen is measured in trajectory frames. How many relevant frames exist
in the trajectory file?
Viscosity in g_energy is computed from the pressure tensor which is
stored in the energy file. Unfortunately is seems like the acflen is
ignored in this case. Please submit a bugzilla.
On the other hand you need all the data you can get for this quantity,
because it converges extremely slowly.
Mark
Vitaly
When using g_energy to calculate viscosity ("g_energy -vis"), in
enecorr.xvg the ACF of stress tensor is output. How can one specify
the "acflen" of this ACF?
I try the key "g_energy -vis -acflen XXX" but the length of the ACF is
still equal to the half of the trajectory length.
If you read g_energy -h, you will need to specify the ACF length in a
particular way. If you'd given us your actual command line and some
information about your trajectory, you might have gotten a definitive
answer... as it is, we can only assume you made an error of some
sort! :-)
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the www
interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
--
David van der Spoel, Ph.D., Professor of Biology
Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone: +46184714205. Fax: +4618511755.
sp...@xray.bmc.uu.se sp...@gromacs.org http://folding.bmc.uu.se
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php