--- trunk/iceWater/iceWater.tex	2013/09/09 12:45:25	3950
+++ trunk/iceWater/iceWater.tex	2013/10/23 17:29:18	3966
@@ -1,47 +1,81 @@
-\documentclass[journal = jpccck, manuscript = article]{achemso}
-\setkeys{acs}{usetitle = true}
-\usepackage{achemso}
-\usepackage{natbib}
+\documentclass[11pt]{article}
+\usepackage{amsmath}
+\usepackage{amssymb}
+\usepackage{setspace}
+%\usepackage{endfloat}
+\usepackage{caption}
+%\usepackage{epsf}
+%\usepackage{tabularx}
+\usepackage{graphicx}
 \usepackage{multirow}
 \usepackage{wrapfig}
-\usepackage{fixltx2e}
-%\mciteErrorOnUnknownfalse
+%\usepackage{booktabs}
+%\usepackage{bibentry}
+%\usepackage{mathrsfs}
+%\usepackage[ref]{overcite}
+\usepackage[square, comma, sort&compress]{natbib}
+\usepackage{url}
+\pagestyle{plain} \pagenumbering{arabic} \oddsidemargin 0.0cm
+\evensidemargin 0.0cm \topmargin -21pt \headsep 10pt \textheight
+9.0in \textwidth 6.5in \brokenpenalty=10000
 
+% double space list of tables and figures
+%\AtBeginDelayedFloats{\renewcommand{\baselinestretch}{1.66}}
+\setlength{\abovecaptionskip}{20 pt}
+\setlength{\belowcaptionskip}{30 pt}
+
+%\renewcommand\citemid{\ } % no comma in optional referenc note
+\bibpunct{}{}{,}{s}{}{;}
+\bibliographystyle{aip}
+
+
+% \documentclass[journal = jpccck, manuscript = article]{achemso}
+% \setkeys{acs}{usetitle = true}
+% \usepackage{achemso}
+% \usepackage{natbib}
+% \usepackage{multirow}
+% \usepackage{wrapfig}
+% \usepackage{fixltx2e}
+
 \usepackage[version=3]{mhchem}  % this is a great package for formatting chemical reactions
 \usepackage{url}
 
+
+\begin{document}
+
 \title{Simulations of solid-liquid friction at ice-I$_\mathrm{h}$ /
   water interfaces}
 
-\author{P. B. Louden}
-\author{J. Daniel Gezelter}
-\email{gezelter@nd.edu}
-\affiliation[University of Notre Dame]{251 Nieuwland Science Hall\\
-  Department of Chemistry and Biochemistry\\ University of Notre
-  Dame\\ Notre Dame, Indiana 46556}
+\author{Patrick B. Louden and J. Daniel 
+Gezelter\footnote{Corresponding author. \ Electronic mail:
+  gezelter@nd.edu} \\
+Department of Chemistry and Biochemistry,\\
+University of Notre Dame\\
+Notre Dame, Indiana 46556}
 
-\keywords{}
+\date{\today}
+\maketitle
+\begin{doublespace}
 
-\begin{document}
 
 \begin{abstract}
   We have investigated the structural and dynamic properties of the
-  basal and prismatic facets of an ice I$_\mathrm{h}$ / water
-  interface when the solid phase is being drawn through the liquid
+  basal and prismatic facets of the ice I$_\mathrm{h}$ / water
+  interface when the solid phase is drawn through the liquid
   (i.e. sheared relative to the fluid phase). To impose the shear, we
   utilized a velocity-shearing and scaling (VSS) approach to reverse
   non-equilibrium molecular dynamics (RNEMD).  This method can create
   simultaneous temperature and velocity gradients and allow the
   measurement of transport properties at interfaces.  The interfacial
-  width was found to be independent of relative velocity of the ice
-  and liquid layers over a wide range of shear rates.  Decays of
-  molecular orientational time correlation functions for gave very
-  similar estimates for the width of the interfaces, although the
-  short- and longer-time decay components of the orientational
-  correlation functions behave differently closer to the interface.
-  Although both facets of ice are in ``stick'' boundary conditions in
-  liquid water, the solid-liquid friction coefficient was found to be
-  different for the basal and prismatic facets of ice.
+  width was found to be independent of the relative velocity of the
+  ice and liquid layers over a wide range of shear rates.  Decays of
+  molecular orientational time correlation functions gave similar
+  estimates for the width of the interfaces, although the short- and
+  longer-time decay components behave differently closer to the
+  interface.  Although both facets of ice are in ``stick'' boundary
+  conditions in liquid water, the solid-liquid friction coefficients
+  were found to be significantly different for the basal and prismatic
+  facets of ice.
 \end{abstract}
 
 \newpage
