g_msd_d

Langue: en

Autres versions - même langue

Version: 264864 (debian - 07/07/09)

Section: 1 (Commandes utilisateur)

NAME

g_msd - calculates mean square displacements

SYNOPSIS

g_msd -f traj.xtc -s topol.tpr -n index.ndx -o msd.xvg -mol diff_mol.xvg -[no]h -nice int -b time -e time -dt time -tu enum -[no]w -[no]xvgr -type enum -lateral enum -ngroup int -[no]mw -trestart time -beginfit time -endfit time

DESCRIPTION

g_msd computes the mean square displacement (MSD) of atoms from their initial positions. This provides an easy way to compute the diffusion constant using the Einstein relation. The time between additional starting points for the MSD calculation is set with -trestart The diffusion constant is calculated by least squares fitting a straight line through the MSD from -beginfit to

-endfit . An error estimate given, which is the difference of the diffusion coefficients obtained from fits over the two halfs of the fit interval.

Option -mol plots the MSD for molecules, this implies

-mw , i.e. for each inidividual molecule an diffusion constant is computed. When using an index file, it should contain molecule numbers instead of atom numbers. Using this option one also gets an accurate error estimate based on the statistics between individual molecules. Since one usually is interested in self-diffusion at infinite dilution this is probably the most useful number.

FILES

-f traj.xtc Input
 Generic trajectory: xtc trr trj gro g96 pdb 

-s topol.tpr Input
 Structure+mass(db): tpr tpb tpa gro g96 pdb xml 

-n index.ndx Input, Opt.
 Index file 

-o msd.xvg Output
 xvgr/xmgr file 

-mol diff_mol.xvg Output, Opt.
 xvgr/xmgr file 

OTHER OPTIONS

-[no]h no
 Print help info and quit

-nice int 19
 Set the nicelevel

-b time 0
 First frame (ps) to read from trajectory

-e time 0
 Last frame (ps) to read from trajectory

-dt time 0
 Only use frame when t MOD dt = first time (ps)

-tu enum ps
 Time unit:  ps , fs , ns , us , ms , s , m or h

-[no]w no
 View output xvg, xpm, eps and pdb files

-[no]xvgr yes
 Add specific codes (legends etc.) in the output xvg files for the xmgrace program

-type enum no
 Compute diffusion coefficient in one direction:  no , x , y or z

-lateral enum no
 Calculate the lateral diffusion in a plane perpendicular to:  no , x , y or z

-ngroup int 1
 Number of groups to calculate MSD for

-[no]mw yes
 Mass weighted MSD

-trestart time 10
 Time between restarting points in trajectory (ps)

-beginfit time -1
 Start time for fitting the MSD (ps), -1 is 10%

-endfit time -1
 End time for fitting the MSD (ps), -1 is 90%

SEE ALSO

gromacs(7)

More information about the GROMACS suite is available in /usr/share/doc/gromacs or at <http://www.gromacs.org/>.