23 |
|
|
24 |
|
\begin{document} |
25 |
|
|
26 |
< |
\title{On the temperature dependent structural and transport properties of the soft sticky dipole (SSD) and related single point water models} |
26 |
> |
\title{On the temperature dependent properties of the soft sticky dipole (SSD) and related single point water models} |
27 |
|
|
28 |
|
\author{Christopher J. Fennell and J. Daniel Gezelter{\thefootnote} |
29 |
|
\footnote[1]{Corresponding author. \ Electronic mail: gezelter@nd.edu}} |
160 |
|
systems, the Ewald summation and even particle-mesh Ewald become |
161 |
|
computational burdens with their respective ideal $N^\frac{3}{2}$ and |
162 |
|
$N\log N$ calculation scaling orders for $N$ particles.\cite{Darden99} |
163 |
< |
|
164 |
< |
Up to this point, a detailed look at the model's structure and ion |
165 |
< |
solvation abilities has been performed.\cite{Ichiye96} In addition, a |
166 |
< |
thorough investigation of the dynamic properties of SSD was performed |
167 |
< |
by Chandra and Ichiye focusing on translational and orientational |
168 |
< |
properties at 298 K.\cite{Ichiye99} This study focuses on determining |
169 |
< |
the density maximum for SSD utilizing both microcanonical and |
170 |
< |
isobaric-isothermal ensemble molecular dynamics, while using the |
171 |
< |
reaction field method for handling long-ranged dipolar interactions. A |
172 |
< |
reaction field method has been previously implemented in Monte Carlo |
173 |
< |
simulations by Liu and Ichiye in order to study the static dielectric |
174 |
< |
constant for the model.\cite{Ichiye96b} This paper will expand the |
175 |
< |
scope of these original simulations to look on how the reaction field |
176 |
< |
affects the physical and dynamic properties of SSD systems. |
163 |
> |
In applying this water model in these types of systems, it would be |
164 |
> |
useful to know its properties and behavior with the more |
165 |
> |
computationally efficient reaction field (RF) technique, and even with |
166 |
> |
a cutoff that lacks any form of long range correction. This study |
167 |
> |
addresses these issues by looking at the structural and transport |
168 |
> |
behavior of SSD over a variety of temperatures, with the purpose of |
169 |
> |
utilizing the RF correction technique. Towards the end, we suggest |
170 |
> |
alterations to the parameters that result in more water-like |
171 |
> |
behavior. It should be noted that in a recent publication, some the |
172 |
> |
original investigators of the SSD water model have put forth |
173 |
> |
adjustments to the original SSD water model to address abnormal |
174 |
> |
density behavior (also observed here), calling the corrected model |
175 |
> |
SSD1.\cite{Ichiye03} This study will consider this new model's |
176 |
> |
behavior as well, and hopefully improve upon its depiction of water |
177 |
> |
under conditions without the Ewald Sum. |
178 |
|
|
179 |
|
\section{Methods} |
180 |
|
|
233 |
|
nine elements long as opposed to 3 or 4 elements for Euler angles and |
234 |
|
quaternions respectively. System memory has become much less of an |
235 |
|
issue in recent times, and this has resulted in substantial benefits |
236 |
< |
in energy conservation. There is still the issue of an additional 5 or |
237 |
< |
6 additional elements for describing the orientation of each particle, |
238 |
< |
which will increase dump files substantially. Simply translating the |
239 |
< |
rotation matrix into its component Euler angles or quaternions for |
240 |
< |
storage purposes relieves this burden. |
236 |
> |
in energy conservation. There is still the issue of 5 or 6 additional |
237 |
> |
elements for describing the orientation of each particle, which will |
238 |
> |
increase dump files substantially. Simply translating the rotation |
239 |
> |
matrix into its component Euler angles or quaternions for storage |
240 |
> |
purposes relieves this burden. |
241 |
|
|
242 |
|
The symplectic splitting method allows for Verlet style integration of |
243 |
|
both linear and angular motion of rigid bodies. In the integration |
245 |
|
evaluations to update the rotation matrix.\cite{Dullweber1997} These |
246 |
|
matrix rotations end up being more costly computationally than the |
247 |
|
simpler arithmetic quaternion propagation. On average, a 1000 SSD |
248 |
< |
particle simulation shows a 7\% increase in simulation time using the |
248 |
> |
particle simulation shows a 7\% increase in computation time using the |
249 |
|
symplectic step method in place of quaternions. This cost is more than |
250 |
|
justified when comparing the energy conservation of the two methods as |
251 |
< |
illustrated in figure \ref{timestep}. |
251 |
> |
illustrated in figure \ref{timestep}. |
252 |
|
|
253 |
|
\begin{figure} |
254 |
|
\includegraphics[width=61mm, angle=-90]{timeStep.epsi} |
270 |
|
with the quaternion method showing a slight energy drift over time in |
271 |
|
the 0.5 fs time step simulation. At time steps of 1 and 2 fs, the |
272 |
|
energy conservation benefits of the symplectic step method are clearly |
273 |
< |
demonstrated. |
273 |
> |
demonstrated. Thus, while maintaining the same degree of energy |
274 |
> |
conservation, one can take considerably longer time steps, leading to |
275 |
> |
an overall reduction in computation time. |
276 |
|
|
277 |
|
Energy drift in these SSD particle simulations was unnoticeable for |
278 |
|
time steps up to three femtoseconds. A slight energy drift on the |