1GMX-VELACC(1)                       GROMACS                      GMX-VELACC(1)
2
3
4

NAME

6       gmx-velacc - Calculate velocity autocorrelation functions
7

SYNOPSIS

9          gmx velacc [-f [<.trr/.cpt/...>]] [-s [<.tpr/.gro/...>]] [-n [<.ndx>]]
10                     [-o [<.xvg>]] [-os [<.xvg>]] [-b <time>] [-e <time>]
11                     [-dt <time>] [-[no]w] [-xvg <enum>] [-[no]m] [-[no]recip]
12                     [-[no]mol] [-acflen <int>] [-[no]normalize] [-P <enum>]
13                     [-fitfn <enum>] [-beginfit <real>] [-endfit <real>]
14

DESCRIPTION

16       gmx velacc computes the velocity autocorrelation function.  When the -m
17       option is used, the momentum autocorrelation function is calculated.
18
19       With option -mol the velocity autocorrelation function of molecules  is
20       calculated.  In  this  case  the index group should consist of molecule
21       numbers instead of atom numbers.
22
23       By using option -os you can also extract  the  estimated  (vibrational)
24       power spectrum, which is the Fourier transform of the velocity autocor‐
25       relation function.  Be sure that your trajectory contains  frames  with
26       velocity information (i.e. nstvout was set in your original .mdp file),
27       and that the time interval  between  data  collection  points  is  much
28       shorter than the time scale of the autocorrelation.
29

OPTIONS

31       Options to specify input files:
32
33       -f [<.trr/.cpt/…>] (traj.trr)
34              Full precision trajectory: trr cpt tng
35
36       -s [<.tpr/.gro/…>] (topol.tpr) (Optional)
37              Structure+mass(db): tpr gro g96 pdb brk ent
38
39       -n [<.ndx>] (index.ndx) (Optional)
40              Index file
41
42       Options to specify output files:
43
44       -o [<.xvg>] (vac.xvg)
45              xvgr/xmgr file
46
47       -os [<.xvg>] (spectrum.xvg) (Optional)
48              xvgr/xmgr file
49
50       Other options:
51
52       -b <time> (0)
53              Time of first frame to read from trajectory (default unit ps)
54
55       -e <time> (0)
56              Time of last frame to read from trajectory (default unit ps)
57
58       -dt <time> (0)
59              Only use frame when t MOD dt = first time (default unit ps)
60
61       -[no]w (no)
62              View output .xvg, .xpm, .eps and .pdb files
63
64       -xvg <enum> (xmgrace)
65              xvg plot formatting: xmgrace, xmgr, none
66
67       -[no]m (no)
68              Calculate the momentum autocorrelation function
69
70       -[no]recip (yes)
71              Use cm^-1 on X-axis instead of 1/ps for spectra.
72
73       -[no]mol (no)
74              Calculate the velocity acf of molecules
75
76       -acflen <int> (-1)
77              Length of the ACF, default is half the number of frames
78
79       -[no]normalize (yes)
80              Normalize ACF
81
82       -P <enum> (0)
83              Order  of  Legendre polynomial for ACF (0 indicates none): 0, 1,
84              2, 3
85
86       -fitfn <enum> (none)
87              Fit function: none, exp, aexp, exp_exp, exp5, exp7, exp9
88
89       -beginfit <real> (0)
90              Time where to begin the exponential fit of the correlation func‐
91              tion
92
93       -endfit <real> (-1)
94              Time  where  to end the exponential fit of the correlation func‐
95              tion, -1 is until the end
96

SEE ALSO

98       gmx(1)
99
100       More    information    about    GROMACS    is    available    at     <‐
101       http://www.gromacs.org/>.
102
104       2019, GROMACS development team
105
106
107
108
1092019.2                           Apr 16, 2019                    GMX-VELACC(1)
Impressum