--- trunk/langevinHull/langevinHull.tex	2011/01/17 21:50:02	3715
+++ trunk/langevinHull/langevinHull.tex	2011/01/27 00:00:20	3716
@@ -17,7 +17,7 @@
 \setlength{\abovecaptionskip}{20 pt}
 \setlength{\belowcaptionskip}{30 pt}
 
-\bibpunct{[}{]}{,}{s}{}{;}
+\bibpunct{}{}{,}{s}{}{;}
 \bibliographystyle{achemso}
 
 \begin{document}
@@ -378,7 +378,7 @@ integrator in our code, OpenMD.\cite{Meineke2005,openm
 configurations, so this appears to be a reasonably good approximation.
 
 We have implemented this method by extending the Langevin dynamics
-integrator in our code, OpenMD.\cite{Meineke2005,openmd}  At each
+integrator in our code, OpenMD.\cite{Meineke2005,open_md}  At each
 molecular dynamics time step, the following process is carried out:
 \begin{enumerate}
 \item The standard inter-atomic forces ($\nabla_iU$) are computed.
@@ -400,7 +400,7 @@ using calls to the qhull library.\cite{Qhull} There is
 \item Atomic positions and velocities are propagated.
 \end{enumerate}
 The Delaunay triangulation and computation of the convex hull are done
-using calls to the qhull library.\cite{Qhull} There is a minimal
+using calls to the qhull library.\cite{Q_hull} There is a minimal
 penalty for computing the convex hull and resistance tensors at each
 step in the molecular dynamics simulation (roughly 0.02 $\times$ cost
 of a single force evaluation), and the convex hull is remarkably easy
@@ -432,10 +432,10 @@ the middle of the cluster. $N$ is the average number o
 \label{eq:BMN}
 \end{equation}
 The region we used is a spherical volume of 20 \AA\ radius centered in
-the middle of the cluster. $N$ is the average number of molecules
+the middle of the cluster with a roughly 25 \AA\ radius. $N$ is the average number of molecules
 found within this region throughout a given simulation. The geometry
-and size of the region is arbitrary, and any bulk-like portion of the
-cluster can be used to compute the compressibility.
+of the region is arbitrary, and any bulk-like portion of the
+cluster can be used to compute the compressibility. 
 
 One might assume that the volume of the convex hull could simply be
 taken as the system volume $V$ in the compressibility expression
@@ -489,19 +489,16 @@ metals.\cite{PhysRevB.59.3527,QSC}
 energy, and elastic moduli for FCC transition metals. The quantum
 Sutton-Chen (QSC) formulation matches these properties while including
 zero-point quantum corrections for different transition
-metals.\cite{PhysRevB.59.3527,QSC}
+metals.\cite{PhysRevB.59.3527,QSC2}
 
 In bulk gold, the experimentally-measured value for the bulk modulus
 is 180.32 GPa, while previous calculations on the QSC potential in
 periodic-boundary simulations of the bulk crystal have yielded values
-of 175.53 GPa.\cite{QSC} Using the same force field, we have performed
-a series of 1 ns simulations on 40 \AA~ radius
-nanoparticles under the Langevin Hull at a variety of applied
-pressures ranging from 0 -- 10 GPa.  We obtain a value of 177.55 GPa
-for the bulk modulus of gold using this technique, in close agreement
-with both previous simulations and the experimental bulk modulus of
-gold.
+of 175.53 GPa.\cite{QSC2} Using the same force field, we have performed
+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.
 
+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}}$.
+
 \begin{figure}
 \includegraphics[width=\linewidth]{stacked}
 \caption{The response of the internal pressure and temperature of gold
@@ -530,7 +527,7 @@ Langevin Hull simulations for pressures between 1 and 
 ensembles) have yielded values for the isothermal compressibility that
 agree well with experiment.\cite{Fine1973} The results of two
 different approaches for computing the isothermal compressibility from
-Langevin Hull simulations for pressures between 1 and 6500 atm are
+Langevin Hull simulations for pressures between 1 and 3000 atm are
 shown in Fig. \ref{fig:compWater} along with compressibility values
 obtained from both other SPC/E simulations and experiment.
 
@@ -594,8 +591,10 @@ compressibilities.
 effects of the empty space due to the vapor phase; for this reason, we
 recommend using the number density (Eq. \ref{eq:BMN}) or number
 density fluctuations (Eq. \ref{eq:BMNfluct}) for computing
-compressibilities.
+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.
 
+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.
+
 \subsection{Molecular orientation distribution at cluster boundary}
 
 In order for a non-periodic boundary method to be widely applicable,
@@ -604,7 +603,7 @@ fixing and characterizing the effects of artifical bou
 methods (e.g. hydrophobic boundary potentials) induced spurious
 orientational correlations deep within the simulated
 system.\cite{Lee1984,Belch1985} This behavior spawned many methods for
-fixing and characterizing the effects of artifical boundaries
+fixing and characterizing the effects of artificial boundaries
 including methods which fix the orientations of a set of edge
 molecules.\cite{Warshel1978,King1989}
 
@@ -680,7 +679,7 @@ To further test the method, we simulated gold nanopart
 
 \subsection{Heterogeneous nanoparticle / water mixtures}
 
-To further test the method, we simulated gold nanopartices ($r = 18$
+To further test the method, we simulated gold nanoparticles ($r = 18$
 \AA) solvated by explicit SPC/E water clusters using a model for the
 gold / water interactions that has been used by Dou {\it et. al.} for
 investigating the separation of water films near hot metal