| 251 |
|
\caption{The external temperature and pressure bath interacts only |
| 252 |
|
with those atoms on the convex hull (grey surface). The hull is |
| 253 |
|
computed dynamically at each time step, and molecules can move |
| 254 |
< |
between the interior (Newtonian) region and the Langevin hull.} |
| 254 |
> |
between the interior (Newtonian) region and the Langevin Hull.} |
| 255 |
|
\label{fig:hullSample} |
| 256 |
|
\end{figure} |
| 257 |
|
|
| 646 |
|
\includegraphics[width=\linewidth]{pAngle} |
| 647 |
|
\caption{Distribution of $\cos{\theta}$ values for molecules on the |
| 648 |
|
interior of the cluster (squares) and for those participating in the |
| 649 |
< |
convex hull (circles) at a variety of pressures. The Langevin hull |
| 649 |
> |
convex hull (circles) at a variety of pressures. The Langevin Hull |
| 650 |
|
exhibits minor dewetting behavior with exposed oxygen sites on the |
| 651 |
|
hull water molecules. The orientational preference for exposed |
| 652 |
|
oxygen appears to be independent of applied pressure. } |
| 661 |
|
forming a dangling hydrogen bond acceptor site. |
| 662 |
|
|
| 663 |
|
Previous molecular dynamics simulations |
| 664 |
< |
of SPC/E water that use periodic boundary conditions have shown that molecules on the liquid side of the liquid/vapor interface favor a similar orientation where oxygen is directed away from the bulk.\cite{Taylor1996} These simulations had both a liquid phase and a well-defined vapor phase in equilibrium and showed that vapor molecules generally had one hydrogen protruding from the surface, forming a dangling hydrogen bond donor. Our water cluster simulations do not have a true lasting vapor phase, but rather a few transient molecules that leave the liquid droplet. Thus we are unable to comment on the orientational preference of vapor phase molecules in a Langevin Hull simulation. |
| 664 |
> |
of SPC/E water using periodic boundary conditions have shown that molecules on the liquid side of the liquid/vapor interface favor a similar orientation where oxygen is directed away from the bulk.\cite{Taylor1996} These simulations had both a liquid phase and a well-defined vapor phase in equilibrium and showed that vapor molecules generally had one hydrogen protruding from the surface, forming a dangling hydrogen bond donor. Our water cluster simulations do not have a true lasting vapor phase, but rather a few transient molecules that leave the liquid droplet. Thus we are unable to comment on the orientational preference of vapor phase molecules in a Langevin Hull simulation. |
| 665 |
|
|
| 666 |
< |
However, the orientational preference exhibited by liquid phase hull molecules in the Langevin hull is significantly weaker than the preference caused by an explicit hydrophobic bounding potential. Additionally, the Langevin Hull does not require that the orientation of any molecules be fixed in order to maintain bulk-like structure, even at the cluster surface. |
| 666 |
> |
However, the orientational preference exhibited by liquid phase hull molecules in the Langevin Hull is significantly weaker than the preference caused by an explicit hydrophobic bounding potential. Additionally, the Langevin Hull does not require that the orientation of any molecules be fixed in order to maintain bulk-like structure, even at the cluster surface. |
| 667 |
|
|
| 668 |
|
\subsection{Heterogeneous nanoparticle / water mixtures} |
| 669 |
|
|
| 670 |
|
To further test the method, we simulated gold nanopartices ($r = 18$ |
| 671 |
|
\AA) solvated by explicit SPC/E water clusters using the Langevin |
| 672 |
< |
hull. This was done at pressures of 1, 2, 5, 10, 20, 50, 100 and 200 atm |
| 672 |
> |
Hull. This was done at pressures of 1, 2, 5, 10, 20, 50, 100 and 200 atm |
| 673 |
|
in order to observe the effects of pressure on the ordering of water |
| 674 |
|
ordering at the surface. In Fig. \ref{fig:RhoR} we show the density |
| 675 |
|
of water adjacent to the surface as a function of pressure, as well as |
| 766 |
|
|
| 767 |
|
For a large number of atoms on a moderately parallel machine, the |
| 768 |
|
total costs are dominated by the computations of the individual hulls, |
| 769 |
< |
and communication of these hulls to create the Langevin hull sees roughly |
| 769 |
> |
and communication of these hulls to create the Langevin Hull sees roughly |
| 770 |
|
linear speed-up with increasing processor counts. |
| 771 |
|
|
| 772 |
|
\section*{Acknowledgments} |
