| 572 |
|
when using the reference method (SPME). |
| 573 |
|
|
| 574 |
|
\subsection{Short-time Dynamics} |
| 575 |
< |
|
| 576 |
< |
\subsection{Long-Time and Collective Motion}\label{sec:LongTimeMethods} |
| 577 |
< |
Evaluation of the long-time dynamics of charged systems was performed |
| 578 |
< |
by considering the NaCl crystal system while using a subset of the |
| 575 |
> |
Evaluation of the short-time dynamics of charged systems was performed |
| 576 |
> |
by considering the 1000 K NaCl crystal system while using a subset of the |
| 577 |
|
best performing pairwise methods. The NaCl crystal was chosen to |
| 578 |
|
avoid possible complications involving the propagation techniques of |
| 579 |
< |
orientational motion in molecular systems. To enhance the atomic |
| 580 |
< |
motion, these crystals were equilibrated at 1000 K, near the |
| 581 |
< |
experimental $T_m$ for NaCl. Simulations were performed under the |
| 582 |
< |
microcanonical ensemble, and velocity autocorrelation functions |
| 583 |
< |
(Eq. \ref{eq:vCorr}) were computed for each of the trajectories, |
| 579 |
> |
orientational motion in molecular systems. All systems were started |
| 580 |
> |
with the same initial positions and velocities. Simulations were |
| 581 |
> |
performed under the microcanonical ensemble, and velocity |
| 582 |
> |
autocorrelation functions (Eq. \ref{eq:vCorr}) were computed for each |
| 583 |
> |
of the trajectories, |
| 584 |
|
\begin{equation} |
| 585 |
< |
C_v(t) = \langle v_i(0)\cdot v_i(t)\rangle. |
| 585 |
> |
C_v(t) = \frac{\langle v_i(0)\cdot v_i(t)\rangle}{\langle v_i(0)\cdot v_i(0)\rangle}. |
| 586 |
|
\label{eq:vCorr} |
| 587 |
|
\end{equation} |
| 588 |
< |
Velocity autocorrelation functions require detailed short time data |
| 589 |
< |
and long trajectories for good statistics, thus velocity information |
| 590 |
< |
was saved every 5 fs over 100 ps trajectories. The power spectrum |
| 591 |
< |
($I(\omega)$) is obtained via Fourier transform of the autocorrelation |
| 592 |
< |
function |
| 588 |
> |
Velocity autocorrelation functions require detailed short time data, |
| 589 |
> |
thus velocity information was saved every 2 fs over 10 ps |
| 590 |
> |
trajectories. Because the NaCl crystal is composed of two different |
| 591 |
> |
atom types, the average of the two resulting velocity autocorrelation |
| 592 |
> |
functions was used for comparisons. |
| 593 |
> |
|
| 594 |
> |
\subsection{Long-Time and Collective Motion}\label{sec:LongTimeMethods} |
| 595 |
> |
Evaluation of the long-time dynamics of charged systems was performed |
| 596 |
> |
by considering the NaCl crystal system, again while using a subset of |
| 597 |
> |
the best performing pairwise methods. To enhance the atomic motion, |
| 598 |
> |
these crystals were equilibrated at 1000 K, near the experimental |
| 599 |
> |
$T_m$ for NaCl. Simulations were performed under the microcanonical |
| 600 |
> |
ensemble, and velocity information was saved every 5 fs over 100 ps |
| 601 |
> |
trajectories. The power spectrum ($I(\omega)$) was obtained via |
| 602 |
> |
Fourier transform of the velocity autocorrelation function |
| 603 |
|
\begin{equation} |
| 604 |
|
I(\omega) = \frac{1}{2\pi}\int^{\infty}_{-\infty}C_v(t)e^{-i\omega t}dt, |
| 605 |
|
\label{eq:powerSpec} |
| 606 |
|
\end{equation} |
| 607 |
< |
where the frequency, $\omega=0,\ 1,\ ...,\ N-1$. |
| 607 |
> |
where the frequency, $\omega=0,\ 1,\ ...,\ N-1$. Again, because the |
| 608 |
> |
NaCl crystal is composed of two different atom types, the average of |
| 609 |
> |
the two resulting power spectra was used for comparisons. |
| 610 |
|
|
| 611 |
|
\subsection{Representative Simulations}\label{sec:RepSims} |
| 612 |
|
A variety of common and representative simulations were analyzed to |
| 892 |
|
up to 0.2 \AA$^{-1}$ proves to be beneficial, but damping is arguably |
| 893 |
|
unnecessary when using the {\sc sf} method. |
| 894 |
|
|
| 895 |
< |
\subsection{Collective Motion: Power Spectra of NaCl Crystals} |
| 895 |
> |
\subsection{Short-Time Dynamics: Velocity Autocorrelation Functions of NaCl Crystals} |
| 896 |
|
|
| 897 |
|
In the previous studies using a {\sc sf} variant of the damped |
| 898 |
|
Wolf coulomb potential, the structure and dynamics of water were |
| 907 |
|
|
| 908 |
|
\begin{figure} |
| 909 |
|
\centering |
| 910 |
+ |
\includegraphics[width = \linewidth]{./vCorrPlot.pdf} |
| 911 |
+ |
\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.} |
| 912 |
+ |
\label{fig:vCorrPlot} |
| 913 |
+ |
\end{figure} |
| 914 |
+ |
|
| 915 |
+ |
The short-time decays through the first collision are nearly identical |
| 916 |
+ |
in figure \ref{fig:vCorrPlot}, but the peaks and troughs of the |
| 917 |
+ |
functions show how the methods differ. The undamped {\sc sf} method |
| 918 |
+ |
has deeper troughs (see inset in Fig. \ref{fig:vCorrPlot}) and higher |
| 919 |
+ |
peaks than any of the other methods. As the damping function is |
| 920 |
+ |
increased, these peaks are smoothed out, and approach the SPME |
| 921 |
+ |
curve. The damping acts as a distance dependent Gaussian screening of |
| 922 |
+ |
the point charges for the pairwise summation methods; thus, the |
| 923 |
+ |
collisions are more elastic in the undamped {\sc sf} potental, and the |
| 924 |
+ |
stiffness of the potential is diminished as the electrostatic |
| 925 |
+ |
interactions are softened by the damping function. With $\alpha$ |
| 926 |
+ |
values of 0.2 \AA$^{-1}$, the {\sc sf} and {\sc sp} functions are |
| 927 |
+ |
nearly identical and track the SPME features quite well. This is not |
| 928 |
+ |
too surprising in that the differences between the {\sc sf} and {\sc |
| 929 |
+ |
sp} potentials are mitigated with increased damping. However, this |
| 930 |
+ |
appears to indicate that once damping is utilized, the form of the |
| 931 |
+ |
potential seems to play a lesser role in the crystal dynamics. |
| 932 |
+ |
|
| 933 |
+ |
\subsection{Collective Motion: Power Spectra of NaCl Crystals} |
| 934 |
+ |
|
| 935 |
+ |
The short time dynamics were extended to evaluate how the differences |
| 936 |
+ |
between the methods affect the collective long-time motion. The same |
| 937 |
+ |
electrostatic summation methods were used as in the short time |
| 938 |
+ |
velocity autocorrelation function evaluation, but the trajectories |
| 939 |
+ |
were sampled over a much longer time. The power spectra of the |
| 940 |
+ |
resulting velocity autocorrelation functions were calculated and are |
| 941 |
+ |
displayed in figure \ref{fig:methodPS}. |
| 942 |
+ |
|
| 943 |
+ |
\begin{figure} |
| 944 |
+ |
\centering |
| 945 |
|
\includegraphics[width = \linewidth]{./spectraSquare.pdf} |
| 946 |
|
\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.} |
| 947 |
|
\label{fig:methodPS} |
| 948 |
|
\end{figure} |
| 949 |
|
|
| 950 |
< |
Figure \ref{fig:methodPS} shows the power spectra for the NaCl |
| 951 |
< |
crystals (from averaged Na and Cl ion velocity autocorrelation |
| 952 |
< |
functions) using the stated electrostatic summation methods. While |
| 953 |
< |
high frequency peaks of all the spectra overlap, showing the same |
| 954 |
< |
general features, the low frequency region shows how the summation |
| 955 |
< |
methods differ. Considering the low-frequency inset (expanded in the |
| 956 |
< |
upper frame of figure \ref{fig:dampInc}), at frequencies below 100 |
| 957 |
< |
cm$^{-1}$, the correlated motions are blue-shifted when using undamped |
| 958 |
< |
or weakly damped {\sc sf}. When using moderate damping ($\alpha |
| 959 |
< |
= 0.2$ \AA$^{-1}$) both the {\sc sf} and {\sc sp} |
| 960 |
< |
methods give near identical correlated motion behavior as the Ewald |
| 961 |
< |
method (which has a damping value of 0.3119). The damping acts as a |
| 962 |
< |
distance dependent Gaussian screening of the point charges for the |
| 963 |
< |
pairwise summation methods. This weakening of the electrostatic |
| 964 |
< |
interaction with distance explains why the long-ranged correlated |
| 920 |
< |
motions are at lower frequencies for the moderately damped methods |
| 921 |
< |
than for undamped or weakly damped methods. To see this effect more |
| 922 |
< |
clearly, we show how damping strength affects a simple real-space |
| 923 |
< |
electrostatic potential, |
| 950 |
> |
While high frequency peaks of the spectra in this figure overlap, |
| 951 |
> |
showing the same general features, the low frequency region shows how |
| 952 |
> |
the summation methods differ. Considering the low-frequency inset |
| 953 |
> |
(expanded in the upper frame of figure \ref{fig:dampInc}), at |
| 954 |
> |
frequencies below 100 cm$^{-1}$, the correlated motions are |
| 955 |
> |
blue-shifted when using undamped or weakly damped {\sc sf}. When |
| 956 |
> |
using moderate damping ($\alpha = 0.2$ \AA$^{-1}$) both the {\sc sf} |
| 957 |
> |
and {\sc sp} methods give near identical correlated motion behavior as |
| 958 |
> |
the Ewald method (which has a damping value of 0.3119). This |
| 959 |
> |
weakening of the electrostatic interaction with increased damping |
| 960 |
> |
explains why the long-ranged correlated motions are at lower |
| 961 |
> |
frequencies for the moderately damped methods than for undamped or |
| 962 |
> |
weakly damped methods. To see this effect more clearly, we show how |
| 963 |
> |
damping strength alone affects a simple real-space electrostatic |
| 964 |
> |
potential, |
| 965 |
|
\begin{equation} |
| 966 |
|
V_\textrm{damped}=\sum^N_i\sum^N_{j\ne i}q_iq_j\left[\frac{\textrm{erfc}({\alpha r})}{r}\right]S(r), |
| 967 |
|
\end{equation} |
| 976 |
|
shift to higher frequency in exponential fashion. Though not shown, |
| 977 |
|
the spectrum for the simple undamped electrostatic potential is |
| 978 |
|
blue-shifted such that the lowest frequency peak resides near 325 |
| 979 |
< |
cm$^{-1}$. In light of these results, the undamped {\sc sf} |
| 980 |
< |
method producing low-lying motion peaks within 10 cm$^{-1}$ of SPME is |
| 981 |
< |
quite respectable; however, it appears as though moderate damping is |
| 982 |
< |
required for accurate reproduction of crystal dynamics. |
| 979 |
> |
cm$^{-1}$. In light of these results, the undamped {\sc sf} method |
| 980 |
> |
producing low-lying motion peaks within 10 cm$^{-1}$ of SPME is quite |
| 981 |
> |
respectable and shows that the shifted force procedure accounts for |
| 982 |
> |
most of the effect afforded through use of the Ewald summation. |
| 983 |
> |
However, it appears as though moderate damping is required for |
| 984 |
> |
accurate reproduction of crystal dynamics. |
| 985 |
|
\begin{figure} |
| 986 |
|
\centering |
| 987 |
|
\includegraphics[width = \linewidth]{./comboSquare.pdf} |
| 1025 |
|
standard by which these simple pairwise sums are judged. However, |
| 1026 |
|
these results do suggest that in the typical simulations performed |
| 1027 |
|
today, the Ewald summation may no longer be required to obtain the |
| 1028 |
< |
level of accuracy most researcher have come to expect |
| 1028 |
> |
level of accuracy most researchers have come to expect |
| 1029 |
|
|
| 1030 |
|
\section{Acknowledgments} |
| 1031 |
|
\newpage |
