next up previous contents
Next: Pair list generation Up: Non-bonded interaction parameters and Previous: Non-bonded van der Waals

Non-bonded electrostatic interactions

The handling of electrostatics is slightly more complicated due to the incorporation of DPMTA to provide full electrostatic interactions. There are two cases to consider, one where DPMTA is used to provide full electrostatics and the other where electrostatics are truncated at a given distance.

First let us consider the latter case, where electrostatics are truncated at the cutoff distance. Using this scheme, all electrostatic interactions beyond a specified distance are ignored, that is, assumed to be zero. If switching is set to on, rather than having a discontinuity in the potential at the cutoff distance, a shifting function is applied to the electrostatic potential as shown in Figure 2. As this figure shows, the shifting function shifts the entire potential curve so that the curve intersects the x-axis at the cutoff distance. This shifting function is based on the shifting function used by X-PLOR. See section 5.1.4 for more details.


 \begin{figure}% latex2html id marker 3753\begin{center}
\epsfig{file=figures/s...
...tential with and without the
application of the shifting function.}}\end{figure}

Next, consider the case where full electrostatics are calculated. In this case, the electrostatic interactions are not truncated at any distance. In this scheme, the cutoff parameter has a slightly different meaning for the electrostatic interactions -- it represents the local interaction distance, or distance within which electrostatic pairs will be directly calculated every timestep. Outside of this distance, interactions will be calculated only periodically using DPMTA. These forces will be applied using a multiple timestep integration scheme as described in section 4.5.4.


 \begin{figure}% latex2html id marker 3762\begin{center}
\epsfig{file=figures/f...
...portion calculated directly
and the other calculated using DPMTA.}}\end{figure}


next up previous contents
Next: Pair list generation Up: Non-bonded interaction parameters and Previous: Non-bonded van der Waals
David Hardy
1998-09-06