50 |
|
\date{\today} |
51 |
|
|
52 |
|
\begin{abstract} |
53 |
< |
Abstract abstract abstract... |
53 |
> |
|
54 |
> |
|
55 |
|
\end{abstract} |
56 |
|
|
57 |
|
\maketitle |
203 |
|
loses memory of its former orientation. Observing the rate at which this decay |
204 |
|
occurs can provide insight to the mechanism and timescales for the relaxation. |
205 |
|
In eq. \eqref{C(t)1}, $P_{2}$ is the second-order Legendre polynomial, and |
206 |
< |
$\mathbf{u}$ is the the bisecting HOH vector. The angle brackets indicate |
206 |
> |
$\mathbf{u}$ is the bisecting HOH vector. The angle brackets indicate |
207 |
|
an ensemble average over all the water molecules in a given spatial region. |
208 |
|
|
209 |
|
To investigate the dynamics of the water molecules across the interface, the |
273 |
|
dependence |
274 |
|
arrises from the sharp discontinuity of the viscosity for the SPC/E model |
275 |
|
at temperatures approaching 200 K\cite{kuang12}. Due to this, we propose |
276 |
< |
a weighting to the structural interfacial parameter, $\kappa$ by the |
277 |
< |
viscosity at $225$ K, the temperature of the interface. $\kappa$ is |
278 |
< |
traditionally defined as |
279 |
< |
\begin{equation}\label{kappa} |
280 |
< |
\kappa = \eta/\delta |
281 |
< |
\end{equation} |
281 |
< |
where $\eta$ is the viscosity and $\delta$ is the slip length. |
282 |
< |
In our ice/water shearing simulations, the system has reached a steady state |
283 |
< |
when the applied force, |
284 |
< |
|
285 |
< |
\begin{equation} |
286 |
< |
f_{applied} = \mathbf{j}_z(\mathbf{p})L_x L_y |
287 |
< |
\end{equation} |
288 |
< |
is equal to the frictional force resisting the motion of the ice block |
289 |
< |
\begin{equation} |
290 |
< |
f_{friction} = \frac{\eta \mathbf{v} L_x L_y}{\delta} |
276 |
> |
a weighting to the interfacial friction coefficient, $\kappa$ by the |
277 |
> |
shear viscosity at 225 K. The interfacial friction coefficient relates |
278 |
> |
the shear stress with the relative velocity of the fluid normal to the |
279 |
> |
interface: |
280 |
> |
\begin{equation}\label{Shenyu-13} |
281 |
> |
j_{z}(p_{x}) = \kappa[v_{x}(fluid)-v_{x}(solid)] |
282 |
|
\end{equation} |
283 |
< |
where $\mathbf{v}$ is the relative velocity of the liquid from the ice. |
284 |
< |
When this condition is met, we are able to solve the resulting expression to |
285 |
< |
obtain, |
286 |
< |
\begin{equation}\label{force_equality} |
287 |
< |
\frac{\eta}{\delta} = \frac{\mathbf{j}_z(\mathbf{p})}{\mathbf{v}} |
283 |
> |
where $j_{z}(p_{x})$ is the applied momentum flux (shear stress) across $z$ |
284 |
> |
in the |
285 |
> |
$x$-dimension, and $v_{x}$(fluid) and $v_{x}$(solid) are the velocities |
286 |
> |
directly adjacent to the interface. The shear viscosity, $\eta(T)$, of the |
287 |
> |
fluid can be determined if we assume a linear response of the momentum |
288 |
> |
gradient to the applied shear stress by |
289 |
> |
\begin{equation}\label{Shenyu-11} |
290 |
> |
j_{z}(p_{x}) = \eta(T) \frac{\partial v_{x}}{\partial z}. |
291 |
> |
\end{equation} |
292 |
> |
Using eqs \eqref{Shenyu-13} and \eqref{Shenyu-11}, we can find the following |
293 |
> |
expression for $\kappa$, |
294 |
> |
\begin{equation}\label{kappa-1} |
295 |
> |
\kappa = \eta(T) \frac{\partial v_{x}}{\partial z}\frac{1}{[v_{x}(fluid)-v_{x}(solid)]}. |
296 |
|
\end{equation} |
297 |
< |
From \eqref{kappa}, \eqref{force_equality} becomes |
298 |
< |
\begin{equation} |
299 |
< |
\kappa = \frac{\mathbf{j}_z(\mathbf{p})}{\mathbf{v}} |
297 |
> |
Here is where we will introduce the weighting term of $\eta(225)/\eta(T)$ |
298 |
> |
giving us |
299 |
> |
\begin{equation}\label{kappa-2} |
300 |
> |
\kappa = \frac{\eta(225)}{[v_{x}(fluid)-v_{x}(solid)]}\frac{\partial v_{x}}{\partial z}. |
301 |
|
\end{equation} |
302 |
– |
which we will multiply by a viscosity weighting term to reach |
303 |
– |
\begin{equation} \label{kappa2} |
304 |
– |
\kappa = \frac{\mathbf{j}_z(\mathbf{p})}{\mathbf{v}} \frac{\eta(225)}{\eta(T)} |
305 |
– |
\end{equation} |
306 |
– |
Assuming linear response theory is valid, an expression for ($\eta$) can |
307 |
– |
be found from the imposed momentum flux and the measured velocity gradient. |
308 |
– |
\begin{equation}\label{eta_eq} |
309 |
– |
\eta = \frac{\mathbf{j}_z(\mathbf{p})}{\frac{\partial v_x}{\partial z}} |
310 |
– |
\end{equation} |
311 |
– |
Substituting eq \eqref{eta_eq} into eq \eqref{kappa2} we arrive at |
312 |
– |
\begin{equation} |
313 |
– |
\kappa = \frac{\frac{\partial v_x}{\partial z}}{\mathbf{v}}\eta(225) |
314 |
– |
\end{equation} |
302 |
|
|
316 |
– |
\begin{table}[h] |
317 |
– |
\centering |
318 |
– |
\caption{$\kappa$ values for the basal, prismatic, pyramidal, and secondary prismatic facets of Ice-I$_\mathrm{h}$} |
319 |
– |
\label{tab:kappa} |
320 |
– |
\begin{tabular}{|ccc|} \hline |
321 |
– |
& \multicolumn{2}{c|}{$\kappa_{Drag direction}$} \\ |
322 |
– |
Interface & $\kappa_{x}$ & $\kappa_{y}$ \\ \hline |
323 |
– |
basal & $0.00059 \pm 0.00003$ & $0.00065 \pm 0.00008$ \\ |
324 |
– |
prismatic & $0.00030 \pm 0.00002$ & $0.00030 \pm 0.00001$ \\ |
325 |
– |
pyramidal & $0.00058 \pm 0.00004$ & $0.00061 \pm 0.00005$ \\ |
326 |
– |
secondary prism & $0.00035 \pm 0.00001$ & $0.00033 \pm 0.00002$ \\ \hline |
327 |
– |
\end{tabular} |
328 |
– |
\end{table} |
329 |
– |
|
330 |
– |
|
303 |
|
To obtain the value of $\eta(225)$ for the SPC/E model, a $31.09 \times 29.38 |
304 |
|
\times 124.39$ \AA\ box with 3744 SPC/E water molecules was equilibrated to |
305 |
|
225K, |
310 |
|
roughly ten times larger than the value found for 280 K SPC/E bulk water by |
311 |
|
Kuang\cite{kuang12}. |
312 |
|
|
313 |
< |
The resulting $\kappa$ values found for the four crystal |
314 |
< |
facets of Ice-I$_\mathrm{h}$ investigated are shown in Table \ref{tab:kappa}. |
315 |
< |
The basal and pyramidal facets were found to have similar values of |
316 |
< |
$\kappa \approx$ 0.0006, while $\kappa \approx$ 0.0003 were found for the |
317 |
< |
prismatic and secondary prismatic facets. |
313 |
> |
The interfacial friction coefficient, $\kappa$, can equivalently be expressed |
314 |
> |
as the ratio of the viscosity of the fluid to the slip length, $\delta$, which |
315 |
> |
is an indication of how 'slippery' the interface is. |
316 |
> |
\begin{equation}\label{kappa-3} |
317 |
> |
\kappa = \frac{\eta}{\delta} |
318 |
> |
\end{equation} |
319 |
> |
In each of the systems, the interfacial temperature was kept fixed to 225K, |
320 |
> |
which ensured the viscosity of the fluid at the |
321 |
> |
interace was approximately the same. Thus, any significant variation in |
322 |
> |
$\kappa$ between |
323 |
> |
the systems indicates differences in the 'slipperiness' of the interfaces. |
324 |
> |
As each of the ice systems are sheared relative to liquid water, the |
325 |
> |
'slipperiness' of the interface can be taken as an indication of how |
326 |
> |
hydrophobic or hydrophilic the interface is. The calculated $\kappa$ values |
327 |
> |
found for the four crystal facets of Ice-I$_\mathrm{h}$ investigated are shown |
328 |
> |
in Table \ref{tab:kapa}. The basal and pyramidal facets were found to have |
329 |
> |
similar values of $\kappa \approx$ 0.0006 (units), while $\kappa \approx$ |
330 |
> |
0.0003 (units) were found for the prismatic and secondary prismatic systems. |
331 |
> |
These results indicate that the prismatic and secondary prismatic facets are |
332 |
> |
more hydrophobic than the basal and pyramidal facets. |
333 |
|
%This indicates something about the similarity between the two facets that |
334 |
|
%share similar values... |
335 |
|
%Maybe find values for kappa for other materials to compare against? |
336 |
|
|
337 |
+ |
\begin{table}[h] |
338 |
+ |
\centering |
339 |
+ |
\caption{$\kappa$ values for the basal, prismatic, pyramidal, and secondary \ |
340 |
+ |
prismatic facets of Ice-I$_\mathrm{h}$} |
341 |
+ |
\label{tab:kappa} |
342 |
+ |
\begin{tabular}{|ccc|} \hline |
343 |
+ |
& \multicolumn{2}{c|}{$\kappa_{Drag direction}$ (units)} \\ |
344 |
+ |
Interface & $\kappa_{x}$ & $\kappa_{y}$ \\ \hline |
345 |
+ |
basal & $0.00059 \pm 0.00003$ & $0.00065 \pm 0.00008$ \\ |
346 |
+ |
prismatic & $0.00030 \pm 0.00002$ & $0.00030 \pm 0.00001$ \\ |
347 |
+ |
pyramidal & $0.00058 \pm 0.00004$ & $0.00061 \pm 0.00005$ \\ |
348 |
+ |
secondary prism & $0.00035 \pm 0.00001$ & $0.00033 \pm 0.00002$ \\ \hline |
349 |
+ |
\end{tabular} |
350 |
+ |
\end{table} |
351 |
|
|
352 |
|
|
353 |
+ |
|
354 |
+ |
|
355 |
|
%\begin{table}[h] |
356 |
|
%\centering |
357 |
|
%\caption{Solid-liquid friction coefficients (measured in amu~fs\textsuperscript\ |
402 |
|
interface by using the minimal perturbing VSS RNEMD method. In agreement with |
403 |
|
our previous findings for the basal and prismatic facets, the interfacial |
404 |
|
width was found to be independent of shear rate as measured by the local |
405 |
< |
order tetrahedral ordering parameter. This width was found to be |
406 |
< |
3.2~$\pm$0.2~\AA\ for both the pyramidal and the secondary prismatic systems. |
405 |
> |
tetrahedral order parameter. This width was found to be |
406 |
> |
3.2~$\pm$0.2~\AA\ for both the pyramidal and the secondary prismatic systems. |
407 |
> |
These values are in good agreement with our previously calculated interfacial |
408 |
> |
widths for the basal (3.2~$\pm$0.4~\AA\ ) and prismatic (3.6~$\pm$0.2~\AA\ ) |
409 |
> |
systems. |
410 |
|
|
411 |
+ |
Orientational dynamics of the Ice-I$_\mathrm{h}$/water interfaces were studied |
412 |
+ |
by calculation of the orientational time correlation function at varying |
413 |
+ |
displacements normal to the interface. The decays were fit |
414 |
+ |
to a tri-exponential decay, where the three decay constants correspond to |
415 |
+ |
the librational motion of the molecules driven by the restoring forces of |
416 |
+ |
existing hydrogen bonds ($\tau_{short}$ $\mathcal{O}$(10 fs)), jumps between |
417 |
+ |
two different hydrogen bonds ($\tau_{middle}$ $\mathcal{O}$(1 ps)), and |
418 |
+ |
translational motion of the molecules ($\tau_{long}$ $\mathcal{O}$(100 ps)). |
419 |
+ |
$\tau_{short}$ was found to decrease approaching the interface due to the |
420 |
+ |
constrained motion of the molecules as the local environment becomes more |
421 |
+ |
ice-like. Conversely, the two longer-time decay constants were found to |
422 |
+ |
increase at small displacements from the interface. As seen in our previous |
423 |
+ |
work on the basal and prismatic facets, there appears to be a dynamic |
424 |
+ |
interface width at which deviations from the bulk liquid values occur. |
425 |
+ |
We had previously found $d_{basal}$ and $d_{prismatic}$ to be approximately |
426 |
+ |
2.8~\AA\ and 3.5~\AA\~. We found good agreement of these values for the |
427 |
+ |
pyramidal and secondary prismatic systems with $d_{pyramidal}$ and |
428 |
+ |
$d_{secondary prism}$ to be 2.7~\AA\ and 3~\AA\~. For all of the facets, there |
429 |
+ |
was found to be no apparent dependence of the dynamic width on the shear rate. |
430 |
+ |
|
431 |
+ |
%Paragraph summarizing the Kappa values |
432 |
+ |
The interfacial friction coefficient, $\kappa$, was determined for each of the |
433 |
+ |
interfaces. We were able to reach an expression for $\kappa$ as a function of |
434 |
+ |
the velocity profile of the system and is scaled by the viscosity of the liquid |
435 |
+ |
at 225 K. In doing so, we have obtained an expression for $\kappa$ which is |
436 |
+ |
independent of temperature differences of the liquid water at far displacements |
437 |
+ |
from the interface. We found the basal and pyramidal facets to have |
438 |
+ |
similar $\kappa$ values, of $\kappa \approx$ 0.0006 (units). However, the |
439 |
+ |
prismatic and secondary prismatic facets were found to have $\kappa$ values of |
440 |
+ |
$\kappa \approx$ 0.0003 (units). As these ice facets are being sheared |
441 |
+ |
relative to liquid water, with the structural and dynamic width of all four |
442 |
+ |
interfaces being approximately the same, the difference in the coefficient of |
443 |
+ |
friction indicates the hydrophilicity of the crystal facets are not |
444 |
+ |
equivalent. Namely, that the basal and pyramidal facets of Ice-I$_\mathrm{h}$ |
445 |
+ |
are more hydrophilic than the prismatic and secondary prismatic facets. |
446 |
|
|
447 |
|
|
407 |
– |
|
448 |
|
\begin{acknowledgments} |
449 |
|
Support for this project was provided by the National |
450 |
|
Science Foundation under grant CHE-1362211. Computational time was |