next up previous contents index
Next: Frozen Atoms Up: The Intramolecular Potential Functions Previous: Inversion Angle Potentials   Contents   Index

Tethering Forces

DL_POLY_2 also allows atomic sites to be tethered to a fixed point in space, $\mbox{$\underline{r}$}_0$ taken as their position at the beginning of the simulation. This is also known as position restraining. The specification, which comes as part of the molecular description, requires a tether potential type and the associated interaction parameters.

Note, firstly, that application of tethering potentials means that momentum will no longer be a conserved quantity of the simulation. Secondly, in constant pressure simulations, where the MD cell changes size or shape, the reference position is scaled with the cell vectors.

The potential functions available in DL_POLY_2 are as follows, in each case $r_{i0}$ is the distance of the atom from its position at $t=0$:

  1. harmonic potential: (harm)
    \begin{displaymath}
U(r_{i0}) = \frac{1}{2}k(r_{i0})^2;
\end{displaymath} (2.63)

  2. restrained harmonic :(rhrm)
    $\displaystyle U(r_{i0})$ $\textstyle =$ $\displaystyle \frac{1}{2}k(r_{i0})^2~~~~~~r_{i0}\le r_{c};$ (2.64)
    $\displaystyle U(r_{i0})$ $\textstyle =$ $\displaystyle \frac{1}{2}kr_{c}^2+kr_{c}(r_{i0}-r_{c})~~~~~~r_{i0}>r_{c};$ (2.65)

  3. Quartic potential: (quar)
    \begin{displaymath}
U(r_{i0})=\frac{k}{2}(r_{i0})^2+\frac{k'}{3}(r_{i0})^3+\frac{k''}{4}(r_{i0})^4.
\end{displaymath} (2.66)

The force on the atom $i$ arising from a tether potential is obtained using the general formula:

\begin{displaymath}
\mbox{$\underline{f}$}_{i}=-\frac{1}{r_{i0}}\left[
\frac{\pa...
...}{\partial r_{i0}}U(r_{i0})\right]\mbox{$\underline{r}$}_{i0},
\end{displaymath} (2.67)

The contribution to be added to the atomic virial is given by

\begin{displaymath}
{\cal W}=\mbox{$\underline{r}$}_{i0}\cdot \mbox{$\underline{f}$}_{i},
\end{displaymath} (2.68)

The contribution to be added to the atomic stress tensor is given by

\begin{displaymath}
\sigma^{\alpha \beta}=-r_{i0}^{\alpha}f_{i}^{\beta},
\end{displaymath} (2.69)

where $\alpha$ and $\beta$ indicate the $x,y,z$ components. The atomic stress tensor derived in this way is symmetric.

In DL_POLY_2 bond forces are handled by the routine TETHFRC.


next up previous contents index
Next: Frozen Atoms Up: The Intramolecular Potential Functions Previous: Inversion Angle Potentials   Contents   Index
W Smith 2003-05-12