--- trunk/electrostaticMethodsPaper/electrostaticMethods.tex	2006/03/19 19:34:53	2638
+++ trunk/electrostaticMethodsPaper/electrostaticMethods.tex	2006/03/20 02:00:26	2639
@@ -77,25 +77,56 @@ accurately incorporate their effect, and since the com
 leading to an effect excluded from the pair interactions within a unit
 box.  In large systems, excessively large cutoffs need to be used to
 accurately incorporate their effect, and since the computational cost
-increases proportionally with the cutoff sphere, it quickly becomes an
-impractical task to perform these calculations.
+increases proportionally with the cutoff sphere, it quickly becomes
+very time-consuming to perform these calculations.
 
+There have been many efforts to address this issue of both proper and
+practical handling of electrostatic interactions, and these have
+resulted in the availability of a variety of
+techniques.\cite{Roux99,Sagui99,Tobias01} These are typically
+classified as implicit methods (i.e., continuum dielectrics, static
+dipolar fields),\cite{Born20,Grossfield00} explicit methods (i.e.,
+Ewald summations, interaction shifting or
+trucation),\cite{Ewald21,Steinbach94} or a mixture of the two (i.e.,
+reaction field type methods, fast multipole
+methods).\cite{Onsager36,Rokhlin85} The explicit or mixed methods are
+often preferred because they incorporate dynamic solvent molecules in
+the system of interest, but these methods are sometimes difficult to
+utilize because of their high computational cost.\cite{Roux99} In
+addition to this cost, there has been some question of the inherent
+periodicity of the explicit Ewald summation artificially influencing
+systems dynamics.\cite{Tobias01}
+
+In this paper, we focus on the common mixed and explicit methods of
+reaction filed and smooth particle mesh
+Ewald\cite{Onsager36,Essmann99} and a new set of shifted methods
+devised by Wolf {\it et al.} which we further extend.\cite{Wolf99}
+These new methods for handling electrostatics are quite
+computationally efficient, since they involve only a simple
+modification to the direct pairwise sum, and they lack the added
+periodicity of the Ewald sum. Below, these methods are evaluated using
+a variety of model systems and comparison methodologies to establish
+their useability in molecular simulations.
+
 \subsection{The Ewald Sum}
-The complete accumulation electrostatic interactions in a system with periodic boundary conditions (PBC) requires the consideration of the effect of all charges within a simulation box, as well as those in the periodic replicas,
+The complete accumulation electrostatic interactions in a system with
+periodic boundary conditions (PBC) requires the consideration of the
+effect of all charges within a simulation box, as well as those in the
+periodic replicas,
 \begin{equation}
 V_\textrm{elec} = \frac{1}{2} {\sum_{\mathbf{n}}}^\prime \left[ \sum_{i=1}^N\sum_{j=1}^N \phi\left( \mathbf{r}_{ij} + L\mathbf{n},\bm{\Omega}_i,\bm{\Omega}_j\right) \right],
 \label{eq:PBCSum}
 \end{equation}
-where the sum over $\mathbf{n}$ is a sum over all periodic box replicas
-with integer coordinates $\mathbf{n} = (l,m,n)$, and the prime indicates
-$i = j$ are neglected for $\mathbf{n} = 0$.\cite{deLeeuw80} Within the
-sum, $N$ is the number of electrostatic particles, $\mathbf{r}_{ij}$ is
-$\mathbf{r}_j - \mathbf{r}_i$, $L$ is the cell length, $\bm{\Omega}_{i,j}$ are
-the Euler angles for $i$ and $j$, and $\phi$ is Poisson's equation
-($\phi(\mathbf{r}_{ij}) = q_i q_j |\mathbf{r}_{ij}|^{-1}$ for charge-charge
-interactions). In the case of monopole electrostatics,
-eq. (\ref{eq:PBCSum}) is conditionally convergent and is discontiuous
-for non-neutral systems.
+where the sum over $\mathbf{n}$ is a sum over all periodic box
+replicas with integer coordinates $\mathbf{n} = (l,m,n)$, and the
+prime indicates $i = j$ are neglected for $\mathbf{n} =
+0$.\cite{deLeeuw80} Within the sum, $N$ is the number of electrostatic
+particles, $\mathbf{r}_{ij}$ is $\mathbf{r}_j - \mathbf{r}_i$, $L$ is
+the cell length, $\bm{\Omega}_{i,j}$ are the Euler angles for $i$ and
+$j$, and $\phi$ is Poisson's equation ($\phi(\mathbf{r}_{ij}) = q_i
+q_j |\mathbf{r}_{ij}|^{-1}$ for charge-charge interactions). In the
+case of monopole electrostatics, eq. (\ref{eq:PBCSum}) is
+conditionally convergent and is discontiuous for non-neutral systems.
 
 This electrostatic summation problem was originally studied by Ewald
 for the case of an infinite crystal.\cite{Ewald21}. The approach he
@@ -542,7 +573,7 @@ between those computed from the particular method and 
 investigated through measurement of the angle ($\theta$) formed
 between those computed from the particular method and those from SPME,
 \begin{equation}
-\theta_f = \cos^{-1} \hat{f}_\textrm{SPME} \cdot \hat{f}_\textrm{Method},
+\theta_f = \cos^{-1} \left(\hat{f}_\textrm{SPME} \cdot \hat{f}_\textrm{Method}\right),
 \end{equation}
 where $\hat{f}_\textrm{M}$ is the unit vector pointing along the
 force vector computed using method $M$.