next up previous contents index
Next: Generalized Born Implicit Solvent Up: Force Field Parameters Previous: MARTINI Residue-Based Coarse-Grain Forcefield   Contents   Index

Subsections

Constraints and Restraints


Bond constraint parameters

Position restraint parameters

The following describes the parameters for the position restraints feature of NAMD. For historical reasons the term ``constraints'' has been carried over from X-PLOR. This feature allows a restraining potential to each atom of an arbitrary set during the simulation.

Fixed atoms parameters

Atoms may be held fixed during a simulation. NAMD avoids calculating most interactions in which all affected atoms are fixed unless fixedAtomsForces is specified.

Extra bond, angle, and dihedral restraints

Additional bond, angle, and dihedral energy terms may be applied to system, allowing secondary or tertiary structure to be restrained, for example. Extra bonded terms are not considered part of the molecular structure and hence do not alter nonbonded exclusions. The energies from extra bonded terms are included with the normal bond, angle, and dihedral energies in NAMD output.

All extra bonded terms are harmonic potentials of the form $ U(x) = k (x-x_{ref})^2$ except dihedrals and impropers with a non-zero periodicity specified, which use $ U(x) = k (1 + cos(n x - x_{ref}))$ . The only difference between dihedrals and impropers is the output field that their potential energy is added to.

Due to a very old bug all NAMD releases prior to 2.13 have used the MARTINI cosine-based angle potential function for all extra angles. Since workflows may unknowingly depend on this undocumented behavior, cosine-based angles remain the default, but a warning is printed unless the desired behavior is specified via the new option extraBondsCosAngles (defaults to ``on'', set to ``off'' to use the normal harmonic angle potential function for all extra angles).

The extra bonded term implementation shares the parallel implementation of regular bonded terms in NAMD, allowing large numbers of extra terms to be specified with minimal impact on parallel scalability. Extra bonded terms do not have to duplicate normal bonds/angles/dihedrals, but each extra bond/angle/dihedral should only involve nearby atoms. If the atoms involved are too far apart a bad global bond count will be reported in parallel runs.

Extra bonded terms are enabled via the following options:

The extra bonds file(s) should contain lines of the following formats:

In all cases <atom> is a zero-based atom index (the first atom has index 0), <ref> is a reference distance in Å (bond) or angle in degrees (others), and <k> is a spring constant in the potential energy function $ U(x) = k (x-x_{ref})^2$ or, for dihedrals and impropers with periodicity <n> specified and not 0, $ U(x) = k (1 + cos(n x - x_{ref}))$ . Note that $ x_{ref}$ is only a minimum for the harmonic potential; the sinusoidal potential has minima at $ (x_{ref} + 180)/n + i \times 360/n$ .

Use of wall implements a harmonic wall potential similar to the Colvars harmonic wall restraint. The potential function is

\begin{displaymath}
U(x) =
\begin{cases}
k(x-x_{\text{upper}})^2, &\text{if $x ...
...ext{lower}})^2, &\text{if $x < x_{\text{lower}}$}
\end{cases}.
\end{displaymath}


next up previous contents index
Next: Generalized Born Implicit Solvent Up: Force Field Parameters Previous: MARTINI Residue-Based Coarse-Grain Forcefield   Contents   Index
http://www.ks.uiuc.edu/Research/namd/