@@ -67,15 +101,15 @@ interface for the SPC\cite{Karim90}, SPC/E\cite{Gay02,
 of quiescent ice/water interfaces utilizing both theory and
 experiment.  Haymet \emph{et al.} have done extensive work on ice I$_\mathrm{h}$,
 including characterizing and determining the width of the ice/water
-interface for the SPC\cite{Karim90}, SPC/E\cite{Gay02,Bryk02}, CF1\cite{Hayward01,Hayward02}, and TIP4P\cite{Karim88} models for
+interface for the SPC,\cite{Karim90} SPC/E,\cite{Gay02,Bryk02} CF1,\cite{Hayward01,Hayward02} and TIP4P~\cite{Karim88} models for
 water.
 More recently, Haymet \emph{et al.} have investigated the effects
 cations and anions have on crystal
-nucleation\cite{Bryk04,Smith05,Wilson08,Wilson10}. Nada \emph{et al.}
+nucleation.\cite{Bryk04,Smith05,Wilson08,Wilson10} Nada \emph{et al.}
 have also studied ice/water
 interfaces,\cite{Nada95,Nada00,Nada03,Nada12} and have found that the
 differential growth rates of the facets of ice I$_\mathrm{h}$ are due to the the
-reordering of the hydrogen bonding network\cite{Nada05}.
+reordering of the hydrogen bonding network.\cite{Nada05}
 
 The movement of liquid water over the facets of ice has been less
 thoroughly studied than the quiescent surfaces. This process is
@@ -121,10 +155,11 @@ and Oxygen lone pairs.\cite{Buch:2008fk}
 have a preference for proton ordering along strips of dangling H-atoms
 and Oxygen lone pairs.\cite{Buch:2008fk}
 
-\begin{wraptable}{r}{3.5in}
+\begin{table}[h]
+\centering
   \caption{Mapping between the Miller indices of four facets of ice in
     the $P6_3/mmc$ crystal system to the orthorhombic $P2_12_12_1$
-    system in reference \protect\cite{Hirsch04}}
+    system in reference \bibpunct{}{}{,}{n}{}{,} \protect\citep{Hirsch04}.}
 \label{tab:equiv}
 \begin{tabular}{|ccc|} \hline
  & hexagonal & orthorhombic \\
@@ -135,8 +170,7 @@ pyramidal & $\{2~0~\bar{2}~1\}$ & $\{2~0~1\}$ \\ \hlin
 secondary prism & $\{1~1~\bar{2}~0\}$ & $\{1~3~0\}$ \\
 pyramidal & $\{2~0~\bar{2}~1\}$ & $\{2~0~1\}$ \\ \hline
 \end{tabular}
-\end{wraptable}
-
+\end{table}
 For small simulated ice interfaces, it is useful to work with
 proton-ordered, but zero-dipole crystal that exposes these strips of
 dangling H-atoms and lone pairs.  When placing another material in
@@ -154,53 +188,53 @@ parameters $a = 4.49225$, $b = 7.78080$, $c = 7.33581$
 $P2_12_12_1$ system.
 
 Structure 6 from the Hirsch and Ojam\"{a}e paper has lattice
-parameters $a = 4.49225$, $b = 7.78080$, $c = 7.33581$ and two water
-molecules whose atoms reside at fractional coordinates given in table
-\ref{tab:p212121}.  To construct the basal and prismatic interfaces,
-these crystallographic coordinates were used to construct an
-orthorhombic unit cell which was then replicated in all three
+parameters $a = 4.49225$ \AA\ , $b = 7.78080$ \AA\ , $c = 7.33581$ \AA\
+and two water molecules whose atoms reside at fractional coordinates
+given in table \ref{tab:p212121}. To construct the basal and prismatic
+interfaces, these crystallographic coordinates were used to construct
+an orthorhombic unit cell which was then replicated in all three
 dimensions yielding a proton-ordered block of ice I$_{h}$. To expose
 the desired face, the orthorhombic representation was then cut along
 the ($001$) or ($100$) planes for the basal and prismatic faces
-respectively.  The resulting block was rotated so that the exposed
-faces were aligned with the $z$-axis normal to the exposed face.  The
+respectively. The resulting block was rotated so that the exposed
+faces were aligned with the $z$-axis normal to the exposed face. The
 block was then cut along two perpendicular directions in a way that
 allowed for perfect periodic replication in the $x$ and $y$ axes,
 creating a slab with either the basal or prismatic faces exposed along
-the $z$ axis.  The slab was then replicated in the $x$ and $y$
+the $z$ axis. The slab was then replicated in the $x$ and $y$
 dimensions until a desired sample size was obtained.
 
