389 |
|
approach than they are in the Taylor-shifted method. |
390 |
|
|
391 |
|
For the gradient shifted (GSF) method with the undamped kernel, |
392 |
< |
$v_{21}(r) = -\frac{1}{r^3} + \frac{3 r}{r_c^4} + \frac{4}{r_c^3}$ and |
392 |
> |
$v_{21}(r) = -\frac{1}{r^3} - \frac{3 r}{r_c^4} + \frac{4}{r_c^3}$ and |
393 |
|
$v_{22}(r) = \frac{3}{r^3} + \frac{9 r}{r_c^4} - \frac{12}{r_c^3}$. |
394 |
|
Again, these functions go smoothly to zero as $r\rightarrow r_c$, and |
395 |
|
because the Taylor expansion retains only one term, they are |
896 |
|
\begin{figure} |
897 |
|
\centering |
898 |
|
\includegraphics[width=\textwidth]{newDrift_12.eps} |
899 |
< |
\caption{Analysis of the energy conservation of the real-space methods |
900 |
< |
for the SSDQ water/ion system. $\delta \mathrm{E}_1$ is the linear |
901 |
< |
drift in energy over time (in kcal/mol/particle/ns) and $\delta |
902 |
< |
\mathrm{E}_0$ is the standard deviation of energy fluctuations |
903 |
< |
around this drift (in kcal/mol/particle). Points that appear in the |
904 |
< |
green region at the bottom exhibit better energy conservation than |
905 |
< |
would be obtained using common parameters for Ewald-based |
906 |
< |
electrostatics.\label{fig:energyDrift}} |
899 |
> |
\caption{Energy conservation of the real-space methods for the SSDQ |
900 |
> |
water/ion system. $\delta \mathrm{E}_1$ is the linear drift in |
901 |
> |
energy over time (in kcal/mol/particle/ns) and $\delta |
902 |
> |
\mathrm{E}_0$ is the standard deviation of energy fluctuations |
903 |
> |
around this drift (in kcal/mol/particle). Points that appear in |
904 |
> |
the green region at the bottom exhibit better energy conservation |
905 |
> |
than would be obtained using common parameters for Ewald-based |
906 |
> |
electrostatics.\label{fig:energyDrift}} |
907 |
|
\end{figure} |
908 |
|
|
909 |
|
\subsection{Reproduction of Structural \& Dynamical Features\label{sec:structure}} |
931 |
|
treatment.\label{fig:gofr}} |
932 |
|
\end{figure} |
933 |
|
|
934 |
< |
There is a very slight overstructuring of the first solvation shell |
935 |
< |
when using when using TSF at lower values of the damping coefficient |
936 |
< |
($\alpha \le 0.1$) or when using undamped GSF. With moderate damping, |
937 |
< |
GSF and SP produce pair distributions that are identical (within |
938 |
< |
numerical noise) to their Ewald counterparts. |
934 |
> |
There is a minor overstructuring of the first solvation shell when |
935 |
> |
using TSF or when overdamping with any of the real-space methods. |
936 |
> |
With moderate damping, GSF and SP produce pair distributions that are |
937 |
> |
identical (within numerical noise) to their Ewald counterparts. The |
938 |
> |
degree of overstructuring can be measured most easily using the |
939 |
> |
coordination number, |
940 |
> |
\begin{equation} |
941 |
> |
n_C = 4\pi\rho \int_{0}^{a}r^2\text{g}(r)dr, |
942 |
> |
\end{equation} |
943 |
> |
where $\rho$ is the number density of the site-site pair interactions, |
944 |
> |
$a$ and is the radial location of the minima following the first peak |
945 |
> |
in $g(r)$ ($a = 4.2$ \AA for the SSDQ water/ion system). The |
946 |
> |
coordination number is shown as a function of the damping coefficient |
947 |
> |
for all of the real space methods in Fig. \ref{fig:Props}. |
948 |
|
|
949 |
< |
A structural property that is a more demanding test of modified |
950 |
< |
electrostatics is the mean value of the electrostatic energy $\langle |
951 |
< |
U_\mathrm{elect} \rangle / N$ which is obtained by sampling the |
952 |
< |
liquid-state configurations experienced by a liquid evolving entirely |
953 |
< |
under the influence of each of the methods. In table \ref{tab:Props} |
954 |
< |
we demonstrate how $\langle U_\mathrm{elect} \rangle / N$ varies with |
955 |
< |
the damping parameter, $\alpha$, for each of the methods. |
949 |
> |
A more demanding test of modified electrostatics is the average value |
950 |
> |
of the electrostatic energy $\langle U_\mathrm{elect} \rangle / N$ |
951 |
> |
which is obtained by sampling the liquid-state configurations |
952 |
> |
experienced by a liquid evolving entirely under the influence of each |
953 |
> |
of the methods. In fig \ref{fig:Props} we demonstrate how $\langle |
954 |
> |
U_\mathrm{elect} \rangle / N$ varies with the damping parameter, |
955 |
> |
$\alpha$, for each of the methods. |
956 |
|
|
957 |
|
As in the crystals studied in the first paper, damping is important |
958 |
|
for converging the mean electrostatic energy values, particularly for |
959 |
|
the two shifted force methods (GSF and TSF). A value of $\alpha |
960 |
|
\approx 0.2$ \AA$^{-1}$ is sufficient to converge the SP and GSF |
961 |
|
energies with a cutoff of 12 \AA, while shorter cutoffs require more |
962 |
< |
dramatic damping ($\alpha \approx 0.3$ \AA$^{-1}$ for $r_c = 9$ \AA). |
962 |
> |
dramatic damping ($\alpha \approx 0.28$ \AA$^{-1}$ for $r_c = 9$ \AA). |
963 |
|
Overdamping the real-space electrostatic methods occurs with $\alpha > |
964 |
< |
0.4$, causing the estimate of the energy to drop below the Ewald |
965 |
< |
results. |
964 |
> |
0.3$, causing the estimate of the electrostatic energy to drop below |
965 |
> |
the Ewald results. |
966 |
|
|
967 |
|
These ``optimal'' values of the damping coefficient are slightly |
968 |
|
larger than what were observed for DSF electrostatics for purely |
969 |
< |
point-charge systems, although a value of $\alpha=0.18$ \AA$^{-1}$ for |
970 |
< |
$r_c = 12$\AA\ appears to be an excellent compromise for mixed |
971 |
< |
charge/multipolar systems. |
969 |
> |
point-charge systems, although the range $\alpha= 0.175 \rightarrow |
970 |
> |
0.225$ \AA$^{-1}$ for $r_c = 12$\AA\ appears to be an excellent |
971 |
> |
compromise for mixed charge/multipolar systems. |
972 |
|
|
973 |
|
To test the fidelity of the electrostatic methods at reproducing |
974 |
< |
dynamics in a multipolar liquid, it is also useful to look at |
974 |
> |
\textit{dynamics} in a multipolar liquid, it is also useful to look at |
975 |
|
transport properties, particularly the diffusion constant, |
976 |
|
\begin{equation} |
977 |
|
D = \lim_{t \rightarrow \infty} \frac{1}{6 t} \langle \left| |
979 |
|
\label{eq:diff} |
980 |
|
\end{equation} |
981 |
|
which measures long-time behavior and is sensitive to the forces on |
982 |
< |
the multipoles. For the soft dipolar fluid and the SSDQ liquid |
983 |
< |
systems, the self-diffusion constants (D) were calculated from linear |
984 |
< |
fits to the long-time portion of the mean square displacement, |
985 |
< |
$\langle r^{2}(t) \rangle$.\cite{Allen87} |
982 |
> |
the multipoles. The self-diffusion constants (D) were calculated from |
983 |
> |
linear fits to the long-time portion of the mean square displacement, |
984 |
> |
$\langle r^{2}(t) \rangle$.\cite{Allen87} In fig. \ref{fig:Props} we |
985 |
> |
demonstrate how the diffusion constant depends on the choice of |
986 |
> |
real-space methods and the damping coefficient. Both the SP and GSF |
987 |
> |
methods can obtain excellent agreement with Ewald again using moderate |
988 |
> |
damping. |
989 |
|
|
990 |
|
In addition to translational diffusion, orientational relaxation times |
991 |
|
were calculated for comparisons with the Ewald simulations and with |
1015 |
|
constant.\cite{Impey82} |
1016 |
|
|
1017 |
|
Results for the diffusion constants and orientational relaxation times |
1018 |
< |
are shown in figure \ref{tab:Props}. From this data, it is apparent |
1018 |
> |
are shown in figure \ref{fig:Props}. From this data, it is apparent |
1019 |
|
that the values for both $D$ and $\tau_2$ using the Ewald sum are |
1020 |
|
reproduced with reasonable fidelity by the GSF method. |
1021 |
|
|
1022 |
< |
The $\tau_2$ results in \ref{tab:Props} show a much greater difference |
1023 |
< |
between the real-space and the Ewald results. |
1022 |
> |
\begin{figure} |
1023 |
> |
\caption{Comparison of the structural and dynamic properties for the |
1024 |
> |
combined multipolar liquid (SSDQ water + ions) for all of the |
1025 |
> |
real-space methods with $r_c = 12$\AA. Electrostatic energies, |
1026 |
> |
$\langle U_\mathrm{elect} \rangle / N$ (in kcal / mol), |
1027 |
> |
coordination numbers, $n_C$, diffusion constants (in cm$^2$ |
1028 |
> |
s$^{-1}$), and rotational correlation times (in fs) all show |
1029 |
> |
excellent agreement with Ewald results for damping coefficients in |
1030 |
> |
the range $\alpha= 0.175 \rightarrow 0.225$ |
1031 |
> |
\AA$^{-1}$. \label{fig:Props}} |
1032 |
> |
\includegraphics[width=\textwidth]{properties.eps} |
1033 |
> |
\end{figure} |
1034 |
|
|
1013 |
– |
\begin{table} |
1014 |
– |
\caption{Comparison of the structural and dynamic properties for the |
1015 |
– |
soft dipolar liquid test for all of the real-space methods.\label{tab:Props}} |
1016 |
– |
\begin{tabular}{l|c|cccc|cccc|cccc} |
1017 |
– |
& Ewald & \multicolumn{4}{c|}{SP} & \multicolumn{4}{c|}{GSF} & \multicolumn{4}{c|}{TSF} \\ |
1018 |
– |
$\alpha$ (\AA$^{-1}$) & & |
1019 |
– |
0.0 & 0.1 & 0.2 & 0.3 & |
1020 |
– |
0.0 & 0.1 & 0.2 & 0.3 & |
1021 |
– |
0.0 & 0.1 & 0.2 & 0.3 \\ \cline{2-6}\cline{6-10}\cline{10-14} |
1035 |
|
|
1023 |
– |
$\langle U_\mathrm{elect} \rangle /N$ &&&&&&&&&&&&&\\ |
1024 |
– |
D ($10^{-4}~\mathrm{cm}^2/\mathrm{s}$)& |
1025 |
– |
470.2(6) & |
1026 |
– |
416.6(5) & |
1027 |
– |
379.6(5) & |
1028 |
– |
438.6(5) & |
1029 |
– |
476.0(6) & |
1030 |
– |
412.8(5) & |
1031 |
– |
421.1(5) & |
1032 |
– |
400.5(5) & |
1033 |
– |
437.5(6) & |
1034 |
– |
434.6(5) & |
1035 |
– |
411.4(5) & |
1036 |
– |
545.3(7) & |
1037 |
– |
459.6(6) \\ |
1038 |
– |
$\tau_2$ (fs) & |
1039 |
– |
1.136 & |
1040 |
– |
1.041 & |
1041 |
– |
1.064 & |
1042 |
– |
1.109 & |
1043 |
– |
1.211 & |
1044 |
– |
1.119 & |
1045 |
– |
1.039 & |
1046 |
– |
1.058 & |
1047 |
– |
1.21 & |
1048 |
– |
1.15 & |
1049 |
– |
1.172 & |
1050 |
– |
1.153 & |
1051 |
– |
1.125 \\ |
1052 |
– |
\end{tabular} |
1053 |
– |
\end{table} |
1054 |
– |
|
1055 |
– |
|
1036 |
|
\section{CONCLUSION} |
1037 |
|
In the first paper in this series, we generalized the |
1038 |
|
charge-neutralized electrostatic energy originally developed by Wolf |