310 |
|
tensor that includes translational and rotational drag as well as |
311 |
|
translational-rotational coupling. The computation of resistance |
312 |
|
tensors for rigid bodies has been detailed |
313 |
< |
elsewhere,\cite{JoseGarciadelaTorre02012000,Garcia-de-la-Torre:2001wd,GarciadelaTorreJ2002,Sun2008} |
313 |
> |
elsewhere,\cite{JoseGarciadelaTorre02012000,Garcia-de-la-Torre:2001wd,GarciadelaTorreJ2002,Sun:2008fk} |
314 |
|
but the standard approach involving bead approximations would be |
315 |
|
prohibitively expensive if it were recomputed at each step in a |
316 |
|
molecular dynamics simulation. |
353 |
|
We have implemented this method by extending the Langevin dynamics |
354 |
|
integrator in our code, OpenMD.\cite{Meineke2005,openmd} The Delaunay |
355 |
|
triangulation and computation of the convex hull are done using calls |
356 |
< |
to the qhull library.\cite{qhull} There is a moderate penalty for |
356 |
> |
to the qhull library.\cite{Qhull} There is a moderate penalty for |
357 |
|
computing the convex hull at each step in the molecular dynamics |
358 |
|
simulation (HOW MUCH?), but the convex hull is remarkably easy to |
359 |
|
parallelize on distributed memory machines (see Appendix A). |
361 |
|
\section{Tests \& Applications} |
362 |
|
\label{sec:tests} |
363 |
|
|
364 |
< |
In order to test this method, we have carried out simulations using |
365 |
< |
the Langevin Hull on a crystalline system (gold nanoparticles), a |
366 |
< |
liquid droplet (SPC/E water), and a heterogeneous mixture (gold |
367 |
< |
nanoparticles in a water droplet). In each case, we have computed |
368 |
< |
properties that depend on the external applied pressure. Of |
364 |
> |
To test the new method, we have carried out simulations using the |
365 |
> |
Langevin Hull on: 1) a crystalline system (gold nanoparticles), 2) a |
366 |
> |
liquid droplet (SPC/E water),\cite{SPCE} and 3) a heterogeneous |
367 |
> |
mixture (gold nanoparticles in a water droplet). In each case, we have |
368 |
> |
computed properties that depend on the external applied pressure. Of |
369 |
|
particular interest for the single-phase systems is the bulk modulus, |
370 |
|
\begin{equation} |
371 |
|
\kappa_{T} = -\frac{1}{V} \left ( \frac{\partial V}{\partial P} \right |
385 |
|
)_{T} |
386 |
|
\label{eq:BMN} |
387 |
|
\end{equation} |
388 |
< |
The region we pick is a spherical volume of 10 \AA radius centered in |
389 |
< |
the middle of the cluster. The geometry and size of the region is |
390 |
< |
arbitrary, and any bulk-like portion of the cluster can be used to |
391 |
< |
compute the bulk modulus. |
388 |
> |
The region we used is a spherical volume of 10 \AA\ radius centered in |
389 |
> |
the middle of the cluster. $N$ is the average number of molecules |
390 |
> |
found within this region throughout a given simulation. The geometry |
391 |
> |
and size of the region is arbitrary, and any bulk-like portion of the |
392 |
> |
cluster can be used to compute the bulk modulus. |
393 |
|
|
394 |
< |
One might also assume that the volume of the convex hull could be |
395 |
< |
taken as the system volume in the compressibility expression |
396 |
< |
(Eq. \ref{eq:BM}), but this has implications at lower pressures (which |
397 |
< |
are explored in detail in the section on water droplets). |
394 |
> |
One might assume that the volume of the convex hull could be taken as |
395 |
> |
the system volume in the compressibility expression (Eq. \ref{eq:BM}), |
396 |
> |
but this has implications at lower pressures (which are explored in |
397 |
> |
detail in the section on water droplets). |
398 |
|
|
399 |
+ |
The metallic force field in use for the gold nanoparticles is the |
400 |
+ |
quantum Sutton-Chen (QSC) model.\cite{PhysRevB.59.3527} In all |
401 |
+ |
simulations involving point charges, we utilized damped shifted-force |
402 |
+ |
(DSF) electrostatics\cite{Fennell06} which is a variant of the Wolf |
403 |
+ |
summation\cite{wolf:8254} that has been shown to provide good forces |
404 |
+ |
and torques on molecular models for water in a computationally |
405 |
+ |
efficient manner.\cite{Fennell06} The damping parameter ($\alpha$) was |
406 |
+ |
set to 0.18 \AA$^{-1}$, and the cutoff radius was set to 12 \AA. The |
407 |
+ |
Spohr potential was adopted in depicting the interaction between metal |
408 |
+ |
atoms and the SPC/E water molecules.\cite{ISI:000167766600035} |
409 |
+ |
|
410 |
|
\subsection{Bulk modulus of gold nanoparticles} |
411 |
|
|
412 |
+ |
The bulk modulus is well-known for gold, and it provides a good first |
413 |
+ |
test of how the method compares to other similar methods. |
414 |
+ |
|
415 |
+ |
|
416 |
|
\begin{figure} |
417 |
|
\includegraphics[width=\linewidth]{pressure_tb} |
418 |
|
\caption{Pressure response is rapid (18 \AA gold nanoparticle), target |
541 |
|
|
542 |
|
The orientational preference exhibited by hull molecules 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. |
543 |
|
|
528 |
– |
|
544 |
|
\subsection{Heterogeneous nanoparticle / water mixtures} |
545 |
|
|
546 |
|
|
547 |
< |
\section{Appendix A: Hydrodynamic tensor for triangular facets} |
547 |
> |
\section*{Appendix A: Computing Convex Hulls on Parallel Computers} |
548 |
|
|
549 |
< |
\section{Appendix B: Computing Convex Hulls on Parallel Computers} |
535 |
< |
|
536 |
< |
\section{Acknowledgments} |
549 |
> |
\section*{Acknowledgments} |
550 |
|
Support for this project was provided by the |
551 |
|
National Science Foundation under grant CHE-0848243. Computational |
552 |
|
time was provided by the Center for Research Computing (CRC) at the |