235 |
|
|
236 |
|
We have also performed a companion set of simulations {\it without} a |
237 |
|
surrounding dielectric (i.e. using a simple cubic switching function |
238 |
< |
at the cutoff radius) and as a result we have two reparamaterizations |
239 |
< |
of SSD which could be used either with or without the Reaction Field |
238 |
> |
at the cutoff radius), and as a result we have two reparamaterizations |
239 |
> |
of SSD which could be used either with or without the reaction field |
240 |
|
turned on. |
241 |
|
|
242 |
|
Simulations to obtain the preferred density were performed in the |
254 |
|
symplectic splitting method proposed by Dullweber {\it et |
255 |
|
al.}\cite{Dullweber1997} Our reason for selecting this integrator |
256 |
|
centers on poor energy conservation of rigid body dynamics using |
257 |
< |
traditional quaternion integration.\cite{Evans77,Evans77b} While quaternions |
258 |
< |
may work well for orientational motion under NVT or NPT integrators, |
259 |
< |
our limits on energy drift in the microcanonical ensemble were quite |
260 |
< |
strict, and the drift under quaternions was substantially greater than |
261 |
< |
in the symplectic splitting method. This steady drift in the total |
262 |
< |
energy has also been observed by Kol {\it et al.}\cite{Laird97} |
257 |
> |
traditional quaternion integration.\cite{Evans77,Evans77b} In typical |
258 |
> |
microcanonical ensemble simulations, the energy drift when using |
259 |
> |
quaternions was substantially greater than when using the symplectic |
260 |
> |
splitting method (fig. \ref{timestep}). This steady drift in the |
261 |
> |
total energy has also been observed by Kol {\it et al.}\cite{Laird97} |
262 |
|
|
263 |
|
The key difference in the integration method proposed by Dullweber |
264 |
|
\emph{et al.} is that the entire rotation matrix is propagated from |
448 |
|
mean-square displacement as a function of time. The averaged results |
449 |
|
from five sets of NVE simulations are displayed in figure |
450 |
|
\ref{diffuse}, alongside experimental, SPC/E, and TIP5P |
451 |
< |
results.\cite{Gillen72,Mills73,Clancy94,Jorgensen01} |
451 |
> |
results.\cite{Gillen72,Holz00,Clancy94,Jorgensen01} |
452 |
|
|
453 |
|
\begin{figure} |
454 |
|
\begin{center} |
456 |
|
\epsfbox{betterDiffuse.epsi} |
457 |
|
\caption{Average self-diffusion constant as a function of temperature for |
458 |
|
SSD, SPC/E [Ref. \citen{Clancy94}], TIP5P [Ref. \citen{Jorgensen01}], |
459 |
< |
and Experimental data [Refs. \citen{Gillen72} and \citen{Mills73}]. Of |
459 |
> |
and Experimental data [Refs. \citen{Gillen72} and \citen{Holz00}]. Of |
460 |
|
the three water models shown, SSD has the least deviation from the |
461 |
|
experimental values. The rapidly increasing diffusion constants for |
462 |
|
TIP5P and SSD correspond to significant decrease in density at the |
918 |
|
vector can be calculated from an exponential fit in the long-time |
919 |
|
regime ($t > \tau_l^\mu$).\cite{Rothschild84} Calculation of these |
920 |
|
time constants were averaged from five detailed NVE simulations |
921 |
< |
performed at the STP density for each of the respective models. Again, |
922 |
< |
SSD/E and SSD/RF show improved behavior over SSD1 both with and |
923 |
< |
without an active reaction field. Numbers published from the original |
924 |
< |
SSD dynamics studies appear closer to the experimental values, and we |
925 |
< |
attribute this discrepancy to the implimentation of an Ewald sum |
926 |
< |
versus a reaction field. |
921 |
> |
performed at the STP density for each of the respective models. It |
922 |
> |
should be noted that the commonly cited value for $\tau_2$ of 1.9 ps |
923 |
> |
was determined from the NMR data in reference \citen{Krynicki66} at a |
924 |
> |
temperature near 34$^\circ$C.\cite{Rahman73} Because of the strong |
925 |
> |
temperature dependence of $\tau_2$, it is necessary to recalculate it |
926 |
> |
at 298 K to make proper comparisons. The value shown in Table |
927 |
> |
\ref{liquidproperties} was calculated from the same NMR data in the |
928 |
> |
fashion described in reference \citen{Krynicki66}. Again, SSD/E and |
929 |
> |
SSD/RF show improved behavior over SSD1, both with and without an |
930 |
> |
active reaction field. Turning on the reaction field leads to much |
931 |
> |
improved time constants for SSD1; however, these results also include |
932 |
> |
a corresponding decrease in system density. Numbers published from the |
933 |
> |
original SSD dynamics studies appear closer to the experimental |
934 |
> |
values, and this difference can be attributed to the use of the Ewald |
935 |
> |
sum technique versus a reaction field.\cite{Ichiye99} |
936 |
|
|
937 |
|
\subsection{Additional Observations} |
938 |
|
|