| VERSION 4.5 |
g_covar calculates and diagonalizes the (mass-weighted) covariance matrix. All structures are fitted to the structure in the structure file. When this is not a run input file periodicity will not be taken into account. When the fit and analysis groups are identical and the analysis is non mass-weighted, the fit will also be non mass-weighted.
The eigenvectors are written to a trajectory file (-v). When the same atoms are used for the fit and the covariance analysis, the reference structure for the fit is written first with t=-1. The average (or reference when -ref is used) structure is written with t=0, the eigenvectors are written as frames with the eigenvector number as timestamp.
The eigenvectors can be analyzed with g_anaeig.
Option -ascii writes the whole covariance matrix to an ASCII file. The order of the elements is: x1x1, x1y1, x1z1, x1x2, ...
Option -xpm writes the whole covariance matrix to an xpm file.
Option -xpma writes the atomic covariance matrix to an xpm file, i.e. for each atom pair the sum of the xx, yy and zz covariances is written.
option | filename | type | description |
---|---|---|---|
-f | traj.xtc | Input | Trajectory: xtc trr trj gro g96 pdb cpt |
-s | topol.tpr | Input | Structure+mass(db): tpr tpb tpa gro g96 pdb |
-n | index.ndx | Input, Opt. | Index file |
-o | eigenval.xvg | Output | xvgr/xmgr file |
-v | eigenvec.trr | Output | Full precision trajectory: trr trj cpt |
-av | average.pdb | Output | Structure file: gro g96 pdb etc. |
-l | covar.log | Output | Log file |
-ascii | covar.dat | Output, Opt. | Generic data file |
-xpm | covar.xpm | Output, Opt. | X PixMap compatible matrix file |
-xpma | covara.xpm | Output, Opt. | X PixMap compatible matrix file |
option | type | default | description |
---|---|---|---|
-[no]h | gmx_bool | no | Print help info and quit |
-[no]version | gmx_bool | no | Print version 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) |
-tu | enum | ps | Time unit: fs, ps, ns, us, ms or s |
-xvg | enum | xmgrace | xvg plot formatting: xmgrace, xmgr or none |
-[no]fit | gmx_bool | yes | Fit to a reference structure |
-[no]ref | gmx_bool | no | Use the deviation from the conformation in the structure file instead of from the average |
-[no]mwa | gmx_bool | no | Mass-weighted covariance analysis |
-last | int | -1 | Last eigenvector to write away (-1 is till the last) |
-[no]pbc | gmx_bool | yes | Apply corrections for periodic boundary conditions |