| 42 |
|
hull surrounding the system. A Langevin thermostat is also applied |
| 43 |
|
to the facets to mimic contact with an external heat bath. This new |
| 44 |
|
method, the ``Langevin Hull'', can handle heterogeneous mixtures of |
| 45 |
< |
materials with different compressibilities. These are systems that |
| 46 |
< |
are problematic for traditional affine transform methods. The |
| 47 |
< |
Langevin Hull does not suffer from the edge effects of boundary |
| 48 |
< |
potential methods, and allows realistic treatment of both external |
| 49 |
< |
pressure and thermal conductivity due to the presence of an implicit |
| 50 |
< |
solvent. We apply this method to several different systems |
| 51 |
< |
including bare metal nanoparticles, nanoparticles in an explicit |
| 52 |
< |
solvent, as well as clusters of liquid water. The predicted |
| 53 |
< |
mechanical properties of these systems are in good agreement with |
| 54 |
< |
experimental data and previous simulation work. |
| 45 |
> |
materials with different compressibilities. These systems are |
| 46 |
> |
problematic for traditional affine transform methods. The Langevin |
| 47 |
> |
Hull does not suffer from the edge effects of boundary potential |
| 48 |
> |
methods, and allows realistic treatment of both external pressure |
| 49 |
> |
and thermal conductivity due to the presence of an implicit solvent. |
| 50 |
> |
We apply this method to several different systems including bare |
| 51 |
> |
metal nanoparticles, nanoparticles in an explicit solvent, as well |
| 52 |
> |
as clusters of liquid water. The predicted mechanical properties of |
| 53 |
> |
these systems are in good agreement with experimental data and |
| 54 |
> |
previous simulation work. |
| 55 |
|
\end{abstract} |
| 56 |
|
|
| 57 |
|
\newpage |
| 126 |
|
have concentrations orders of magnitude lower than this in the |
| 127 |
|
cellular environment. The effective concentrations of single proteins |
| 128 |
|
in simulations may have significant effects on the structure and |
| 129 |
< |
dynamics of simulated structures. |
| 129 |
> |
dynamics of simulated systems. |
| 130 |
|
|
| 131 |
|
\subsection*{Boundary Methods} |
| 132 |
|
There have been a number of approaches to handle simulations of |
| 481 |
|
interactions between the pseudoatom cores. The $\sqrt{\rho_i}$ term in |
| 482 |
|
Eq. (\ref{eq:SCP1}) is an attractive many-body potential that models |
| 483 |
|
the interactions between the valence electrons and the cores of the |
| 484 |
< |
pseudo-atoms. $D_{ij}$, $D_{ii}$ set the appropriate overall energy |
| 484 |
> |
pseudo-atoms. $D_{ij}$ and $D_{ii}$ set the appropriate overall energy |
| 485 |
|
scale, $c_i$ scales the attractive portion of the potential relative |
| 486 |
|
to the repulsive interaction and $\alpha_{ij}$ is a length parameter |
| 487 |
|
that assures a dimensionless form for $\rho$. These parameters are |
| 494 |
|
In bulk gold, the experimentally-measured value for the bulk modulus |
| 495 |
|
is 180.32 GPa, while previous calculations on the QSC potential in |
| 496 |
|
periodic-boundary simulations of the bulk crystal have yielded values |
| 497 |
< |
of 175.53 GPa.\cite{QSC2} Using the same force field, we have performed |
| 498 |
< |
a series of 1 ns simulations on gold nanoparticles of three different radii under the Langevin Hull at a variety of applied pressures ranging from 0 -- 10 GPa. For the 40 \AA~ radius nanoparticle we obtain a value of 177.55 GPa for the bulk modulus of gold, in close agreement with both previous simulations and the experimental bulk modulus reported for gold single crystals.\cite{Collard1991} Polycrystalline gold has a reported bulk modulus of 220 GPa. The smaller gold nanoparticles (30 and 20 \AA~ radii) have calculated bulk moduli of 215.58 and 208.86 GPa, respectively, indicating that smaller nanoparticles approach the polycrystalline bulk modulus value while larger nanoparticles approach the single crystal value. As nanoparticle size decreases, the bulk modulus becomes larger and the nanoparticle is less compressible. This stiffening of the small nanoparticles may be related to their high degree of surface curvature, resulting in a lower coordination number of surface atoms relative to the the surface atoms in the 40 \AA~ radius particle. |
| 497 |
> |
of 175.53 GPa.\cite{QSC2} Using the same force field, we have |
| 498 |
> |
performed a series of 1 ns simulations on gold nanoparticles of three |
| 499 |
> |
different radii under the Langevin Hull at a variety of applied |
| 500 |
> |
pressures ranging from 0 -- 10 GPa. For the 40 \AA~ radius |
| 501 |
> |
nanoparticle we obtain a value of 177.55 GPa for the bulk modulus of |
| 502 |
> |
gold, in close agreement with both previous simulations and the |
| 503 |
> |
experimental bulk modulus reported for gold single |
| 504 |
> |
crystals.\cite{Collard1991} The smaller gold nanoparticles (30 and 20 |
| 505 |
> |
\AA~ radii) have calculated bulk moduli of 215.58 and 208.86 GPa, |
| 506 |
> |
respectively, indicating that smaller nanoparticles are somewhat |
| 507 |
> |
stiffer (less compressible) than the larger nanoparticles. This |
| 508 |
> |
stiffening of the small nanoparticles may be related to their high |
| 509 |
> |
degree of surface curvature, resulting in a lower coordination number |
| 510 |
> |
of surface atoms relative to the the surface atoms in the 40 \AA~ |
| 511 |
> |
radius particle. |
| 512 |
|
|
| 513 |
< |
We measure a gold lattice constant of 4.051 \AA~ using the Langevin Hull at 1 atm, close to the experimentally-determined value for bulk gold and the value for gold simulated using the QSC potential and periodic boundary conditions (4.079 \AA~ and 4.088\AA~, respectively).\cite{QSC2} The slightly smaller calculated lattice constant is most likely due to the presence of surface tension in the non-periodic Langevin Hull cluster, an effect absent from a bulk simulation. The specific heat of a 40 \AA~ gold nanoparticle under the Langevin Hull at 1 atm is 24.914 $\mathrm {\frac{J}{mol \, K}}$, which compares very well with the experimental value of 25.42 $\mathrm {\frac{J}{mol \, K}}$. |
| 513 |
> |
We obtain a gold lattice constant of 4.051 \AA~ using the Langevin |
| 514 |
> |
Hull at 1 atm, close to the experimentally-determined value for bulk |
| 515 |
> |
gold and the value for gold simulated using the QSC potential and |
| 516 |
> |
periodic boundary conditions (4.079 \AA~ and 4.088\AA~, |
| 517 |
> |
respectively).\cite{QSC2} The slightly smaller calculated lattice |
| 518 |
> |
constant is most likely due to the presence of surface tension in the |
| 519 |
> |
non-periodic Langevin Hull cluster, an effect absent from a bulk |
| 520 |
> |
simulation. The specific heat of a 40 \AA~ gold nanoparticle under the |
| 521 |
> |
Langevin Hull at 1 atm is 24.914 $\mathrm {\frac{J}{mol \, K}}$, which |
| 522 |
> |
compares very well with the experimental value of 25.42 $\mathrm |
| 523 |
> |
{\frac{J}{mol \, K}}$. |
| 524 |
|
|
| 525 |
|
\begin{figure} |
| 526 |
|
\includegraphics[width=\linewidth]{stacked} |
| 614 |
|
effects of the empty space due to the vapor phase; for this reason, we |
| 615 |
|
recommend using the number density (Eq. \ref{eq:BMN}) or number |
| 616 |
|
density fluctuations (Eq. \ref{eq:BMNfluct}) for computing |
| 617 |
< |
compressibilities. We achieved the best results using a sampling radius approximately 80\% of the cluster radius. This ratio of sampling radius to cluster radius excludes the problematic vapor phase on the outside of the cluster while including enough of the liquid phase to avoid poor statistics due to fluctuating local densities. |
| 617 |
> |
compressibilities. We achieved the best results using a sampling |
| 618 |
> |
radius approximately 80\% of the cluster radius. This ratio of |
| 619 |
> |
sampling radius to cluster radius excludes the problematic vapor phase |
| 620 |
> |
on the outside of the cluster while including enough of the liquid |
| 621 |
> |
phase to avoid poor statistics due to fluctuating local densities. |
| 622 |
|
|
| 623 |
< |
A comparison of the oxygen-oxygen radial distribution functions for SPC/E water simulated using the Langevin Hull and bulk SPC/E using periodic boundary conditions -- both at 1 atm and 300K -- reveals a slight understructuring of water in the Langevin Hull that manifests as a minor broadening of the solvation shells. This effect may be related to the introduction of surface tension around the entire cluster, an effect absent in bulk systems. As a result, molecules on the hull may experience an increased inward force, slightly compressing the solvation shell structure. |
| 623 |
> |
A comparison of the oxygen-oxygen radial distribution functions for |
| 624 |
> |
SPC/E water simulated using the Langevin Hull and bulk SPC/E using |
| 625 |
> |
periodic boundary conditions -- both at 1 atm and 300K -- reveals an |
| 626 |
> |
understructuring of water in the Langevin Hull that manifests as a |
| 627 |
> |
slight broadening of the solvation shells. This effect may be related |
| 628 |
> |
to the introduction of surface tension around the entire cluster, an |
| 629 |
> |
effect absent in bulk systems. As a result, molecules on the hull may |
| 630 |
> |
experience an increased inward force, slightly compressing the |
| 631 |
> |
solvation shell for these molecules. |
| 632 |
|
|
| 633 |
|
\subsection{Molecular orientation distribution at cluster boundary} |
| 634 |
|
|
| 736 |
|
\includegraphics[width=\linewidth]{RhoR} |
| 737 |
|
\caption{Density profiles of gold and water at the nanoparticle |
| 738 |
|
surface. Each curve has been normalized by the average density in |
| 739 |
< |
the bulk-like region available to the corresponding material. Higher applied pressures |
| 740 |
< |
de-structure both the gold nanoparticle surface and water at the |
| 741 |
< |
metal/water interface.} |
| 739 |
> |
the bulk-like region available to the corresponding material. |
| 740 |
> |
Higher applied pressures de-structure both the gold nanoparticle |
| 741 |
> |
surface and water at the metal/water interface.} |
| 742 |
|
\label{fig:RhoR} |
| 743 |
|
\end{figure} |
| 744 |
|
|
| 751 |
|
to prevent water from diffusing into the center of the gold |
| 752 |
|
nanoparticles. This behavior is likely not a realistic description of |
| 753 |
|
the real physics of the situation. A better model of the gold-water |
| 754 |
< |
adsorption behavior appears to require harder repulsive walls to |
| 755 |
< |
prevent this behavior. |
| 754 |
> |
adsorption behavior would require harder repulsive walls to prevent |
| 755 |
> |
this behavior. |
| 756 |
|
|
| 757 |
|
\section{Discussion} |
| 758 |
|
\label{sec:discussion} |
