| 575 |
|
\theta_f = \cos^{-1} \left(\hat{F}_\textrm{SPME} \cdot \hat{F}_\textrm{M}\right), |
| 576 |
|
\end{equation} |
| 577 |
|
where $\hat{f}_\textrm{M}$ is the unit vector pointing along the force |
| 578 |
< |
vector computed using method M. |
| 579 |
< |
|
| 580 |
< |
Each of these $\theta$ values was accumulated in a distribution |
| 581 |
< |
function and weighted by the area on the unit sphere. Non-linear |
| 582 |
< |
Gaussian fits were used to measure the width of the resulting |
| 583 |
< |
distributions. |
| 584 |
< |
|
| 585 |
< |
\begin{figure} |
| 586 |
< |
\centering |
| 587 |
< |
\includegraphics[width = \linewidth]{./gaussFit.pdf} |
| 588 |
< |
\caption{Sample fit of the angular distribution of the force vectors |
| 589 |
< |
accumulated using all of the studied systems. Gaussian fits were used |
| 590 |
< |
to obtain values for the variance in force and torque vectors.} |
| 591 |
< |
\label{fig:gaussian} |
| 592 |
< |
\end{figure} |
| 593 |
< |
|
| 594 |
< |
Figure \ref{fig:gaussian} shows an example distribution with applied |
| 595 |
< |
non-linear fits. The solid line is a Gaussian profile, while the |
| 596 |
< |
dotted line is a Voigt profile, a convolution of a Gaussian and a |
| 597 |
< |
Lorentzian. Since this distribution is a measure of angular error |
| 578 |
> |
vector computed using method M. Each of these $\theta$ values was |
| 579 |
> |
accumulated in a distribution function and weighted by the area on the |
| 580 |
> |
unit sphere. Since this distribution is a measure of angular error |
| 581 |
|
between two different electrostatic summation methods, there is no |
| 582 |
< |
{\it a priori} reason for the profile to adhere to any specific shape. |
| 583 |
< |
Gaussian fits was used to compare all the tested methods. The |
| 584 |
< |
variance ($\sigma^2$) was extracted from each of these fits and was |
| 585 |
< |
used to compare distribution widths. Values of $\sigma^2$ near zero |
| 586 |
< |
indicate vector directions indistinguishable from those calculated |
| 587 |
< |
when using the reference method (SPME). |
| 582 |
> |
{\it a priori} reason for the profile to adhere to any specific |
| 583 |
> |
shape. Thus, gaussian fits were used to measure the width of the |
| 584 |
> |
resulting distributions. |
| 585 |
> |
% |
| 586 |
> |
%\begin{figure} |
| 587 |
> |
%\centering |
| 588 |
> |
%\includegraphics[width = \linewidth]{./gaussFit.pdf} |
| 589 |
> |
%\caption{Sample fit of the angular distribution of the force vectors |
| 590 |
> |
%accumulated using all of the studied systems. Gaussian fits were used |
| 591 |
> |
%to obtain values for the variance in force and torque vectors.} |
| 592 |
> |
%\label{fig:gaussian} |
| 593 |
> |
%\end{figure} |
| 594 |
> |
% |
| 595 |
> |
%Figure \ref{fig:gaussian} shows an example distribution with applied |
| 596 |
> |
%non-linear fits. The solid line is a Gaussian profile, while the |
| 597 |
> |
%dotted line is a Voigt profile, a convolution of a Gaussian and a |
| 598 |
> |
%Lorentzian. |
| 599 |
> |
%Since this distribution is a measure of angular error between two |
| 600 |
> |
%different electrostatic summation methods, there is no {\it a priori} |
| 601 |
> |
%reason for the profile to adhere to any specific shape. |
| 602 |
> |
%Gaussian fits was used to compare all the tested methods. |
| 603 |
> |
The variance ($\sigma^2$) was extracted from each of these fits and |
| 604 |
> |
was used to compare distribution widths. Values of $\sigma^2$ near |
| 605 |
> |
zero indicate vector directions indistinguishable from those |
| 606 |
> |
calculated when using the reference method (SPME). |
| 607 |
|
|
| 608 |
|
\subsection{Short-time Dynamics} |
| 609 |
|
|
| 678 |
|
the resulting configurations were again equilibrated individually. |
| 679 |
|
Finally, for the Argon / Water ``charge void'' systems, the identities |
| 680 |
|
of all the SPC/E waters within 6 \AA\ of the center of the |
| 681 |
< |
equilibrated water configurations were converted to argon |
| 682 |
< |
(Fig. \ref{fig:argonSlice}). |
| 681 |
> |
equilibrated water configurations were converted to argon. |
| 682 |
> |
%(Fig. \ref{fig:argonSlice}). |
| 683 |
|
|
| 684 |
|
These procedures guaranteed us a set of representative configurations |
| 685 |
|
from chemically-relevant systems sampled from an appropriate |
| 686 |
|
ensemble. Force field parameters for the ions and Argon were taken |
| 687 |
|
from the force field utilized by {\sc oopse}.\cite{Meineke05} |
| 688 |
|
|
| 689 |
< |
\begin{figure} |
| 690 |
< |
\centering |
| 691 |
< |
\includegraphics[width = \linewidth]{./slice.pdf} |
| 692 |
< |
\caption{A slice from the center of a water box used in a charge void |
| 693 |
< |
simulation. The darkened region represents the boundary sphere within |
| 694 |
< |
which the water molecules were converted to argon atoms.} |
| 695 |
< |
\label{fig:argonSlice} |
| 696 |
< |
\end{figure} |
| 689 |
> |
%\begin{figure} |
| 690 |
> |
%\centering |
| 691 |
> |
%\includegraphics[width = \linewidth]{./slice.pdf} |
| 692 |
> |
%\caption{A slice from the center of a water box used in a charge void |
| 693 |
> |
%simulation. The darkened region represents the boundary sphere within |
| 694 |
> |
%which the water molecules were converted to argon atoms.} |
| 695 |
> |
%\label{fig:argonSlice} |
| 696 |
> |
%\end{figure} |
| 697 |
|
|
| 698 |
|
\subsection{Comparison of Summation Methods}\label{sec:ESMethods} |
| 699 |
|
We compared the following alternative summation methods with results |
| 803 |
|
\begin{figure} |
| 804 |
|
\centering |
| 805 |
|
\includegraphics[width=5.5in]{./frcMagplot.pdf} |
| 806 |
< |
\caption{Statistical analysis of the quality of the force vector magnitudes for a given electrostatic method compared with the reference Ewald sum. Results with a value equal to 1 (dashed line) indicate force magnitude values indistinguishable from those obtained using SPME. Different values of the cutoff radius are indicated with different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = inverted triangles).} |
| 806 |
> |
\caption{Statistical analysis of the quality of the force vector |
| 807 |
> |
magnitudes for a given electrostatic method compared with the |
| 808 |
> |
reference Ewald sum. Results with a value equal to 1 (dashed line) |
| 809 |
> |
indicate force magnitude values indistinguishable from those obtained |
| 810 |
> |
using SPME. Different values of the cutoff radius are indicated with |
| 811 |
> |
different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = |
| 812 |
> |
inverted triangles).} |
| 813 |
|
\label{fig:frcMag} |
| 814 |
|
\end{figure} |
| 815 |
|
|
| 835 |
|
\begin{figure} |
| 836 |
|
\centering |
| 837 |
|
\includegraphics[width=5.5in]{./trqMagplot.pdf} |
| 838 |
< |
\caption{Statistical analysis of the quality of the torque vector magnitudes for a given electrostatic method compared with the reference Ewald sum. Results with a value equal to 1 (dashed line) indicate torque magnitude values indistinguishable from those obtained using SPME. Different values of the cutoff radius are indicated with different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = inverted triangles).} |
| 838 |
> |
\caption{Statistical analysis of the quality of the torque vector |
| 839 |
> |
magnitudes for a given electrostatic method compared with the |
| 840 |
> |
reference Ewald sum. Results with a value equal to 1 (dashed line) |
| 841 |
> |
indicate torque magnitude values indistinguishable from those obtained |
| 842 |
> |
using SPME. Different values of the cutoff radius are indicated with |
| 843 |
> |
different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = |
| 844 |
> |
inverted triangles).} |
| 845 |
|
\label{fig:trqMag} |
| 846 |
|
\end{figure} |
| 847 |
|
|
| 871 |
|
\begin{figure} |
| 872 |
|
\centering |
| 873 |
|
\includegraphics[width=5.5in]{./frcTrqAngplot.pdf} |
| 874 |
< |
\caption{Statistical analysis of the quality of the Gaussian fit of the force and torque vector angular distributions for a given electrostatic method compared with the reference Ewald sum. Results with a variance ($\sigma^2$) equal to zero (dashed line) indicate force and torque directions indistinguishable from those obtained using SPME. Different values of the cutoff radius are indicated with different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = inverted triangles).} |
| 874 |
> |
\caption{Statistical analysis of the quality of the Gaussian fit of |
| 875 |
> |
the force and torque vector angular distributions for a given |
| 876 |
> |
electrostatic method compared with the reference Ewald sum. Results |
| 877 |
> |
with a variance ($\sigma^2$) equal to zero (dashed line) indicate |
| 878 |
> |
force and torque directions indistinguishable from those obtained |
| 879 |
> |
using SPME. Different values of the cutoff radius are indicated with |
| 880 |
> |
different symbols (9\AA\ = circles, 12\AA\ = squares, and 15\AA\ = |
| 881 |
> |
inverted triangles).} |
| 882 |
|
\label{fig:frcTrqAng} |
| 883 |
|
\end{figure} |
| 884 |
|
|
| 902 |
|
|
| 903 |
|
\begin{table}[htbp] |
| 904 |
|
\centering |
| 905 |
< |
\caption{Variance ($\sigma^2$) of the force (top set) and torque (bottom set) vector angle difference distributions for the Shifted Potential and Shifted Force methods. Calculations were performed both with (Y) and without (N) group based cutoffs and a switching function. The $\alpha$ values have units of \AA$^{-1}$ and the variance values have units of degrees$^2$.} |
| 905 |
> |
\caption{Variance ($\sigma^2$) of the force (top set) and torque |
| 906 |
> |
(bottom set) vector angle difference distributions for the Shifted Potential and Shifted Force methods. Calculations were performed both with (Y) and without (N) group based cutoffs and a switching function. The $\alpha$ values have units of \AA$^{-1}$ and the variance values have units of degrees$^2$.} |
| 907 |
|
\begin{tabular}{@{} ccrrrrrrrr @{}} |
| 908 |
|
\\ |
| 909 |
|
\toprule |
| 991 |
|
\begin{figure} |
| 992 |
|
\centering |
| 993 |
|
\includegraphics[width = \linewidth]{./vCorrPlot.pdf} |
| 994 |
< |
\caption{Velocity auto-correlation functions of NaCl crystals at 1000 K while using SPME, {\sc sf} ($\alpha$ = 0.0, 0.1, \& 0.2), and {\sc sp} ($\alpha$ = 0.2). The inset is a magnification of the first trough. The times to first collision are nearly identical, but the differences can be seen in the peaks and troughs, where the undamped to weakly damped methods are stiffer than the moderately damped and SPME methods.} |
| 994 |
> |
\caption{Velocity auto-correlation functions of NaCl crystals at |
| 995 |
> |
1000 K while using SPME, {\sc sf} ($\alpha$ = 0.0, 0.1, \& 0.2), and |
| 996 |
> |
{\sc sp} ($\alpha$ = 0.2). The inset is a magnification of the first |
| 997 |
> |
trough. The times to first collision are nearly identical, but the |
| 998 |
> |
differences can be seen in the peaks and troughs, where the undamped |
| 999 |
> |
to weakly damped methods are stiffer than the moderately damped and |
| 1000 |
> |
SPME methods.} |
| 1001 |
|
\label{fig:vCorrPlot} |
| 1002 |
|
\end{figure} |
| 1003 |
|
|
| 1032 |
|
\begin{figure} |
| 1033 |
|
\centering |
| 1034 |
|
\includegraphics[width = \linewidth]{./spectraSquare.pdf} |
| 1035 |
< |
\caption{Power spectra obtained from the velocity auto-correlation functions of NaCl crystals at 1000 K while using SPME, {\sc sf} ($\alpha$ = 0, 0.1, \& 0.2), and {\sc sp} ($\alpha$ = 0.2). Apodization of the correlation functions via a cubic switching function between 40 and 50 ps was used to clear up the spectral noise resulting from data truncation, and had no noticeable effect on peak location or magnitude. The inset shows the frequency region below 100 cm$^{-1}$ to highlight where the spectra begin to differ.} |
| 1035 |
> |
\caption{Power spectra obtained from the velocity auto-correlation |
| 1036 |
> |
functions of NaCl crystals at 1000 K while using SPME, {\sc sf} |
| 1037 |
> |
($\alpha$ = 0, 0.1, \& 0.2), and {\sc sp} ($\alpha$ = 0.2). |
| 1038 |
> |
Apodization of the correlation functions via a cubic switching |
| 1039 |
> |
function between 40 and 50 ps was used to clear up the spectral noise |
| 1040 |
> |
resulting from data truncation, and had no noticeable effect on peak |
| 1041 |
> |
location or magnitude. The inset shows the frequency region below 100 |
| 1042 |
> |
cm$^{-1}$ to highlight where the spectra begin to differ.} |
| 1043 |
|
\label{fig:methodPS} |
| 1044 |
|
\end{figure} |
| 1045 |
|
|
| 1081 |
|
\begin{figure} |
| 1082 |
|
\centering |
| 1083 |
|
\includegraphics[width = \linewidth]{./comboSquare.pdf} |
| 1084 |
< |
\caption{Regions of spectra showing the low-frequency correlated motions for NaCl crystals at 1000 K using various electrostatic summation methods. The upper plot is a zoomed inset from figure \ref{fig:methodPS}. As the damping value for the {\sc sf} potential increases, the low-frequency peaks red-shift. The lower plot is of spectra when using SPME and a simple damped Coulombic sum with damping coefficients ($\alpha$) ranging from 0.15 to 0.3 \AA$^{-1}$. As $\alpha$ increases, the peaks are red-shifted toward and eventually beyond the values given by SPME. The larger $\alpha$ values weaken the real-space electrostatics, explaining this shift towards less strongly correlated motions in the crystal.} |
| 1084 |
> |
\caption{Regions of spectra showing the low-frequency correlated |
| 1085 |
> |
motions for NaCl crystals at 1000 K using various electrostatic |
| 1086 |
> |
summation methods. The upper plot is a zoomed inset from figure |
| 1087 |
> |
\ref{fig:methodPS}. As the damping value for the {\sc sf} potential |
| 1088 |
> |
increases, the low-frequency peaks red-shift. The lower plot is of |
| 1089 |
> |
spectra when using SPME and a simple damped Coulombic sum with damping |
| 1090 |
> |
coefficients ($\alpha$) ranging from 0.15 to 0.3 \AA$^{-1}$. As |
| 1091 |
> |
$\alpha$ increases, the peaks are red-shifted toward and eventually |
| 1092 |
> |
beyond the values given by SPME. The larger $\alpha$ values weaken |
| 1093 |
> |
the real-space electrostatics, explaining this shift towards less |
| 1094 |
> |
strongly correlated motions in the crystal.} |
| 1095 |
|
\label{fig:dampInc} |
| 1096 |
|
\end{figure} |
| 1097 |
|
|
