share/html/online/editconf.html

   1 <HTML>
   2 <HEAD>
   3 <TITLE>editconf</TITLE>
   4 <LINK rel=stylesheet href="style.css" type="text/css">
   5 <BODY text="#000000" bgcolor="#FFFFFF" link="#0000FF" vlink="#990000" alink="#FF0000">
   6 <TABLE WIDTH="98%" NOBORDER >
   7 <TR><TD WIDTH=400>
   8 <TABLE WIDTH=400 NOBORDER>
   9 <TD WIDTH=116>
  10 <a href="http://www.gromacs.org/"><img SRC="../images/gmxlogo_small.png"BORDER=0 </a></td>
  11 <td ALIGN=LEFT VALIGN=TOP WIDTH=280><br><h2>editconf</h2><font size=-1><A HREF="../online.html">Main Table of Contents</A></font><br><br></td>
  12 </TABLE></TD><TD WIDTH="*" ALIGN=RIGHT VALIGN=BOTTOM><p><B>VERSION 4.6<br>
  13 Sat 19 Jan 2013</B></td></tr></TABLE>
  14 <HR>
  15 <H3>Description</H3>
  16 <p>
  17 <tt>editconf</tt> converts generic structure format to <tt>.<a href="gro.html">gro</a></tt>, <tt>.<a href="g96.html">g96</a></tt>
  18 or <tt>.<a href="pdb.html">pdb</a></tt>.
  19 <p>
  20 The box can be modified with options <tt>-box</tt>, <tt>-d</tt> and
  21 <tt>-angles</tt>. Both <tt>-box</tt> and <tt>-d</tt>
  22 will center the system in the box, unless <tt>-noc</tt> is used.
  23 <p>
  24 Option <tt>-bt</tt> determines the box type: <tt>triclinic</tt> is a
  25 triclinic box, <tt>cubic</tt> is a rectangular box with all sides equal
  26 <tt>dodecahedron</tt> represents a rhombic dodecahedron and
  27 <tt>octahedron</tt> is a truncated octahedron.
  28 The last two are special cases of a triclinic box.
  29 The length of the three box vectors of the truncated octahedron is the
  30 shortest distance between two opposite hexagons.
  31 Relative to a cubic box with some periodic image distance, the volume of a
  32 dodecahedron with this same periodic distance is 0.71 times that of the cube,
  33 and that of a truncated octahedron is 0.77 times.
  34 <p>
  35 Option <tt>-box</tt> requires only
  36 one value for a cubic, rhombic dodecahedral, or truncated octahedral box.
  37 <p>
  38 With <tt>-d</tt> and a <tt>triclinic</tt> box the size of the system in the <it>x</it>-, <it>y</it>-,
  39 and <it>z</it>-directions is used. With <tt>-d</tt> and <tt>cubic</tt>,
  40 <tt>dodecahedron</tt> or <tt>octahedron</tt> boxes, the dimensions are set
  41 to the diameter of the system (largest distance between atoms) plus twice
  42 the specified distance.
  43 <p>
  44 Option <tt>-angles</tt> is only meaningful with option <tt>-box</tt> and
  45 a triclinic box and cannot be used with option <tt>-d</tt>.
  46 <p>
  47 When <tt>-n</tt> or <tt>-ndef</tt> is set, a group
  48 can be selected for calculating the size and the geometric center,
  49 otherwise the whole system is used.
  50 <p>
  51 <tt>-rotate</tt> rotates the coordinates and velocities.
  52 <p>
  53 <tt>-princ</tt> aligns the principal axes of the system along the
  54 coordinate axes, with the longest axis aligned with the <it>x</it>-axis.
  55 This may allow you to decrease the box volume,
  56 but beware that molecules can rotate significantly in a nanosecond.
  57 <p>
  58 Scaling is applied before any of the other operations are
  59 performed. Boxes and coordinates can be scaled to give a certain density (option
  60 <tt>-density</tt>). Note that this may be inaccurate in case a <tt>.<a href="gro.html">gro</a></tt>
  61 file is given as input. A special feature of the scaling option is that when the
  62 factor -1 is given in one dimension, one obtains a mirror image,
  63 mirrored in one of the planes. When one uses -1 in three dimensions,
  64 a point-mirror image is obtained.<p>
  65 Groups are selected after all operations have been applied.<p>
  66 Periodicity can be removed in a crude manner.
  67 It is important that the box vectors at the bottom of your input file
  68 are correct when the periodicity is to be removed.
  69 <p>
  70 When writing <tt>.<a href="pdb.html">pdb</a></tt> files, B-factors can be
  71 added with the <tt>-bf</tt> option. B-factors are read
  72 from a file with with following format: first line states number of
  73 entries in the file, next lines state an index
  74 followed by a B-factor. The B-factors will be attached per residue
  75 unless an index is larger than the number of residues or unless the
  76 <tt>-atom</tt> option is set. Obviously, any type of numeric data can
  77 be added instead of B-factors. <tt>-legend</tt> will produce
  78 a row of CA atoms with B-factors ranging from the minimum to the
  79 maximum value found, effectively making a legend for viewing.
  80 <p>
  81 With the option <tt>-mead</tt> a special <tt>.<a href="pdb.html">pdb</a></tt> (<tt>.pqr</tt>)
  82 file for the MEAD electrostatics
  83 program (Poisson-Boltzmann solver) can be made. A further prerequisite
  84 is that the input file is a run input file.
  85 The B-factor field is then filled with the Van der Waals radius
  86 of the atoms while the occupancy field will hold the charge.
  87 <p>
  88 The option <tt>-grasp</tt> is similar, but it puts the charges in the B-factor
  89 and the radius in the occupancy.
  90 <p>
  91 Option <tt>-align</tt> allows alignment
  92 of the principal axis of a specified group against the given vector,
  93 with an optional center of rotation specified by <tt>-aligncenter</tt>.
  94 <p>
  95 Finally, with option <tt>-label</tt>, <tt>editconf</tt> can add a chain identifier
  96 to a <tt>.<a href="pdb.html">pdb</a></tt> file, which can be useful for analysis with e.g. Rasmol.
  97 <p>
  98 To convert a truncated octrahedron file produced by a package which uses
  99 a cubic box with the corners cut off (such as GROMOS), use:<br>
 100 <tt>editconf -f in -rotate 0 45 35.264 -bt o -box veclen -o out</tt><br>
 101 where <tt>veclen</tt> is the size of the cubic box times sqrt(3)/2.
 102 <P>
 103 <H3>Files</H3>
 104 <TABLE BORDER=1 CELLSPACING=0 CELLPADDING=2>
 105 <TR><TH>option</TH><TH>filename</TH><TH>type</TH><TH>description</TH></TR>
 106 <TR><TD ALIGN=RIGHT> <b><tt>-f</tt></b> </TD><TD ALIGN=RIGHT> <tt><a href="files.html">    conf.gro</a></tt> </TD><TD> Input </TD><TD> Structure file: <a href="gro.html">gro</a> <a href="g96.html">g96</a> <a href="pdb.html">pdb</a> <a href="tpr.html">tpr</a> etc. </TD></TR>
 107 <TR><TD ALIGN=RIGHT> <b><tt>-n</tt></b> </TD><TD ALIGN=RIGHT> <tt><a href="ndx.html">   index.ndx</a></tt> </TD><TD> Input, Opt. </TD><TD> Index file </TD></TR>
 108 <TR><TD ALIGN=RIGHT> <b><tt>-o</tt></b> </TD><TD ALIGN=RIGHT> <tt><a href="files.html">     out.gro</a></tt> </TD><TD> Output, Opt. </TD><TD> Structure file: <a href="gro.html">gro</a> <a href="g96.html">g96</a> <a href="pdb.html">pdb</a> etc. </TD></TR>
 109 <TR><TD ALIGN=RIGHT> <b><tt>-mead</tt></b> </TD><TD ALIGN=RIGHT> <tt><a href="pqr.html">    mead.pqr</a></tt> </TD><TD> Output, Opt. </TD><TD> Coordinate file for MEAD </TD></TR>
 110 <TR><TD ALIGN=RIGHT> <b><tt>-bf</tt></b> </TD><TD ALIGN=RIGHT> <tt><a href="dat.html">   bfact.dat</a></tt> </TD><TD> Input, Opt. </TD><TD> Generic data file </TD></TR>
 111 </TABLE>
 112 <P>
 113 <H3>Other options</H3>
 114 <TABLE BORDER=1 CELLSPACING=0 CELLPADDING=2>
 115 <TR><TH>option</TH><TH>type</TH><TH>default</TH><TH>description</TH></TR>
 116 <TR><TD ALIGN=RIGHT> <b><tt>-[no]h</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Print help info and quit </TD></TD>
 117 <TR><TD ALIGN=RIGHT> <b><tt>-[no]version</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Print version info and quit </TD></TD>
 118 <TR><TD ALIGN=RIGHT> <b><tt>-nice</tt></b> </TD><TD ALIGN=RIGHT> int </TD><TD ALIGN=RIGHT> <tt>0</tt> </TD><TD> Set the nicelevel </TD></TD>
 119 <TR><TD ALIGN=RIGHT> <b><tt>-[no]w</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> View output <tt>.<a href="xvg.html">xvg</a></tt>, <tt>.<a href="xpm.html">xpm</a></tt>, <tt>.<a href="eps.html">eps</a></tt> and <tt>.<a href="pdb.html">pdb</a></tt> files </TD></TD>
 120 <TR><TD ALIGN=RIGHT> <b><tt>-[no]ndef</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Choose output from default index groups </TD></TD>
 121 <TR><TD ALIGN=RIGHT> <b><tt>-bt</tt></b> </TD><TD ALIGN=RIGHT> enum </TD><TD ALIGN=RIGHT> <tt>triclinic</tt> </TD><TD> Box type for <tt>-box</tt> and <tt>-d</tt>: <tt>triclinic</tt>, <tt>cubic</tt>, <tt>dodecahedron</tt> or <tt>octahedron</tt> </TD></TD>
 122 <TR><TD ALIGN=RIGHT> <b><tt>-box</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Box vector lengths (a,b,c) </TD></TD>
 123 <TR><TD ALIGN=RIGHT> <b><tt>-angles</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>90 90 90</tt> </TD><TD> Angles between the box vectors (bc,ac,ab) </TD></TD>
 124 <TR><TD ALIGN=RIGHT> <b><tt>-d</tt></b> </TD><TD ALIGN=RIGHT> real </TD><TD ALIGN=RIGHT> <tt>0     </tt> </TD><TD> Distance between the solute and the box </TD></TD>
 125 <TR><TD ALIGN=RIGHT> <b><tt>-[no]c</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Center molecule in box (implied by <tt>-box</tt> and <tt>-d</tt>) </TD></TD>
 126 <TR><TD ALIGN=RIGHT> <b><tt>-center</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Coordinates of geometrical center </TD></TD>
 127 <TR><TD ALIGN=RIGHT> <b><tt>-aligncenter</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Center of rotation for alignment </TD></TD>
 128 <TR><TD ALIGN=RIGHT> <b><tt>-align</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Align to target vector </TD></TD>
 129 <TR><TD ALIGN=RIGHT> <b><tt>-translate</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Translation </TD></TD>
 130 <TR><TD ALIGN=RIGHT> <b><tt>-rotate</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>0 0 0</tt> </TD><TD> Rotation around the X, Y and Z axes in degrees </TD></TD>
 131 <TR><TD ALIGN=RIGHT> <b><tt>-[no]princ</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Orient molecule(s) along their principal axes </TD></TD>
 132 <TR><TD ALIGN=RIGHT> <b><tt>-scale</tt></b> </TD><TD ALIGN=RIGHT> vector </TD><TD ALIGN=RIGHT> <tt>1 1 1</tt> </TD><TD> Scaling factor </TD></TD>
 133 <TR><TD ALIGN=RIGHT> <b><tt>-density</tt></b> </TD><TD ALIGN=RIGHT> real </TD><TD ALIGN=RIGHT> <tt>1000  </tt> </TD><TD> Density (g/L) of the output box achieved by scaling </TD></TD>
 134 <TR><TD ALIGN=RIGHT> <b><tt>-[no]pbc</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Remove the periodicity (make molecule whole again) </TD></TD>
 135 <TR><TD ALIGN=RIGHT> <b><tt>-resnr</tt></b> </TD><TD ALIGN=RIGHT> int </TD><TD ALIGN=RIGHT> <tt>-1</tt> </TD><TD>  Renumber residues starting from resnr </TD></TD>
 136 <TR><TD ALIGN=RIGHT> <b><tt>-[no]grasp</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Store the charge of the atom in the B-factor field and the radius of the atom in the occupancy field </TD></TD>
 137 <TR><TD ALIGN=RIGHT> <b><tt>-rvdw</tt></b> </TD><TD ALIGN=RIGHT> real </TD><TD ALIGN=RIGHT> <tt>0.12  </tt> </TD><TD> Default Van der Waals radius (in nm) if one can not be found in the database or if no parameters are present in the topology file </TD></TD>
 138 <TR><TD ALIGN=RIGHT> <b><tt>-[no]sig56</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Use rmin/2 (minimum in the Van der Waals potential) rather than &sigma;/2  </TD></TD>
 139 <TR><TD ALIGN=RIGHT> <b><tt>-[no]vdwread</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Read the Van der Waals radii from the file <tt>vdwradii.<a href="dat.html">dat</a></tt> rather than computing the radii based on the force field </TD></TD>
 140 <TR><TD ALIGN=RIGHT> <b><tt>-[no]atom</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Force B-factor attachment per atom </TD></TD>
 141 <TR><TD ALIGN=RIGHT> <b><tt>-[no]legend</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Make B-factor legend </TD></TD>
 142 <TR><TD ALIGN=RIGHT> <b><tt>-label</tt></b> </TD><TD ALIGN=RIGHT> string </TD><TD ALIGN=RIGHT> <tt>A</tt> </TD><TD> Add chain label for all residues </TD></TD>
 143 <TR><TD ALIGN=RIGHT> <b><tt>-[no]conect</tt></b> </TD><TD ALIGN=RIGHT> bool </TD><TD ALIGN=RIGHT> <tt>no    </tt> </TD><TD> Add CONECT records to a <tt>.<a href="pdb.html">pdb</a></tt> file when written. Can only be done when a topology is present </TD></TD>
 144 </TABLE>
 145 <P>
 146 <H3>Known problems</H3>
 147 <UL>
 148 <LI>For complex molecules, the periodicity removal routine may break down, in that case you can use <tt><a href="trjconv.html">trjconv</a></tt>.
 149 </UL>
 150 <P>
 151 <hr>
 152 <div ALIGN=RIGHT>
 153 <font size="-1"><a href="http://www.gromacs.org">http://www.gromacs.org</a></font><br>
 154 <font size="-1"><a href="mailto:gromacs@gromacs.org">gromacs@gromacs.org</a></font><br>
 155 </div>
 156 </BODY>