1g_sham(1)                 GROMACS suite, VERSION 4.5                 g_sham(1)
2
3
4

NAME

6       g_sham - read/write xmgr and xvgr data sets
7
8       VERSION 4.5
9

SYNOPSIS

11       g_sham  -f graph.xvg -ge gibbs.xvg -ene esham.xvg -dist ener.xvg -histo
12       edist.xvg -bin bindex.ndx -lp prob.xpm -ls gibbs.xpm -lsh  enthalpy.xpm
13       -lss  entropy.xpm  -map  map.xpm  -ls3 gibbs3.pdb -mdata mapdata.xvg -g
14       shamlog.log -[no]h -[no]version -nice int -[no]w -xvg enum -[no]time -b
15       real  -e  real  -ttol real -n int -[no]d -bw real -[no]sham -tsham real
16       -pmin real -dim vector -ngrid vector -xmin vector  -xmax  vector  -pmax
17       real -gmax real -emin real -emax real -nlevels int -mname string
18

DESCRIPTION

20       g_sham makes multi-dimensional free-energy, enthalpy and entropy plots.
21       g_sham reads one or more xvg files  and  analyzes  data  sets.   g_sham
22       basic purpose is plotting Gibbs free energy landscapes (option  -ls) by
23       Bolzmann inverting multi-dimensional histograms (option   -lp)  but  it
24       can also make enthalpy (option  -lsh) and entropy (option  -lss) plots.
25       The histograms can be made for any quantities  the  user  supplies.   A
26       line  in  the  input file may start with a time (see option  -time) and
27       any number of y values may follow.  Multiple sets can also be read when
28       they  are separated by & (option  -n), in this case only one y value is
29       read from each line.  All lines starting with  and @ are skipped.
30
31
32       Option  -ge can be used to supply a file with free  energies  when  the
33       ensemble  is  not  a Boltzmann ensemble, but needs to be biased by this
34       free energy. One free energy value is required for  each  (multi-dimen‐
35       sional) data point in the  -f input.
36
37
38       Option   -ene  can be used to supply a file with energies.  These ener‐
39       gies are used as a weighting function in the single histogram  analysis
40       method  due  to Kumar et. al. When also temperatures are supplied (as a
41       second column in the file) an experimental weighting scheme is applied.
42       In addition the vales are used for making enthalpy and entropy plots.
43
44
45       With  option   -dim dimensions can be gives for distances.  When a dis‐
46       tance is 2- or 3-dimensional, the circumference or surface  sampled  by
47       two  particles  increases  with increasing distance.  Depending on what
48       one would like to show, one can choose to  correct  the  histogram  and
49       free-energy for this volume effect.  The probability is normalized by r
50       and r2 for a dimension of 2 and 3 respectively.  A value of -1 is  used
51       to  indicate an angle in degrees between two vectors: a sin(angle) nor‐
52       malization will be applied.  Note that for angles between  vectors  the
53       inner-product or cosine is the natural quantity to use, as it will pro‐
54       duce bins of the same volume.
55

FILES

57       -f graph.xvg Input
58        xvgr/xmgr file
59
60       -ge gibbs.xvg Input, Opt.
61        xvgr/xmgr file
62
63       -ene esham.xvg Input, Opt.
64        xvgr/xmgr file
65
66       -dist ener.xvg Output, Opt.
67        xvgr/xmgr file
68
69       -histo edist.xvg Output, Opt.
70        xvgr/xmgr file
71
72       -bin bindex.ndx Output, Opt.
73        Index file
74
75       -lp prob.xpm Output, Opt.
76        X PixMap compatible matrix file
77
78       -ls gibbs.xpm Output, Opt.
79        X PixMap compatible matrix file
80
81       -lsh enthalpy.xpm Output, Opt.
82        X PixMap compatible matrix file
83
84       -lss entropy.xpm Output, Opt.
85        X PixMap compatible matrix file
86
87       -map map.xpm Output, Opt.
88        X PixMap compatible matrix file
89
90       -ls3 gibbs3.pdb Output, Opt.
91        Protein data bank file
92
93       -mdata mapdata.xvg Output, Opt.
94        xvgr/xmgr file
95
96       -g shamlog.log Output, Opt.
97        Log file
98
99

OTHER OPTIONS

101       -[no]hno
102        Print help info and quit
103
104       -[no]versionno
105        Print version info and quit
106
107       -nice int 19
108        Set the nicelevel
109
110       -[no]wno
111        View output xvg, xpm, eps and pdb files
112
113       -xvg enum xmgrace
114        xvg plot formatting:  xmgrace,  xmgr or  none
115
116       -[no]timeyes
117        Expect a time in the input
118
119       -b real -1
120        First time to read from set
121
122       -e real -1
123        Last time to read from set
124
125       -ttol real 0
126        Tolerance on time in appropriate units (usually ps)
127
128       -n int 1
129        Read  sets separated by &
130
131       -[no]dno
132        Use the derivative
133
134       -bw real 0.1
135        Binwidth for the distribution
136
137       -[no]shamyes
138        Turn off energy weighting even if energies are given
139
140       -tsham real 298.15
141        Temperature for single histogram analysis
142
143       -pmin real 0
144        Minimum probability. Anything lower than this will be set to zero
145
146       -dim vector 1 1 1
147        Dimensions for distances, used for volume correction  (max  3  values,
148       dimensions  3 will get the same value as the last)
149
150       -ngrid vector 32 32 32
151        Number of bins for energy landscapes (max 3 values, dimensions  3 will
152       get the same value as the last)
153
154       -xmin vector 0 0 0
155        Minimum for the axes in energy landscape (see above for  3 dimensions)
156
157       -xmax vector 1 1 1
158        Maximum for the axes in energy landscape (see above for  3 dimensions)
159
160       -pmax real 0
161        Maximum probability in output, default is calculate
162
163       -gmax real 0
164        Maximum free energy in output, default is calculate
165
166       -emin real 0
167        Minimum enthalpy in output, default is calculate
168
169       -emax real 0
170        Maximum enthalpy in output, default is calculate
171
172       -nlevels int 25
173        Number of levels for energy landscape
174
175       -mname string
176        Legend label for the custom landscape
177
178

SEE ALSO

180       gromacs(7)
181
182       More  information  about  GROMACS  is  available  at   <http://www.gro
183       macs.org/>.
184
185
186
187                                Thu 26 Aug 2010                      g_sham(1)
Impressum