.TH g_anaeig 1 "Mon 4 Apr 2011" "" "GROMACS suite, VERSION 4.5.4-dev-20110404-bc5695c"
.SH NAME
g_anaeig - analyzes the eigenvectors

.B VERSION 4.5.4-dev-20110404-bc5695c
.SH SYNOPSIS
\f3g_anaeig\fP
.BI "\-v" " eigenvec.trr "
.BI "\-v2" " eigenvec2.trr "
.BI "\-f" " traj.xtc "
.BI "\-s" " topol.tpr "
.BI "\-n" " index.ndx "
.BI "\-eig" " eigenval.xvg "
.BI "\-eig2" " eigenval2.xvg "
.BI "\-comp" " eigcomp.xvg "
.BI "\-rmsf" " eigrmsf.xvg "
.BI "\-proj" " proj.xvg "
.BI "\-2d" " 2dproj.xvg "
.BI "\-3d" " 3dproj.pdb "
.BI "\-filt" " filtered.xtc "
.BI "\-extr" " extreme.pdb "
.BI "\-over" " overlap.xvg "
.BI "\-inpr" " inprod.xpm "
.BI "\-[no]h" ""
.BI "\-[no]version" ""
.BI "\-nice" " int "
.BI "\-b" " time "
.BI "\-e" " time "
.BI "\-dt" " time "
.BI "\-tu" " enum "
.BI "\-[no]w" ""
.BI "\-xvg" " enum "
.BI "\-first" " int "
.BI "\-last" " int "
.BI "\-skip" " int "
.BI "\-max" " real "
.BI "\-nframes" " int "
.BI "\-[no]split" ""
.BI "\-[no]entropy" ""
.BI "\-temp" " real "
.BI "\-nevskip" " int "
.SH DESCRIPTION
\&\fB g_anaeig\fR analyzes eigenvectors. The eigenvectors can be of a
\&covariance matrix (\fB g_covar\fR) or of a Normal Modes analysis
\&(\fB g_nmeig\fR).


\&When a trajectory is projected on eigenvectors, all structures are
\&fitted to the structure in the eigenvector file, if present, otherwise
\&to the structure in the structure file. When no run input file is
\&supplied, periodicity will not be taken into account. Most analyses
\&are performed on eigenvectors \fB \-first\fR to \fB \-last\fR, but when
\&\fB \-first\fR is set to \-1 you will be prompted for a selection.


\&\fB \-comp\fR: plot the vector components per atom of eigenvectors
\&\fB \-first\fR to \fB \-last\fR.


\&\fB \-rmsf\fR: plot the RMS fluctuation per atom of eigenvectors
\&\fB \-first\fR to \fB \-last\fR (requires \fB \-eig\fR).


