| VERSION 3.3_beta_20050823 |
Compute partial densities across the box, using an index file. Densities
in kg/m^3, number densities or electron densities can be
calculated. For electron densities, a file describing the number of
electrons for each type of atom should be provided using -ei.
It should look like:
2
atomname = nrelectrons
atomname = nrelectrons
The first line contains the number of lines to read from the file.
There should be one line for each unique atom name in your system.
The number of electrons for each atom is modified by its atomic
partial charge.
option | filename | type | description |
---|---|---|---|
-f | traj.xtc | Input | Generic trajectory: xtc trr trj gro g96 pdb |
-n | index.ndx | Input, Opt. | Index file |
-s | topol.tpr | Input | Generic run input: tpr tpb tpa xml |
-ei | electrons.dat | Input, Opt. | Generic data file |
-o | density.xvg | Output | xvgr/xmgr file |
option | type | default | description |
---|---|---|---|
-[no]h | bool | no | Print help info and quit |
-nice | int | 19 | Set the nicelevel |
-b | time | 0 | First frame (ps) to read from trajectory |
-e | time | 0 | Last frame (ps) to read from trajectory |
-dt | time | 0 | Only use frame when t MOD dt = first time (ps) |
-[no]w | bool | no | View output xvg, xpm, eps and pdb files |
-[no]xvgr | bool | yes | Add specific codes (legends etc.) in the output xvg files for the xmgrace program |
-d | string | Z | Take the normal on the membrane in direction X, Y or Z. |
-sl | int | 10 | Divide the box in #nr slices. |
-[no]number | bool | no | Calculate number density instead of mass density. Hydrogens are not counted! |
-[no]ed | bool | no | Calculate electron density instead of mass density |
-[no]count | bool | no | Only count atoms in slices, no densities. Hydrogens are not counted |
-ng | int | 0 | Number of groups to compute densities of |
-[no]symm | bool | no | Symmetrize the density along the axis, with respect to the center. Useful for bilayers. |
-[no]center | bool | no | Shift the center of mass along the axis to zero. This means if your axis is Z and your box is bX, bY, bZ, the center of mass will be at bX/2, bY/2, 0. |