Synopsis
gmx potential [-f [<.xtc/.trr/...>]] [-n [<.ndx>]] [-s [<.tpr/.tpb/...>]]
[-o [<.xvg>]] [-oc [<.xvg>]] [-of [<.xvg>]] [-nice ]
[-b ] [-e ] [-dt ] [-[no]w] [-xvg ]
[-d ] [-sl ] [-cb ] [-ce ] [-tz ]
[-[no]spherical] [-ng ] [-[no]correct]
Description
gmx potential computes the electrostatical potential across the box. The potential is calculated by first summing the charges per slice and then integrating twice of this charge distribution. Periodic boundaries are not taken into account. Reference of potential is taken to be the left side of the box. It is also possible to calculate the potential in spherical coordinates as function of r by calculating a charge distribution in spherical slices and twice integrating them. epsilon_r is taken as 1, but 2 is more appropriate in many cases.
Options
Options to specify input and output files:
-f xtc /.trr /...>] (traj.xtc ) (Input)
Trajectory: xtc trr cpt trj gro g96 pdb tng
-n ndx >] (index.ndx ) (Input)
Index file
-s tpr /.tpb /...>] (topol.tpr ) (Input)
Run input file: tpr tpb tpa
-o xvg >] (potential.xvg ) (Output)
xvgr/xmgr file
-oc xvg >] (charge.xvg ) (Output)
xvgr/xmgr file
-of xvg >] (field.xvg ) (Output)
xvgr/xmgr file
Other options:
-nice
Set the nicelevel
-b
First frame (ps) to read from trajectory
-e
Last frame (ps) to read from trajectory
-dt
Only use frame when t MOD dt = first time (ps)
-[no]w
View output .xvg , .xpm , .eps and .pdb files
-xvg
xvg plot formatting: xmgrace, xmgr, none
-d
Take the normal on the membrane in direction X, Y or Z.
-sl
Calculate potential as function of boxlength, dividing the box in this number of slices.
-cb
Discard this number of first slices of box for integration
-ce
Discard this number of last slices of box for integration
-tz
Translate all coordinates by this distance in the direction of the box
-[no]spherical
Calculate spherical thingie
-ng
Number of groups to consider
-[no]correct
Assume net zero charge of groups to improve accuracy
Known Issues
Discarding slices for integration should not be necessary.