-\begin{wraptable}{r}{2.85in}
+\begin{table}[h]
+\centering
   \caption{Fractional coordinates for water in the orthorhombic
-    $P2_12_12_1$ system for ice I$_\mathrm{h}$ in reference
-    \protect\cite{Hirsch04}}
+    $P2_12_12_1$ system for ice I$_\mathrm{h}$ in reference \bibpunct{}{}{,}{n}{}{,} \protect\citep{Hirsch04}.}
 \label{tab:p212121}
 \begin{tabular}{|cccc|}  \hline
 atom type & x & y & z \\ \hline
- O & 0.75 & 0.1667 & 0.4375 \\
+ O & 0.7500 & 0.1667 & 0.4375 \\
  H & 0.5735 & 0.2202 & 0.4836 \\
  H & 0.7420 & 0.0517 & 0.4836 \\
- O & 0.25 & 0.6667 & 0.4375 \\
+ O & 0.2500 & 0.6667 & 0.4375 \\
  H & 0.2580 & 0.6693 & 0.3071 \\
  H & 0.4265 & 0.7255 & 0.4756 \\ \hline
 \end{tabular}
-\end{wraptable}
+\end{table}
 
 Our ice / water interfaces were created using a box of liquid water
 that had the same dimensions (in $x$ and $y$) as the ice block.
 Although the experimental solid/liquid coexistence temperature under
-atmospheric pressure is close to 273K, Haymet \emph{et al.} have done
+atmospheric pressure is close to 273~K, Haymet \emph{et al.} have done
 extensive work on characterizing the ice/water interface, and find
 that the coexistence temperature for simulated water is often quite a
 bit different.\cite{Karim88,Karim90,Hayward01,Bryk02,Hayward02} They
 have found that for the SPC/E water model,\cite{Berendsen87} which is
 also used in this study, the ice/water interface is most stable at