\&\fB \-proj\fR: calculate projections of a trajectory on eigenvectors
\&\fB \-first\fR to \fB \-last\fR.
\&The projections of a trajectory on the eigenvectors of its
\&covariance matrix are called principal components (pc's).
\&It is often useful to check the cosine content of the pc's,
\&since the pc's of random diffusion are cosines with the number
\&of periods equal to half the pc index.
\&The cosine content of the pc's can be calculated with the program
\&\fB g_analyze\fR.


\&\fB \-2d\fR: calculate a 2d projection of a trajectory on eigenvectors
\&\fB \-first\fR and \fB \-last\fR.


\&\fB \-3d\fR: calculate a 3d projection of a trajectory on the first
\&three selected eigenvectors.


\&\fB \-filt\fR: filter the trajectory to show only the motion along
\&eigenvectors \fB \-first\fR to \fB \-last\fR.


\&\fB \-extr\fR: calculate the two extreme projections along a trajectory
\&on the average structure and interpolate \fB \-nframes\fR frames
\&between them, or set your own extremes with \fB \-max\fR. The
\&eigenvector \fB \-first\fR will be written unless \fB \-first\fR and
\&\fB \-last\fR have been set explicitly, in which case all eigenvectors
\&will be written to separate files. Chain identifiers will be added
\&when writing a \fB .pdb\fR file with two or three structures (you
\&can use \fB rasmol \-nmrpdb\fR to view such a \fB .pdb\fR file).


\&  Overlap calculations between covariance analysis:

\&  \fB Note:\fR the analysis should use the same fitting structure


\&\fB \-over\fR: calculate the subspace overlap of the eigenvectors in
\&file \fB \-v2\fR with eigenvectors \fB \-first\fR to \fB \-last\fR
\&in file \fB \-v\fR.


\&\fB \-inpr\fR: calculate a matrix of inner\-products between
\&eigenvectors in files \fB \-v\fR and \fB \-v2\fR. All eigenvectors
\&of both files will be used unless \fB \-first\fR and \fB \-last\fR
\&have been set explicitly.


\&When \fB \-v\fR, \fB \-eig\fR, \fB \-v2\fR and \fB \-eig2\fR are given,
\&a single number for the overlap between the covariance matrices is
\&generated. The formulas are:

\&        difference = sqrt(tr((sqrt(M1) \- sqrt(M2))2))

\&normalized overlap = 1 \- difference/sqrt(tr(M1) + tr(M2))

\&     shape overlap = 1 \- sqrt(tr((sqrt(M1/tr(M1)) \- sqrt(M2/tr(M2)))2))

\&where M1 and M2 are the two covariance matrices and tr is the trace
\&of a matrix. The numbers are proportional to the overlap of the square
\&root of the fluctuations. The normalized overlap is the most useful
\&number, it is 1 for identical matrices and 0 when the sampled
\&subspaces are orthogonal.


\&When the \fB \-entropy\fR flag is given an entropy estimate will be
\&computed based on the Quasiharmonic approach and based on
\&Schlitter's formula.
.SH FILES
.BI "\-v" " eigenvec.trr" 
.B Input
 Full precision trajectory: trr trj cpt 

.BI "\-v2" " eigenvec2.trr" 
.B Input, Opt.
 Full precision trajectory: trr trj cpt 

.BI "\-f" " traj.xtc" 
.B Input, Opt.
 Trajectory: xtc trr trj gro g96 pdb cpt 

.BI "\-s" " topol.tpr" 
.B Input, Opt.
 Structure+mass(db): tpr tpb tpa gro g96 pdb 

.BI "\-n" " index.ndx" 
.B Input, Opt.
 Index file 

.BI "\-eig" " eigenval.xvg" 
.B Input, Opt.
 xvgr/xmgr file 

.BI "\-eig2" " eigenval2.xvg" 
.B Input, Opt.
 xvgr/xmgr file 

.BI "\-comp" " eigcomp.xvg" 
.B Output, Opt.
 xvgr/xmgr file 

.BI "\-rmsf" " eigrmsf.xvg" 
.B Output, Opt.
 xvgr/xmgr file 

.BI "\-proj" " proj.xvg" 
.B Output, Opt.
 xvgr/xmgr file 

.BI "\-2d" " 2dproj.xvg" 
.B Output, Opt.
 xvgr/xmgr file 

.BI "\-3d" " 3dproj.pdb" 
.B Output, Opt.
 Structure file: gro g96 pdb etc. 

.BI "\-filt" " filtered.xtc" 
.B Output, Opt.
 Trajectory: xtc trr trj gro g96 pdb cpt 

.BI "\-extr" " extreme.pdb" 
.B Output, Opt.
 Trajectory: xtc trr trj gro g96 pdb cpt 

.BI "\-over" " overlap.xvg" 
.B Output, Opt.
 xvgr/xmgr file 

.BI "\-inpr" " inprod.xpm" 
.B Output, Opt.
 X PixMap compatible matrix file 

.SH OTHER OPTIONS
.BI "\-[no]h"  "no    "
 Print help info and quit

.BI "\-[no]version"  "no    "
 Print version info and quit

.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 "\-tu"  " enum" " ps" 
 Time unit: \fB fs\fR, \fB ps\fR, \fB ns\fR, \fB us\fR, \fB ms\fR or \fB s\fR

.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: \fB xmgrace\fR, \fB xmgr\fR or \fB none\fR

.BI "\-first"  " int" " 1" 
 First eigenvector for analysis (\-1 is select)

.BI "\-last"  " int" " 8" 
 Last eigenvector for analysis (\-1 is till the last)

.BI "\-skip"  " int" " 1" 
 Only analyse every nr\-th frame

.BI "\-max"  " real" " 0     " 
 Maximum for projection of the eigenvector on the average structure, max=0 gives the extremes

.BI "\-nframes"  " int" " 2" 
 Number of frames for the extremes output

.BI "\-[no]split"  "no    "
 Split eigenvector projections where time is zero

.BI "\-[no]entropy"  "no    "
 Compute entropy according to the Quasiharmonic formula or Schlitter's method.

.BI "\-temp"  " real" " 298.15" 
 Temperature for entropy calculations

.BI "\-nevskip"  " int" " 6" 
 Number of eigenvalues to skip when computing the entropy due to the quasi harmonic approximation. When you do a rotational and/or translational fit prior to the covariance analysis, you get 3 or 6 eigenvalues that are very close to zero, and which should not be taken into account when computing the entropy.

.SH SEE ALSO
.BR gromacs(7)

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