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} |
21 |
|
\begin{document} |
22 |
|
|
23 |
|
\title{Computational free energy studies of a new ice polymorph which |
24 |
< |
exhibits greater stability than Ice $I_h$} |
24 |
> |
exhibits greater stability than Ice I$_h$} |
25 |
|
|
26 |
|
\author{Christopher J. Fennell and J. Daniel Gezelter \\ |
27 |
|
Department of Chemistry and Biochemistry\\ |
96 |
|
tetramers form a crystal structure similar in appearance to a recent |
97 |
|
two-dimensional surface tessellation simulated on silica.\cite{Yang04} |
98 |
|
As expected in an ice crystal constructed of water tetramers, the |
99 |
< |
hydrogen bonds are not as linear as those observed in ice $I_h$, |
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 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$. |
105 |
> |
dense than ice I$_h$. |
106 |
|
|
107 |
|
\begin{figure} |
108 |
|
\centering |
119 |
|
\includegraphics[width=\linewidth]{orderedIcei.eps} |
120 |
|
\caption{A rendering of a proton ordered crystal of Ice-{\it i} looking |
121 |
|
down the (001) crystal face. The presence of large octagonal pores |
122 |
< |
leads to a polymorph that is less dense than ice $I_h$.} |
122 |
> |
leads to a polymorph that is less dense than ice I$_h$.} |
123 |
|
\label{iCrystal} |
124 |
|
\end{figure} |
125 |
|
|
132 |
|
contributions to the overall free energy. To address this issue, we |
133 |
|
have calculated the absolute free energy of this crystal using |
134 |
|
thermodynamic integration and compared it to the free energies of ice |
135 |
< |
$I_c$ and ice $I_h$ (the common low density ice polymorphs) and ice B |
135 |
> |
I$_c$ and ice I$_h$ (the common low density ice polymorphs) and ice B |
136 |
|
(a higher density, but very stable crystal structure observed by |
137 |
|
B\`{a}ez and Clancy in free energy studies of SPC/E).\cite{Baez95b} |
138 |
|
This work includes results for the water model from which Ice-{\it i} |
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 |
261 |
|
\section{Results and Discussion} |
262 |
|
|
263 |
|
The calculated free energies of proton-ordered variants of three low |
264 |
< |
density polymorphs ($I_h$, $I_c$, and Ice-{\it i} or Ice-{\it |
264 |
> |
density polymorphs (I$_h$, I$_c$, and Ice-{\it i} or Ice-{\it |
265 |
|
i}$^\prime$) and the stable higher density ice B are listed in Table |
266 |
|
\ref{freeEnergy}. Ice B was included because it has been |
267 |
|
shown to be a minimum free energy structure for SPC/E at ambient |
276 |
|
All other models have similar structure, although the crossing points |
277 |
|
between the phases move to different temperatures and pressures as |
278 |
|
indicated from the transition temperatures in Table \ref{freeEnergy}. |
279 |
< |
It is interesting to note that ice $I_h$ (and ice $I_c$ for that |
279 |
> |
It is interesting to note that ice I$_h$ (and ice I$_c$ for that |
280 |
|
matter) do not appear in any of the phase diagrams for any of the |
281 |
|
models. For purposes of this study, ice B is representative of the |
282 |
|
dense ice polymorphs. A recent study by Sanz {\it et al.} provides |
294 |
|
Kelvin. Calculated error of the final digits is in parentheses.} |
295 |
|
\begin{tabular}{lccccccc} |
296 |
|
\hline |
297 |
< |
Water Model & $I_h$ & $I_c$ & B & Ice-{\it i} & Ice-{\it i}$^\prime$ & $T_m$ (*$T_s$) & $T_b$\\ |
297 |
> |
Water Model & I$_h$ & I$_c$ & B & Ice-{\it i} & Ice-{\it i}$^\prime$ & $T_m$ (*$T_s$) & $T_b$\\ |
298 |
|
\hline |
299 |
|
TIP3P & -11.41(2) & -11.23(3) & -11.82(3) & -12.30(3) & - & 269(4) & 357(2)\\ |
300 |
|
TIP4P & -11.84(3) & -12.04(2) & -12.08(3) & - & -12.33(3) & 266(5) & 354(2)\\ |
301 |
|
TIP5P & -11.85(3) & -11.86(2) & -11.96(2) & - & -12.29(2) & 271(4) & 337(2)\\ |
302 |
|
SPC/E & -12.87(2) & -13.05(2) & -13.26(3) & - & -13.55(2) & 296(3) & 396(2)\\ |
303 |
|
SSD/E & -11.27(2) & -11.19(4) & -12.09(2) & -12.54(2) & - & *355(2) & -\\ |
304 |
< |
SSD/RF & -11.51(2) & -11.47(2) & -12.08(3) & -12.29(2) & - & 278(4) & 349(2)\\ |
304 |
> |
SSD/RF & -11.96(2) & -11.60(2) & -12.53(3) & -12.79(2) & - & 287(4) & 382(2)\\ |
305 |
|
\end{tabular} |
306 |
|
\label{freeEnergy} |
307 |
|
\end{center} |
321 |
|
Most of the water models have melting points that compare quite |
322 |
|
favorably with the experimental value of 273 K. The unfortunate |
323 |
|
aspect of this result is that this phase change occurs between |
324 |
< |
Ice-{\it i} and the liquid state rather than ice $I_h$ and the liquid |
324 |
> |
Ice-{\it i} and the liquid state rather than ice I$_h$ and the liquid |
325 |
|
state. These results do not contradict other studies. Studies of ice |
326 |
< |
$I_h$ using TIP4P predict a $T_m$ ranging from 214 to 238 K |
326 |
> |
I$_h$ using TIP4P predict a $T_m$ ranging from 214 to 238 K |
327 |
|
(differences being attributed to choice of interaction truncation and |
328 |
|
different ordered and disordered molecular |
329 |
|
arrangements).\cite{Vlot99,Gao00,Sanz04} If the presence of ice B and |
392 |
|
|
393 |
|
\section{Conclusions} |
394 |
|
|
395 |
< |
In this report, thermodynamic integration was used to determine the |
396 |
< |
absolute free energies of several ice polymorphs. Of the studied |
397 |
< |
crystal forms, Ice-{\it i} was observed to be the stable crystalline |
398 |
< |
state for {\it all} the water models when using a 9.0 \AA\ |
399 |
< |
intermolecular interaction cutoff. Through investigation of possible |
400 |
< |
interaction truncation methods, the free energy was shown to be |
401 |
< |
partially dependent on simulation conditions; however, Ice-{\it i} was |
402 |
< |
still observered to be a stable polymorph of the studied water models. |
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 |
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. In general, the |
412 |
< |
above results support the finding that the Ice-{\it i} polymorph is a |
413 |
< |
stable crystal structure that should be considered when studying the |
414 |
< |
phase behavior of water models. |
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 |
418 |
< |
stable than ice $I_h$. However, the calculations presented above seem |
418 |
> |
stable than ice I$_h$. However, the calculations presented above seem |
419 |
|
interesting enough to communicate before the role of polarizability |
420 |
|
(or flexibility) has been thoroughly investigated. |
421 |
|
|
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 |
446 |
|
\includegraphics[width=\linewidth]{iceGofr.eps} |
447 |
< |
\caption{Radial distribution functions of ice $I_h$, $I_c$, and |
447 |
> |
\caption{Radial distribution functions of ice I$_h$, I$_c$, and |
448 |
|
Ice-{\it i} calculated from from simulations of the SSD/RF water model |
449 |
|
at 77 K. The Ice-{\it i} distribution function was obtained from |
450 |
|
simulations composed of TIP4P water.} |
454 |
|
\begin{figure} |
455 |
|
\centering |
456 |
|
\includegraphics[width=\linewidth]{sofq.eps} |
457 |
< |
\caption{Predicted structure factors for ice $I_h$, $I_c$, Ice-{\it i}, |
457 |
> |
\caption{Predicted structure factors for ice I$_h$, I$_c$, Ice-{\it i}, |
458 |
|
and Ice-{\it i}$^\prime$ at 77 K. The raw structure factors have |
459 |
|
been convoluted with a gaussian instrument function (0.075 \AA$^{-1}$ |
460 |
|
width) to compensate for the trunction effects in our finite size |
470 |
|
|
471 |
|
\newpage |
472 |
|
|
473 |
< |
\bibliographystyle{jcp} |
473 |
> |
\bibliographystyle{achemso} |
474 |
|
\bibliography{iceiPaper} |
475 |
|
|
476 |
|
|