1g_polystat(1) GROMACS suite, VERSION 4.5 g_polystat(1)
2
3
4
6 g_polystat - calculates static properties of polymers
7
8 VERSION 4.5
9
11 g_polystat -s topol.tpr -f traj.xtc -n index.ndx -o polystat.xvg -v
12 polyvec.xvg -p persist.xvg -i intdist.xvg -[no]h -[no]version -nice int
13 -b time -e time -dt time -tu enum -[no]w -xvg enum -[no]mw -[no]pc
14
16 g_polystat plots static properties of polymers as a function of time
17 and prints the average.
18
19
20 By default it determines the average end-to-end distance and radii of
21 gyration of polymers. It asks for an index group and split this into
22 molecules. The end-to-end distance is then determined using the first
23 and the last atom in the index group for each molecules. For the
24 radius of gyration the total and the three principal components for the
25 average gyration tensor are written. With option -v the eigenvectors
26 are written. With option -pc also the average eigenvalues of the
27 individual gyration tensors are written. With option -i the mean
28 square internal distances are written.
29
30
31 With option -p the persistence length is determined. The chosen index
32 group should consist of atoms that are consecutively bonded in the
33 polymer mainchains. The persistence length is then determined from the
34 cosine of the angles between bonds with an index difference that is
35 even, the odd pairs are not used, because straight polymer backbones
36 are usually all trans and therefore only every second bond aligns. The
37 persistence length is defined as number of bonds where the average cos
38 reaches a value of 1/e. This point is determined by a linear interpola‐
39 tion of log(cos).
40
42 -s topol.tpr Input
43 Run input file: tpr tpb tpa
44
45 -f traj.xtc Input
46 Trajectory: xtc trr trj gro g96 pdb cpt
47
48 -n index.ndx Input, Opt.
49 Index file
50
51 -o polystat.xvg Output
52 xvgr/xmgr file
53
54 -v polyvec.xvg Output, Opt.
55 xvgr/xmgr file
56
57 -p persist.xvg Output, Opt.
58 xvgr/xmgr file
59
60 -i intdist.xvg Output, Opt.
61 xvgr/xmgr file
62
63
65 -[no]hno
66 Print help info and quit
67
68 -[no]versionno
69 Print version info and quit
70
71 -nice int 19
72 Set the nicelevel
73
74 -b time 0
75 First frame (ps) to read from trajectory
76
77 -e time 0
78 Last frame (ps) to read from trajectory
79
80 -dt time 0
81 Only use frame when t MOD dt = first time (ps)
82
83 -tu enum ps
84 Time unit: fs, ps, ns, us, ms or s
85
86 -[no]wno
87 View output xvg, xpm, eps and pdb files
88
89 -xvg enum xmgrace
90 xvg plot formatting: xmgrace, xmgr or none
91
92 -[no]mwyes
93 Use the mass weighting for radii of gyration
94
95 -[no]pcno
96 Plot average eigenvalues
97
98
100 gromacs(7)
101
102 More information about GROMACS is available at <http://www.gro‐
103 macs.org/>.
104
105
106
107 Thu 26 Aug 2010 g_polystat(1)