1g_rotacf(1) GROMACS suite, VERSION 4.5 g_rotacf(1)
2
3
4
6 g_rotacf - calculates the rotational correlation function for molecules
7
8 VERSION 4.5
9
11 g_rotacf -f traj.xtc -s topol.tpr -n index.ndx -o rotacf.xvg -[no]h
12 -[no]version -nice int -b time -e time -dt time -[no]w -xvg enum -[no]d
13 -[no]aver -acflen int -[no]normalize -P enum -fitfn enum -ncskip int
14 -beginfit real -endfit real
15
17 g_rotacf calculates the rotational correlation function for molecules.
18 Three atoms (i,j,k) must be given in the index file, defining two vec‐
19 tors ij and jk. The rotational acf is calculated as the autocorrelation
20 function of the vector n = ij x jk, i.e. the cross product of the two
21 vectors. Since three atoms span a plane, the order of the three atoms
22 does not matter. Optionally, controlled by the -d switch, you can cal‐
23 culate the rotational correlation function for linear molecules by
24 specifying two atoms (i,j) in the index file.
25
26
27 EXAMPLES
28
29
30 g_rotacf -P 1 -nparm 2 -fft -n index -o rotacf-x-P1 -fa expfit-x-P1
31 -beginfit 2.5 -endfit 20.0
32
33
34 This will calculate the rotational correlation function using a first
35 order Legendre polynomial of the angle of a vector defined by the index
36 file. The correlation function will be fitted from 2.5 ps till 20.0 ps
37 to a two parameter exponential.
38
40 -f traj.xtc Input
41 Trajectory: xtc trr trj gro g96 pdb cpt
42
43 -s topol.tpr Input
44 Run input file: tpr tpb tpa
45
46 -n index.ndx Input
47 Index file
48
49 -o rotacf.xvg Output
50 xvgr/xmgr file
51
52
54 -[no]hno
55 Print help info and quit
56
57 -[no]versionno
58 Print version info and quit
59
60 -nice int 19
61 Set the nicelevel
62
63 -b time 0
64 First frame (ps) to read from trajectory
65
66 -e time 0
67 Last frame (ps) to read from trajectory
68
69 -dt time 0
70 Only use frame when t MOD dt = first time (ps)
71
72 -[no]wno
73 View output xvg, xpm, eps and pdb files
74
75 -xvg enum xmgrace
76 xvg plot formatting: xmgrace, xmgr or none
77
78 -[no]dno
79 Use index doublets (vectors) for correlation function instead of
80 triplets (planes)
81
82 -[no]averyes
83 Average over molecules
84
85 -acflen int -1
86 Length of the ACF, default is half the number of frames
87
88 -[no]normalizeyes
89 Normalize ACF
90
91 -P enum 0
92 Order of Legendre polynomial for ACF (0 indicates none): 0, 1, 2 or
93 3
94
95 -fitfn enum none
96 Fit function: none, exp, aexp, exp_exp, vac, exp5, exp7 or
97 exp9
98
99 -ncskip int 0
100 Skip N points in the output file of correlation functions
101
102 -beginfit real 0
103 Time where to begin the exponential fit of the correlation function
104
105 -endfit real -1
106 Time where to end the exponential fit of the correlation function, -1
107 is until the end
108
109
111 gromacs(7)
112
113 More information about GROMACS is available at <http://www.gro‐
114 macs.org/>.
115
116
117
118 Thu 26 Aug 2010 g_rotacf(1)