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

Harmonic restraint parameters

The following describes the parameters for the harmonic restraints feature of NAMD. For historical reasons the terminology of ``harmonic constraints'' has been carried over from X-PLOR. This feature allows a harmonic restraining force to be applied to any set of atoms in 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.

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$.


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/