-225$\pm$5K.\cite{Bryk02} This liquid box was therefore equilibrated at
-225 K and 1 atm of pressure in the NPAT ensemble (with the $z$ axis
+225~$\pm$5~K.\cite{Bryk02} This liquid box was therefore equilibrated at
+225~K and 1~atm of pressure in the NPAT ensemble (with the $z$ axis
 allowed to fluctuate to equilibrate to the correct pressure).  The
 liquid and solid systems were combined by carving out any water
-molecule from the liquid simulation cell that was within 3 \AA\ of any
-atom in the ice slab.
+molecule from the liquid simulation cell that was within 3~\AA\ of any
+atom in the ice slab. The resulting basal system was 24 \AA\ x 36 \AA\ x 99 \AA\ with 900 SPC/E molecules in the ice slab and 1846 SPC/E molecules in the liquid phase. Similarly, the prismatic system was constructed as 36 \AA\ x 36 \AA\ x 86 \AA\ with 1000 SPC/E molecules in the ice slab and 2684 SPC/E molecules in the liquid phase.
 
 Molecular translation and orientational restraints were applied in the
 early stages of equilibration to prevent melting of the ice slab.
@@ -288,11 +322,11 @@ temperatures (90~K) with a single 1 ns simulation.\cit
 minimal, as is thermal anisotropy.  This ability to generate
 simultaneous thermal and shear fluxes has been previously utilized to
 map out the shear viscosity of SPC/E water over a wide range of
-temperatures (90~K) with a single 1 ns simulation.\cite{Kuang12}
+temperatures (90~K) with a single 1~ns simulation.\cite{Kuang12}
 
 For this work, we are using the VSS-RNEMD method primarily to generate
 a shear between the ice slab and the liquid phase, while using a weak
-thermal gradient to maintain the interface at the 225K target
+thermal gradient to maintain the interface at the 225~K target
 value. This ensures minimal melting of the bulk ice phase and allows
 us to control the exact temperature of the interface.
 
@@ -302,14 +336,14 @@ which were attempted every 50 fs. 
 dimensions.  Electrostatics were handled using the damped-shifted
 force real-space electrostatic kernel.\cite{Ewald} The systems were
 divided into 100 bins along the $z$-axis for the VSS-RNEMD moves,
-which were attempted every 50 fs. 
+which were attempted every 50~fs. 
 
 The interfaces were equilibrated for a total of 10 ns at equilibrium
 conditions before being exposed to either a shear or thermal gradient.
 This consisted of 5 ns under a constant temperature (NVT) integrator
 set to 225K followed by 5 ns under a microcanonical integrator.  Weak
 thermal gradients were allowed to develop using the VSS-RNEMD (NVE)
-integrator using a a small thermal flux ($-2.0\times 10^{-6}$
+integrator using a small thermal flux ($-2.0\times 10^{-6}$
 kcal/mol/\AA$^2$/fs) for a duration of 5 ns to allow the gradient to
 stabilize.  The resulting temperature gradient was $\approx$ 10K over
 the entire 100 \AA\  box length, which was sufficient to keep the
@@ -317,10 +351,10 @@ allowed to stabilize for 1 ns before data collection b
 
 Velocity gradients were then imposed using the VSS-RNEMD (NVE)
 integrator with a range of momentum fluxes.  These gradients were
-allowed to stabilize for 1 ns before data collection began. Once
-established, four successive 0.5 ns runs were performed for each shear
+allowed to stabilize for 1~ns before data collection began. Once
+established, four successive 0.5~ns runs were performed for each shear
 rate.  During these simulations, snapshots of the system were taken
-every 1 ps, and statistics on the structure and dynamics in each bin
+every 1~ps, and statistics on the structure and dynamics in each bin
 were accumulated throughout the simulations.
 
 \section{Results and discussion}
@@ -329,15 +363,15 @@ interface, Haymet {\it et al.}\cite{} have utilized st
 Any order parameter or time correlation function that changes as one
 crosses an interface from a bulk liquid to a solid can be used to
 measure the width of the interface.  In previous work on the ice/water
-interface, Haymet {\it et al.}\cite{} have utilized structural
+interface, Haymet {\it et al.}\cite{Bryk02} have utilized structural
 features (including the density) as well as dynamic properties
 (including the diffusion constant) to estimate the width of the
 interfaces for a number of facets of the ice crystals.  Because
 VSS-RNEMD imposes a lateral flow, parameters that depend on
 translational motion of the molecules (e.g. diffusion) may be
-artifically skewed by the RNEMD moves.  A structural parameter is not
+artificially skewed by the RNEMD moves.  A structural parameter is not
 influenced by the RNEMD perturbations to the same degree. Here, we
-have used the local tetraherdal order parameter as described by
+have used the local tetrahedral order parameter as described by
 Kumar\cite{Kumar09} and Errington\cite{Errington01} as our principal
 measure of the interfacial width.
 
@@ -363,7 +397,7 @@ solvation shell was set to 3.41 \AA\ .  The $q_{z}$ fu
 
 To estimate the interfacial width, the system was divided into 100
 bins along the $z$-dimension, and a cutoff radius for the first
-solvation shell was set to 3.41 \AA\ .  The $q_{z}$ function was
+solvation shell was set to 3.41~\AA\ .  The $q_{z}$ function was
 time-averaged to give yield a tetrahedrality profile of the
 system. The profile was then fit to a hyperbolic tangent that smoothly
 links the liquid and solid states,
@@ -375,12 +409,12 @@ interfaces, respectively.  The last term in equation \
 Here $q_{liq}$ and $q_{ice}$ are the local tetrahedral order parameter
 for the bulk liquid and ice domains, respectively, $w$ is the width of
 the interface.  $l$ and $r$ are the midpoints of the left and right
-interfaces, respectively.  The last term in equation \eqref{tet_fit}
+interfaces, respectively.  The last term in eq. \eqref{tet_fit}
 accounts for the influence that the weak thermal gradient has on the
 tetrahedrality profile in the liquid region.  To estimate the
-10\%-90\% widths commonly used in previous studies,\cite{} it is a
-simple matter to scale the widths obtained from the hyberbolic tangent
-fits to obtain $w_{10-90} = 2.9 w$.\cite{}
+10\%-90\% widths commonly used in previous studies,\cite{Bryk02} it is
+a simple matter to scale the widths obtained from the hyperbolic
+tangent fits to obtain $w_{10-90} = 2.1971 \times w$.\cite{Bryk02}
 
 In Figures \ref{fig:bComic} and \ref{fig:pComic} we see the
 $z$-coordinate profiles for tetrahedrality, temperature, and the
@@ -390,19 +424,19 @@ interfaces ($r$ and $l$ in equation \eqref{tet_fit}). 
 tetrahedral order parameter, $q(z) \approx 0.75$ while in the
 crystalline region, $q(z) \approx 0.94$, indicating a more tetrahedral
 environment.  The vertical dotted lines denote the midpoint of the
-interfaces ($r$ and $l$ in equation \eqref{tet_fit}). The weak thermal
+interfaces ($r$ and $l$ in eq. \eqref{tet_fit}). The weak thermal
 gradient applied to the systems in order to keep the interface at
-225$\pm$5K, can be seen in middle panels.  The tranverse velocity
+225~$\pm$~5K, can be seen in middle panels.  The transverse velocity
 profile is shown in the upper panels.  It is clear from the upper
 panels that water molecules in close proximity to the surface (i.e.
-within 10 \AA\ to 15 \AA\ of the interfaces) have transverse
+within 10~\AA\ to 15~\AA~of the interfaces) have transverse
 velocities quite close to the velocities within the ice block.  There
 is no velocity discontinuity at the interface, which indicates that
 the shearing of ice/water interfaces occurs in the ``stick'' or
 no-slip boundary conditions.
 
 \begin{figure}
-\includegraphics[width=\linewidth]{bComicStrip.pdf}
+\includegraphics[width=\linewidth]{bComicStrip}
 \caption{\label{fig:bComic} The basal interfaces.  Lower panel: the
   local tetrahedral order parameter, $q(z)$, (black circles) and the
   hyperbolic tangent fit (red line).  Middle panel: the imposed
@@ -414,14 +448,14 @@ no-slip boundary conditions.
 \end{figure}
 
 \begin{figure}
-\includegraphics[width=\linewidth]{pComicStrip.pdf}
+\includegraphics[width=\linewidth]{pComicStrip}
 \caption{\label{fig:pComic} The prismatic interfaces.  Panel
   descriptions match those in figure \ref{fig:bComic}}
 \end{figure}
 
-From the fits using equation \eqref{tet_fit}, we find the interfacial
-width for the basal and prismatic systems to be 3.2$\pm$0.4 \AA\ and
-3.6$\pm$0.2 \AA\ , respectively, with no applied momentum flux. Over
+From the fits using eq. \eqref{tet_fit}, we find the interfacial
+width for the basal and prismatic systems to be 3.2~$\pm$~0.4~\AA\ and
+3.6~$\pm$~0.2~\AA\ , respectively, with no applied momentum flux. Over
 the range of shear rates investigated, $0.6 \pm 0.3 \mathrm{ms}^{-1}
 \rightarrow 5.3 \pm 0.5 \mathrm{ms}^{-1}$ for the basal system and
 $0.9 \pm 0.2 \mathrm{ms}^{-1} \rightarrow 4.5 \pm 0.1
@@ -430,7 +464,7 @@ error bars.
 over all shear rates contained the values reported above within their
 error bars.
 
-\subsubsection{Orientational Time Correlation Function}
+\subsubsection{Orientational Dynamics}
 The orientational time correlation function,
 \begin{equation}\label{C(t)1}
   C_{2}(t)=\langle P_{2}(\mathbf{u}(0) \cdot \mathbf{u}(t)) \rangle,
@@ -453,12 +487,13 @@ then divided into 30 bins and $C_2(t)$ was evaluated f
 the 0.5 ns simulations was followed by a shorter 200 ps microcanonical
 (NVE) simulation in which the positions and orientations of every
 molecule in the system were recorded every 0.1 ps. The systems were
-then divided into 30 bins and $C_2(t)$ was evaluated for each bin.
+then divided into 30 bins along the $z$-axis and $C_2(t)$ was
+evaluated for each bin.
 
 In simulations of water at biological interfaces, Furse {\em et al.}
 fit $C_2(t)$ functions for water with triexponential
 functions,\cite{Furse08} where the three components of the decay
-correspond to a fast (<200 fs) reorientational piece driven by the
+correspond to a fast ($<$200 fs) reorientational piece driven by the
 restoring forces of existing hydrogen bonds, a middle (on the order of
 several ps) piece describing the large angle jumps that occur during
 the breaking and formation of new hydrogen bonds,and a slow (on the
@@ -473,7 +508,7 @@ time domains, as well as a constant piece that account
 temperatures, and the water molecules are further perturbed by the
 presence of the ice phase nearby.  We have obtained the most
 reasonable fits using triexponential functions with three distinct
-time domains, as well as a constant piece that accounts for the water
+time domains, as well as a constant piece to account for the water
 stuck in the ice phase that does not experience any long-time
 orientational decay,
 \begin{equation}
@@ -487,14 +522,14 @@ slab.
 slab.
 
 \begin{figure}
-\includegraphics[width=\linewidth]{basal_Tau_comic_strip.pdf}
+\includegraphics[width=\linewidth]{basal_Tau_comic_strip}
 \caption{\label{fig:basal_Tau_comic_strip} The three decay constants
   of the orientational time correlation function, $C_2(t)$, for water
   as a function of distance from the center of the ice slab.  The
   dashed line indicates the location of the basal face (as determined
   from the tetrahedrality order parameter).  The moderate and long
   time contributions slow down close to the interface which would be
-  expected under reorganizations that inolve large motions of the
+  expected under reorganizations that involve large motions of the
   molecules (e.g. frame-reorientations and jumps).  The observed
   speed-up in the short time contribution is surprising, but appears
   to reflect the restricted motion of librations closer to the
@@ -502,7 +537,7 @@ slab.
 \end{figure}
 
 \begin{figure}
-\includegraphics[width=\linewidth]{prismatic_Tau_comic_strip.pdf}
+\includegraphics[width=\linewidth]{prismatic_Tau_comic_strip}
 \caption{\label{fig:prismatic_Tau_comic_strip}
  Decay constants for $C_2(t)$ at the prismatic interface.  Panel
   descriptions match those in figure \ref{fig:basal_Tau_comic_strip}.}
@@ -532,83 +567,160 @@ One remarkable feature, is that for each of the interf
 \ref{fig:prismatic_Tau_comic_strip}, while those obtained while a
 shear was active are shown in red.
 
-One remarkable feature, is that for each of the interfaces, there is
-an apparent fixed liquid-state value for $\tau_{short}$,
+Two notable features deserve clarification.  First, there are
+nearly-constant liquid-state values for $\tau_{short}$,
 $\tau_{middle}$, and $\tau_{long}$ at large displacements from the
-interface. There also appears to be a single distance, $d_{basal}$ or
-$d_{prismatic}$, from the interface at which all three decay times
-begin to deviate from their bulk liquid values. We find $d_{basal}$
-and $d_{prismatic}$ to be roughly 15 \AA\ and 8 \AA\ respectively.
-These two results indicate that the dynamics of the water molecules
-within $d_{basal}$ and $d_{prismatic}$ are being significantly
-perturbed by the interface, even though the structural width of the
-interface via analysis of the tetrahedrality profile indicates that
-bulk liquid structure of water is recovered after about 4 \AA\ from
-the edge of the ice.
+interface. Second, there appears to be a single distance, $d_{basal}$
+or $d_{prismatic}$, from the interface at which all three decay times
+begin to deviate from their bulk liquid values. We find these
+distances to be approximately 15~\AA\ and 8~\AA\, respectively,
+although significantly finer binning of the $C_2(t)$ data would be
+necessary to provide better estimates of a ``dynamic'' interfacial
+thickness.
 
-Beaglehole and Wilson have measured the ice/water interface to have a
-thickness approximately 10 \AA\ for both the basal and prismatic face
-of ice by ellipticity measurements \cite{Beaglehole93}. Structurally,
-we have found the basal and prismatic interfacial width to be
-3.2$\pm$0.4 \AA\ and 3.6$\pm$0.2 \AA\ . However, decomposition of the
-spatial dependence of the decay times of $C_2(t)$ appears to indicate
-that a somewhat thicker interfacial region exists in which the
-orientational dynamics of the water molecules begin to resemble the
-trapped interfacial water more than the surrounding bulk.
+Beaglehole and Wilson have measured the ice/water interface using
+ellipsometry and find a thickness of approximately 10~\AA\ for both
+the basal and prismatic faces.\cite{Beaglehole93} Structurally, we
+have found the basal and prismatic interfacial width to be
+3.2~$\pm$~0.4~\AA\ and 3.6~$\pm$~0.2~\AA. However, decomposition of
+the spatial dependence of the decay times of $C_2(t)$ indicates that a
+somewhat thicker interfacial region exists in which the orientational
+dynamics of the water molecules begin to resemble the trapped
+interfacial water more than the surrounding liquid.
 
+Our results indicate that the dynamics of the water molecules within
+$d_{basal}$ and $d_{prismatic}$ are being significantly perturbed by
+the interface, even though the structural width of the interface via
+analysis of the tetrahedrality profile indicates that bulk liquid
+structure of water is recovered after about 4 \AA\ from the edge of
+the ice.
+
 \subsection{Coefficient of Friction of the Interface}
-As the ice is sheared through the liquid, there will be a friction
-between the solid and the liquid. Pit has shown how to calculate the
-coefficient of friction $\lambda$ for a solid-liquid interface for a
-Newtonian fluid of viscosity $\eta$ and has a slip length of
-$\delta$. \cite{Pit99}
+As liquid water flows over an ice interface, there is a distance from
+the structural interface where bulk-like hydrodynamics are recovered.
+Bocquet and Barrat constructed a theory for the hydrodynamic boundary
+parameters, which include the slipping length
+$\left(\delta_\mathrm{wall}\right)$ of this boundary layer and the
+``hydrodynamic position'' of the boundary
+$\left(z_\mathrm{wall}\right)$.\cite{PhysRevLett.70.2726,PhysRevE.49.3079}
+This last parameter is the location (relative to a solid surface)
+where the bulk-like behavior is recovered.  Work by Mundy {\it et al.}
+has helped to combine these parameters into a liquid-solid friction
+coefficient, which quantifies the resistance to pulling the solid
+interface through a liquid,\cite{Mundy1997305}
+\begin{equation}
+\lambda_\mathrm{wall} = \frac{\eta}{\delta_\mathrm{wall}}.
+\end{equation}
+This expression is nearly identical to one provided by Pit {\it et
+  al.} for the solid-liquid friction of an interface,\cite{Pit99}
 \begin{equation}\label{Pit}
-\lambda=\eta/\delta
+  \lambda=\frac{\eta}{\delta}
 \end{equation}
-From linear response theory, $\eta$ can be obtained from the imposed
-momentum flux and the slope of the velocity about the dimension of the
-imposed flux.\cite{Kuang12}
+where $\delta$ is the slip length for the liquid measured at the
+location of the interface itself.
+
+In both of these expressions, $\eta$ is the shear viscosity of the
+bulk-like region of the liquid, a quantity which is easily obtained in
+VSS-RNEMD simulations by fitting the velocity profile in the region
+far from the surface.\cite{Kuang12} Assuming linear response in the
+bulk-like region,
 \begin{equation}\label{Kuang}
-j_{z}(p_{x})=-\eta\frac{\partial v_{x}}{\partial z}
+j_{z}(p_{x})=-\eta \left(\frac{\partial v_{x}}{\partial z}\right)
 \end{equation}
-Solving eq. \eqref{Kuang} for $\eta$ and substituting the result into
-eq. \eqref{Pit}, we obtain an alternate expression for the coefficient
-of friction.
+Substituting this result into eq. \eqref{Pit}, we can estimate the
+solid-liquid coefficient using the slip length,
 \begin{equation}
-\lambda=-\frac{j_{z}(p_{x})}{\delta \frac{\partial v_{x}}{\partial z}}
+\lambda=-\frac{j_{z}(p_{x})} {\left(\frac{\partial v_{x}}{\partial
+      z}\right) \delta}
 \end{equation}
 
-For our simulations, we obtain $\delta$ from the difference between the structural edge of the ice block determined by the tetrahedrality profile fit, and the intersection of the linear regression of the $v_{x}$ profiles about the $z$-dimension for the ice and liquid. (See Figure \ref{fig:delta_example}) The coefficient of friction for the basal and the prismatic facets were found to be (0.07$\pm$0.01) amu fs\textsuperscript{-1} and (0.032$\pm$0.007) amu fs\textsuperscript{-1}. It is known that the basal and prismatic faces have different surface structures. The basal face is smoother than the prismatic with small alternating valleys and crests, while the prismatic surface has deep corrugating channels. We believe the reason that the prismatic face's coefficient of friction was found to be smaller than the basal's is due to the direction of the shear. The shear of the ice/water was in the same direction of the corrugating channels, allowing water molecules to pass through the channels during the shear. 
+For ice / water interfaces, the boundary conditions are markedly
+no-slip, so projecting the bulk liquid state velocity profile yields a
+negative slip length. This length is the difference between the
+structural edge of the ice (determined by the tetrahedrality profile)
+and the location where the projected velocity of the bulk liquid
+intersects the solid phase velocity (see Figure
+\ref{fig:delta_example}). The coefficients of friction for the basal
+and the prismatic facets are found to be
+0.07~$\pm$~0.01~amu~fs\textsuperscript{-1} and
+0.032~$\pm$~0.007~amu~fs\textsuperscript{-1}, respectively. These
+results may seem surprising as the basal face is smoother than the
+prismatic with only small undulations of the oxygen positions, while
+the prismatic surface has deep corrugated channels. The applied
+momentum flux used in our simulations is parallel to these channels,
+however, and this results in a flow of water in the same direction as
+the corrugations, allowing water molecules to pass through the
+channels during the shear.
 
  \begin{figure}
-\includegraphics[width=\linewidth]{delta_example.pdf}
-\caption{\label{fig:delta_example} A schematic of determining the slip length ($\delta$). The slip length is the difference of the structural starting point of the ice and the point of intersection of the linear regressions of the liquid phase velocity profile (red) and of the solid ice velocity profile (black). The dotted line indicates the location of the ice as determined by the tetrahedrality profile.}
+\includegraphics[width=\linewidth]{delta_example}
+\caption{\label{fig:delta_example} Determining the (negative) slip
+  length ($\delta$) for the ice-water interfaces (which have decidedly
+  non-slip behavior).  This length is the difference between the
+  structural edge of the ice (determined by the tetrahedrality
+  profile) and the location where the projected velocity of the bulk
+  liquid (dashed red line) intersects the solid phase velocity (solid
+  black line).  The dotted line indicates the location of the ice as
+  determined by the tetrahedrality profile.}
 \end{figure}
 
 
 \section{Conclusion}
-Here we have simulated the basal and prismatic facets of an SPC/E model of the ice I$_\mathrm{h}$ / water interface. Using VSS-RNEMD, the ice was sheared relative to the liquid while imposed thermal gradients kept the interface at a stable temperature. Caculation of the local tetrahedrality order parameter has shown an apparent independence of the shear rate on the interfacial width, which was found to be 3.2$\pm$0.4 \AA\ and 3.6$\pm$0.2 \AA\ for the basal and prismatic systems. The orientational time correlation function was calculated from varying displacements from the interface. Decomposition by a triexponential decay also showed an apparent independence of the shear rate. The short time decay due to the restoring forces of existing hydrogen bonds decreased at close displacements from the interface, while the middle and long time decays were found to increase. There is also an apparent displacement, $d_{basal}$ and $d_{prismatic}$, from the interface at which these deviations from bulk liquid values occurs. We found $d_{basal}$ and $d_{prismatic}$ to be approximately 15 \AA\ and 8 \AA\ . This implies that the dynamics of water molecules which are structurally equivalent to bulk phase molecules are being perturbed by the presence of the ice and/or the interface. The coefficient of friction of each of the facets was also determined. They were found to be (0.07$\pm$0.01) amu fs\textsuperscript{-1} and (0.032$\pm$0.007) amu fs\textsuperscript{-1} for the basal and prismatic facets respectively. We believe the large difference between the two friction coefficients is due to the direction of the shear and the surface structure of the crystal facets. 
+We have simulated the basal and prismatic facets of an SPC/E model of
+the ice I$_\mathrm{h}$ / water interface.  Using VSS-RNEMD, the ice
+was sheared relative to the liquid while simultaneously being exposed
+to a weak thermal gradient which kept the interface at a stable
+temperature.  Calculation of the local tetrahedrality order parameter
+has shown an apparent independence of the interfacial width on the
+shear rate.  This width was found to be 3.2~$\pm$0.4~\AA\ and
+3.6~$\pm$0.2~\AA\ for the basal and prismatic systems, respectively.
 
-\begin{acknowledgement}
+Orientational time correlation functions were calculated at varying
+displacements from the interface, and were found to be similarly
+independent of the applied momentum flux. The short decay due to the
+restoring forces of existing hydrogen bonds decreased close to the
+interface, while the longer-time decay constants increased in close
+proximity to the interface.  There is also an apparent dynamic
+interface width, $d_{basal}$ and $d_{prismatic}$, at which these
+deviations from bulk liquid values begin.  We found $d_{basal}$ and
+$d_{prismatic}$ to be approximately 15~\AA\ and 8~\AA\ . This implies
+that the dynamics of water molecules which have similar structural
+environments to liquid phase molecules are dynamically perturbed by
+the presence of the ice interface.
+
+The coefficient of liquid-solid friction for each of the facets was
+also determined. They were found to be
+0.07~$\pm$~0.01~amu~fs\textsuperscript{-1} and
+0.032~$\pm$~0.007~amu~fs\textsuperscript{-1} for the basal and
+prismatic facets respectively. We attribute the large difference
+between the two friction coefficients to the direction of the shear
+and to the surface structure of the crystal facets.
+
+\section{Acknowledgements}
   Support for this project was provided by the National Science
   Foundation under grant CHE-0848243. Computational time was provided
   by the Center for Research Computing (CRC) at the University of
   Notre Dame.
-\end{acknowledgement}
 
 \newpage
-\bibstyle{achemso}
 \bibliography{iceWater}
 
-\begin{tocentry}
-\begin{wrapfigure}{l}{0.5\textwidth}
-\begin{center}
-\includegraphics[width=\linewidth]{SystemImage.png}
-\end{center}
-\end{wrapfigure}
-An image of our system.
-\end{tocentry}
+\end{doublespace}
 
+% \begin{tocentry}
+% \begin{wrapfigure}{l}{0.5\textwidth}
+% \begin{center}
+% \includegraphics[width=\linewidth]{SystemImage.png}
+% \end{center}
+% \end{wrapfigure}
+% The cell used to simulate liquid-solid shear in ice I$_\mathrm{h}$ /
+% water interfaces.  Velocity gradients were applied using the velocity
+% shearing and scaling variant of reverse non-equilibrium molecular
+% dynamics (VSS-RNEMD) with a weak thermal gradient to prevent melting.
+% The interface width is relatively robust in both structual and dynamic
+% measures as a function of the applied shear.
+% \end{tocentry}
+
 \end{document}
 
 %**************************************************************
