.TH gmx-rotmat 1 "" "VERSION 5.0.2" "GROMACS Manual"
.SH NAME
gmx-rotmat \- Plot the rotation matrix for fitting to a reference structure

.SH SYNOPSIS
gmx rotmat [\-f [<.xtc/.trr/...>]] [\-s [<.tpr/.tpb/...>]]
           [\-n [<.ndx>]] [\-o [<.xvg>]] [\-nice <int>] [\-b <time>]
           [\-e <time>] [\-dt <time>] [\-[no]w] [\-xvg <enum>]
           [\-ref <enum>] [\-skip <int>] [\-[no]fitxy] [\-[no]mw]

.SH DESCRIPTION
\fBgmx rotmat\fR plots the rotation matrix required for least squares fitting a conformation onto the reference conformation provided with \fB\-s\fR. Translation is removed before fitting. The output are the three vectors that give the new directions of the x, y and z directions of the reference conformation, for example: (zx,zy,zz) is the orientation of the reference z\-axis in the trajectory frame.

This tool is useful for, for instance, determining the orientation of a molecule at an interface, possibly on a trajectory produced with \fBgmx trjconv \-fit rotxy+transxy\fR to remove the rotation in the \fIx\-y\fR plane.

Option \fB\-ref\fR determines a reference structure for fitting, instead of using the structure from \fB\-s\fR. The structure with the lowest sum of RMSD's to all other structures is used. Since the computational cost of this procedure grows with the square of the number of frames, the \fB\-skip\fR option can be useful. A full fit or only a fit in the \fIx\-y\fR plane can be performed.

Option \fB\-fitxy\fR fits in the \fIx\-y\fR plane before determining the rotation matrix.

.SH OPTIONS
Options to specify input and output files:

.BI "\-f" " [<.xtc/.trr/...>] (traj.xtc) (Input)"
    Trajectory: xtc trr cpt trj gro g96 pdb tng

.BI "\-s" " [<.tpr/.tpb/...>] (topol.tpr) (Input)"
    Structure+mass(db): tpr tpb tpa gro g96 pdb brk ent

.BI "\-n" " [<.ndx>] (index.ndx) (Input, Optional)"
    Index file

.BI "\-o" " [<.xvg>] (rotmat.xvg) (Output)"
    xvgr/xmgr file


Other options:

.BI "\-nice" " <int> (19)"
    Set the nicelevel

.BI "\-b" " <time> (0)"
    First frame (ps) to read from trajectory

.BI "\-e" " <time> (0)"
    Last frame (ps) to read from trajectory

.BI "\-dt" " <time> (0)"
    Only use frame when t MOD dt = first time (ps)

.BI "\-[no]w" "  (no)"
    View output \fB.xvg\fR, \fB.xpm\fR, \fB.eps\fR and \fB.pdb\fR files

.BI "\-xvg" " <enum> (xmgrace)"
    xvg plot formatting: xmgrace, xmgr, none

.BI "\-ref" " <enum> (none)"
    Determine the optimal reference structure: none, xyz, xy

.BI "\-skip" " <int> (1)"
    Use every nr\-th frame for \fB\-ref\fR

.BI "\-[no]fitxy" "  (no)"
    Fit the x/y rotation before determining the rotation

.BI "\-[no]mw" "  (yes)"
    Use mass weighted fitting


.SH SEE ALSO
.BR gromacs(7)

More information about \fBGROMACS\fR is available at <\fIhttp://www.gromacs.org/\fR>.