1 |
%\documentclass[prb,aps,twocolumn,tabularx]{revtex4} |
2 |
\documentclass[preprint,aps,endfloats]{revtex4} |
3 |
%\documentclass[11pt]{article} |
4 |
%\usepackage{endfloat} |
5 |
\usepackage{amsmath} |
6 |
\usepackage{epsf} |
7 |
\usepackage{berkeley} |
8 |
%\usepackage{setspace} |
9 |
%\usepackage{tabularx} |
10 |
\usepackage{graphicx} |
11 |
%\usepackage[ref]{overcite} |
12 |
%\pagestyle{plain} |
13 |
%\pagenumbering{arabic} |
14 |
%\oddsidemargin 0.0cm \evensidemargin 0.0cm |
15 |
%\topmargin -21pt \headsep 10pt |
16 |
%\textheight 9.0in \textwidth 6.5in |
17 |
%\brokenpenalty=10000 |
18 |
|
19 |
%\renewcommand\citemid{\ } % no comma in optional reference note |
20 |
|
21 |
\begin{document} |
22 |
|
23 |
\title{A Free Energy Study of Low Temperature and Anomolous Ice} |
24 |
|
25 |
\author{Christopher J. Fennell and J. Daniel Gezelter{\thefootnote} |
26 |
\footnote[1]{Corresponding author. \ Electronic mail: gezelter@nd.edu}} |
27 |
|
28 |
\address{Department of Chemistry and Biochemistry\\ University of Notre Dame\\ |
29 |
Notre Dame, Indiana 46556} |
30 |
|
31 |
\date{\today} |
32 |
|
33 |
%\maketitle |
34 |
%\doublespacing |
35 |
|
36 |
\begin{abstract} |
37 |
\end{abstract} |
38 |
|
39 |
\maketitle |
40 |
|
41 |
\newpage |
42 |
|
43 |
%\narrowtext |
44 |
|
45 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
46 |
% BODY OF TEXT |
47 |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
48 |
|
49 |
\section{Introduction} |
50 |
|
51 |
\section{Methods} |
52 |
|
53 |
Canonical ensemble (NVT) molecular dynamics calculations were |
54 |
performed using the OOPSE (Object-Oriented Parallel Simulation Engine) |
55 |
molecular mechanics package. All molecules were treated as rigid |
56 |
bodies, with orientational motion propogated using the symplectic DLM |
57 |
integration method. Details about the implementation of these |
58 |
techniques can be found in a recent publication.\cite{Meineke05} |
59 |
|
60 |
Thermodynamic integration was utilized to calculate the free energy of |
61 |
several ice crystals using the TIP3P, TIP4P, TIP5P, SPC/E, and SSD/E |
62 |
water models. Liquid state free energies at 300 and 400 K for all of |
63 |
these water models were also determined using this same technique, in |
64 |
order to determine melting points and generate phase diagrams. |
65 |
|
66 |
For the thermodynamic integration of molecular crystals, the Einstein |
67 |
Crystal is chosen as the reference state that the system is converted |
68 |
to over the course of the simulation. In an Einstein Crystal, the |
69 |
molecules are harmonically restrained at their ideal lattice locations |
70 |
and orientations. The partition function for a molecular crystal |
71 |
restrained in this fashion has been evaluated, and the Helmholtz Free |
72 |
Energy ({\it A}) is given by |
73 |
\begin{eqnarray} |
74 |
A = E_m\ -\ kT\ln \left (\frac{kT}{h\nu}\right )^3&-&kT\ln \left |
75 |
[\pi^\frac{1}{2}\left (\frac{8\pi^2I_\mathrm{A}kT}{h^2}\right |
76 |
)^\frac{1}{2}\left (\frac{8\pi^2I_\mathrm{B}kT}{h^2}\right |
77 |
)^\frac{1}{2}\left (\frac{8\pi^2I_\mathrm{C}kT}{h^2}\right |
78 |
)^\frac{1}{2}\right ] \nonumber \\ &-& kT\ln \left [\frac{kT}{2(\pi |
79 |
K_\omega K_\theta)^{\frac{1}{2}}}\exp\left |
80 |
(-\frac{kT}{2K_\theta}\right)\int_0^{\left (\frac{kT}{2K_\theta}\right |
81 |
)^\frac{1}{2}}\exp(t^2)\mathrm{d}t\right ], |
82 |
\label{ecFreeEnergy} |
83 |
\end{eqnarray} |
84 |
where $2\pi\nu = (K_\mathrm{v}/m)^{1/2}$.\cite{Baez95a} In equation |
85 |
\ref{ecFreeEnergy}, $K_\mathrm{v}$, $K_\mathrm{\theta}$, and |
86 |
$K_\mathrm{\omega}$ are the spring constants restraining translational |
87 |
motion and deflection of and rotation around the principle axis of the |
88 |
molecule respectively (Fig. \ref{waterSpring}), and $E_m$ is the |
89 |
minimum potential energy of the ideal crystal. In the case of |
90 |
molecular liquids, the ideal vapor is chosen as the target reference |
91 |
state. |
92 |
|
93 |
|
94 |
|
95 |
|
96 |
\section{Results and discussion} |
97 |
|
98 |
\section{Conclusions} |
99 |
|
100 |
\section{Acknowledgments} |
101 |
Support for this project was provided by the National Science |
102 |
Foundation under grant CHE-0134881. Computation time was provided by |
103 |
the Notre Dame Bunch-of-Boxes (B.o.B) computer cluster under NSF grant |
104 |
DMR-0079647. |
105 |
|
106 |
\newpage |
107 |
|
108 |
\bibliographystyle{jcp} |
109 |
\bibliography{iceiPaper} |
110 |
|
111 |
|
112 |
\end{document} |