2 |
|
\section{\label{sec:DUFF}Dipolar Unified-Atom Force Field} |
3 |
|
|
4 |
|
The \underline{D}ipolar \underline{U}nified-Atom |
5 |
< |
\underline{F}orce \underline{F}ield (DUFF) was developed to |
5 |
> |
\underline{F}orce \underline{F}ield ({\sc duff}) was developed to |
6 |
|
simulate lipid bilayers. We needed a model capable of forming |
7 |
|
bilayers, while still being sufficiently computationally efficient to |
8 |
|
allow simulations of large systems ($\approx$100's of phospholipids, |
17 |
|
computationally expensive Ewald-Sum. Instead, we can use |
18 |
|
neighbor-lists and cutoff radii for the dipolar interactions. |
19 |
|
|
20 |
< |
As an example, lipid head groups in DUFF are represented as point |
20 |
> |
As an example, lipid head groups in {\sc duff} are represented as point |
21 |
|
dipole interaction sites.PC and PE Lipid head groups are typically |
22 |
|
zwitterionic in nature, with charges separated by as much as |
23 |
|
6~$\mbox{\AA}$. By placing a dipole of 20.6~Debye at the head group |
27 |
|
atom in Fig.~\ref{fig:lipidModel}. |
28 |
|
|
29 |
|
\begin{figure} |
30 |
< |
\includegraphics[angle=-90,width=\linewidth]{lipidModel.epsi} |
30 |
> |
\epsfxsize=6in |
31 |
> |
\epsfbox{lipidModel.epsi} |
32 |
|
\caption{A representation of the lipid model. $\phi$ is the torsion angle, $\theta$ % |
33 |
|
is the bend angle, $\mu$ is the dipole moment of the head group, and n is the chain length.} |
34 |
|
\label{fig:lipidModel} |
67 |
|
used when integrating the equations of motion. |
68 |
|
|
69 |
|
|
70 |
< |
\subsection{\label{subSec:energyFunctions}DUFF Energy Functions} |
70 |
> |
\subsection{\label{subSec:energyFunctions}{\sc duff} Energy Functions} |
71 |
|
|
72 |
< |
The total energy of function in DUFF is given by the following: |
72 |
> |
The total energy of function in {\sc duff} is given by the following: |
73 |
|
\begin{equation} |
74 |
|
V_{\text{Total}} = \sum^{N}_{I=1} V^{I}_{\text{Internal}} |
75 |
|
+ \sum^{N}_{I=1} \sum^{N}_{J=1} V^{IJ}_{\text{Cross}} |