51 |
|
reproduce the behavior exhibited in simulations using SPME with an |
52 |
|
$\mathscr{O}(N)$ computational cost, equivalent to merely the |
53 |
|
real-space portion of the lattice summation. |
54 |
+ |
|
55 |
|
\end{abstract} |
56 |
|
|
57 |
|
\newpage |
267 |
|
|
268 |
|
\begin{figure} |
269 |
|
\centering |
270 |
< |
\includegraphics[width = \linewidth]{./linearFit.pdf} |
271 |
< |
\caption{Example least squares regression of the $\Delta E$ between configurations for the SF method against SPME in the pure water system. } |
270 |
> |
\includegraphics[width = \linewidth]{./dualLinear.pdf} |
271 |
> |
\caption{Example least squares regressions of the configuration energy differences for SPC/E water systems. The upper plot shows a data set with a poor correlation coefficient ($R^2$), while the lower plot shows a data set with a good correlation coefficient.} |
272 |
|
\label{fig:linearFit} |
273 |
|
\end{figure} |
274 |
|
|
275 |
< |
Each system type (detailed in section \ref{sec:RepSims}) studied consisted of 500 independent configurations, each equilibrated from higher temperature trajectories. Thus, 124,750 $\Delta E$ data points are used in a regression of a single system type. Results and discussion for the individual analysis of each of the system types appear in the supporting information, while the cumulative results over all the investigated systems appears below in section \ref{sec:EnergyResults}. |
275 |
> |
Each system type (detailed in section \ref{sec:RepSims}) studied |
276 |
> |
consisted of 500 independent configurations, each equilibrated from |
277 |
> |
higher temperature trajectories. Thus, 124,750 $\Delta E$ data points |
278 |
> |
are used in a regression of a single system type. Results and |
279 |
> |
discussion for the individual analysis of each of the system types |
280 |
> |
appear in the supporting information, while the cumulative results |
281 |
> |
over all the investigated systems appears below in section |
282 |
> |
\ref{sec:EnergyResults}. |
283 |
|
|
284 |
|
\subsection{Molecular Dynamics and the Force and Torque Vectors}\label{sec:MDMethods} |
285 |
|
Evaluation of the pairwise methods (outlined in section \ref{sec:ESMethods}) for use in MD simulations was performed through comparison of the force and torque vectors obtained with those from SPME. Both the magnitude and the direction of these vectors on each of the bodies in the system were analyzed. For the magnitude of these vectors, linear least squares regression analysis can be performed as described previously for comparing $\Delta E$ values. Instead of a single value between two system configurations, there is a value for each particle in each configuration. For a system of 1000 water molecules and 40 ions, there are 1040 force vectors and 1000 torque vectors. With 500 configurations, this results in 520,000 force and 500,000 torque vector comparisons samples for each system type. |