Added information about dihedral restraints.
[alexxy/gromacs.git] / manual / topology.tex
index 35e4e821cede1676c9462f3cf1ce5451b3c000af..87782be2f792f48655dd0239f6203477670c77a3 100644 (file)
@@ -1365,6 +1365,11 @@ Urea         3
    2     1     1000       0    1000 ; Restrain to a line (Y-axis)
    3     1     1000       0       0 ; Restrain to a plane (Y-Z-plane)
 
+[ dihedral_restraints ]
+; ai   aj    ak    al  type  label  phi  dphi  kfac  power
+    3    6     1    2     1      1  180     0     1      2
+    1    4     3    5     1      1  180     0     1      2
+
 ; Include TIP3P water topology
 #include "amber99/tip3p.itp"
 
@@ -1474,6 +1479,19 @@ to reference positions (\ssecref{positionrestraint}).
 The reference positions are read from a 
 separate coordinate file by {\tt \normindex{grompp}}.
 
+
+{\bf {\tt [~dihedral_restraints~]} :} restrain selected dihedrals to a reference value.
+The implementation of dihedral restraints is described in section \ssecref{dihedralrestraint} of the manual.
+The parameters specified in the [dihedral_restraints] directive are as follows:
+\begin{itemize}
+\item {\tt type} has only one possible value which is 1
+\item {\tt label} is unused and has been removed from the code.
+\item {\tt phi} is the value of $\phi_0$ in \eqnref{dphi} and \eqnref{dihre} of the manual.
+\item {\tt dphi} is the value of $\Delta\phi$ in \eqnref{dihre} of the manual.
+\item {\tt kfac} is analogous to {\tt fac} in the implementation of distance restraints. It is the factor by which the force constant is multiplied. By doing so, different restraints can be maintained with different force constants.
+\item {\tt power} is unused and has been removed from the code.
+\end{itemize}
+
 {\bf {\tt \#include "tip3p.itp"} :} includes a topology file that was already
 constructed (see section~\ssecref{molitp}).