1 |
|
%\documentclass[prb,aps,twocolumn,tabularx]{revtex4} |
2 |
< |
\documentclass[11pt]{article} |
2 |
> |
\documentclass[12pt]{article} |
3 |
|
\usepackage{endfloat} |
4 |
|
\usepackage{amsmath} |
5 |
|
\usepackage{epsf} |
6 |
< |
\usepackage{berkeley} |
6 |
> |
\usepackage{times} |
7 |
> |
\usepackage{mathptm} |
8 |
|
\usepackage{setspace} |
9 |
|
\usepackage{tabularx} |
10 |
|
\usepackage{graphicx} |
99 |
|
hydrogen bonds are not as linear as those observed in ice $I_h$, |
100 |
|
however the interlocking of these subunits appears to provide |
101 |
|
significant stabilization to the overall crystal. The arrangement of |
102 |
< |
these tetramers results in surrounding open octagonal cavities that |
103 |
< |
are typically greater than 6.3 \AA\ in diameter |
104 |
< |
(Fig. \ref{iCrystal}). This open structure leads to crystals that |
105 |
< |
are typically 0.07 g/cm$^3$ less dense than ice $I_h$. |
102 |
> |
these tetramers results in octagonal cavities that are typically |
103 |
> |
greater than 6.3 \AA\ in diameter (Fig. \ref{iCrystal}). This open |
104 |
> |
structure leads to crystals that are typically 0.07 g/cm$^3$ less |
105 |
> |
dense than ice $I_h$. |
106 |
|
|
107 |
|
\begin{figure} |
108 |
|
\centering |
192 |
|
of and rotation around the principle axis of the molecule |
193 |
|
respectively. These spring constants are typically calculated from |
194 |
|
the mean-square displacements of water molecules in an unrestrained |
195 |
< |
ice crystal at 200 K. For these studies, $K_\mathrm{r} = 4.29$ kcal |
196 |
< |
mol$^{-1}$, $K_\theta\ = 13.88$ kcal mol$^{-1}$, and $K_\omega\ = |
197 |
< |
17.75$ kcal mol$^{-1}$. It is clear from Fig. \ref{waterSpring} that |
198 |
< |
the values of $\theta$ range from $0$ to $\pi$, while $\omega$ ranges |
199 |
< |
from $-\pi$ to $\pi$. The partition function for a molecular crystal |
200 |
< |
restrained in this fashion can be evaluated analytically, and the |
201 |
< |
Helmholtz Free Energy ({\it A}) is given by |
195 |
> |
ice crystal at 200 K. For these studies, $K_\mathrm{v} = 4.29$ kcal |
196 |
> |
mol$^{-1}$ \AA$^{-2}$, $K_\theta\ = 13.88$ kcal mol$^{-1}$ rad$^{-2}$, |
197 |
> |
and $K_\omega\ = 17.75$ kcal mol$^{-1}$ rad$^{-2}$. It is clear from |
198 |
> |
Fig. \ref{waterSpring} that the values of $\theta$ range from $0$ to |
199 |
> |
$\pi$, while $\omega$ ranges from $-\pi$ to $\pi$. The partition |
200 |
> |
function for a molecular crystal restrained in this fashion can be |
201 |
> |
evaluated analytically, and the Helmholtz Free Energy ({\it A}) is |
202 |
> |
given by |
203 |
|
\begin{eqnarray} |
204 |
|
A = E_m\ -\ kT\ln \left (\frac{kT}{h\nu}\right )^3&-&kT\ln \left |
205 |
|
[\pi^\frac{1}{2}\left (\frac{8\pi^2I_\mathrm{A}kT}{h^2}\right |
390 |
|
conditions, such as the density in fixed-volume simulations, can |
391 |
|
influence the polymorph expressed upon crystallization. |
392 |
|
|
393 |
< |
So what is the preferred solid polymorph for simulated water? The |
392 |
< |
answer appears to be dependent both on the conditions and the model |
393 |
< |
used. In the case of short cutoffs without a long-range interaction |
394 |
< |
correction, Ice-{\it i} and Ice-{\it i}$^\prime$ have the lowest free |
395 |
< |
energy of the studied polymorphs with all the models. Ideally, |
396 |
< |
crystallization of each model under constant pressure conditions, as |
397 |
< |
was done with SSD/E, would aid in the identification of their |
398 |
< |
respective preferred structures. This work, however, helps illustrate |
399 |
< |
how studies involving one specific model can lead to insight about |
400 |
< |
important behavior of others. In general, the above results support |
401 |
< |
the finding that the Ice-{\it i} polymorph is a stable crystal |
402 |
< |
structure that should be considered when studying the phase behavior |
403 |
< |
of water models. |
393 |
> |
\section{Conclusions} |
394 |
|
|
395 |
+ |
In this work, thermodynamic integration was used to determine the |
396 |
+ |
absolute free energies of several ice polymorphs. The new polymorph, |
397 |
+ |
Ice-{\it i} was observed to be the stable crystalline state for {\it |
398 |
+ |
all} the water models when using a 9.0 \AA\ cutoff. However, the free |
399 |
+ |
energy partially depends on simulation conditions (particularly on the |
400 |
+ |
choice of long range correction method). Regardless, Ice-{\it i} was |
401 |
+ |
still observered to be a stable polymorph for all of the studied water |
402 |
+ |
models. |
403 |
+ |
|
404 |
+ |
So what is the preferred solid polymorph for simulated water? As |
405 |
+ |
indicated above, the answer appears to be dependent both on the |
406 |
+ |
conditions and the model used. In the case of short cutoffs without a |
407 |
+ |
long-range interaction correction, Ice-{\it i} and Ice-{\it |
408 |
+ |
i}$^\prime$ have the lowest free energy of the studied polymorphs with |
409 |
+ |
all the models. Ideally, crystallization of each model under constant |
410 |
+ |
pressure conditions, as was done with SSD/E, would aid in the |
411 |
+ |
identification of their respective preferred structures. This work, |
412 |
+ |
however, helps illustrate how studies involving one specific model can |
413 |
+ |
lead to insight about important behavior of others. |
414 |
+ |
|
415 |
|
We also note that none of the water models used in this study are |
416 |
|
polarizable or flexible models. It is entirely possible that the |
417 |
|
polarizability of real water makes Ice-{\it i} substantially less |
430 |
|
situation for possible observation. These include the negative |
431 |
|
pressure or stretched solid regime, small clusters in vacuum |
432 |
|
deposition environments, and in clathrate structures involving small |
433 |
< |
non-polar molecules. For experimental comparison purposes, example |
434 |
< |
$g_{OO}(r)$ and $S(\vec{q})$ plots were generated for the two Ice-{\it |
435 |
< |
i} variants (along with example ice $I_h$ and $I_c$ plots) at 77K, and |
436 |
< |
they are shown in figures \ref{fig:gofr} and \ref{fig:sofq} |
437 |
< |
respectively. |
433 |
> |
non-polar molecules. For the purpose of comparison with experimental |
434 |
> |
results, we have calculated the oxygen-oxygen pair correlation |
435 |
> |
function, $g_{OO}(r)$, and the structure factor, $S(\vec{q})$ for the |
436 |
> |
two Ice-{\it i} variants (along with example ice $I_h$ and $I_c$ |
437 |
> |
plots) at 77K, and they are shown in figures \ref{fig:gofr} and |
438 |
> |
\ref{fig:sofq} respectively. It is interesting to note that the |
439 |
> |
structure factors for Ice-{\it i}$^\prime$ and Ice-I$_c$ are quite similar. |
440 |
> |
The primary differences are small peaks at 1.125, 2.29, and 2.53 |
441 |
> |
\AA${-1}$, so particular attention to these regions would be needed |
442 |
> |
to identify the new {\it i}$^\prime$ variant from the I$_{c}$ variant. |
443 |
|
|
444 |
|
\begin{figure} |
445 |
|
\centering |
470 |
|
|
471 |
|
\newpage |
472 |
|
|
473 |
< |
\bibliographystyle{jcp} |
473 |
> |
\bibliographystyle{achemso} |
474 |
|
\bibliography{iceiPaper} |
475 |
|
|
476 |
|
|