--- trunk/nanoglass/experimental.tex	2007/09/06 20:44:02	3221
+++ trunk/nanoglass/experimental.tex	2007/09/21 21:31:25	3228
@@ -1,3 +1,5 @@
+%!TEX root = /Users/charles/Desktop/nanoglass/nanoglass.tex
+
 \section{Computational Methodology} 
 \label{sec:details} 
 
@@ -5,12 +7,12 @@ $\mathrm{Ag}_6\mathrm{Cu}_4$.  Three different sizes o
 
 Cu-core / Ag-shell and random alloy structures were constructed on an
 underlying FCC lattice (4.09 {\AA}) at the bulk eutectic composition
-$\mathrm{Ag}_6\mathrm{Cu}_4$.  Three different sizes of nanoparticles
+$\mathrm{Ag}_6\mathrm{Cu}_4$.  All three compositions were considered although experimental results suggest that the random structure is the most likely composition after synthesis.\cite{Jiang:2005lr,gonzalo:5163} Three different sizes of nanoparticles
 corresponding to a 20 \AA radius (1961 atoms), 30 {\AA} radius (6603
 atoms) and 40 {\AA} radius (15683 atoms) were constructed.  These
 initial structures were relaxed to their equilibrium structures at 20
 K for 20 ps and again at 300 K for 100 ps sampling from a
-Maxwell-Boltzmann distribution at each temperature.
+Maxwell-Boltzmann distribution at each temperature.  
 
 To mimic the effects of the heating due to laser irradiation, the
 particles were allowed to melt by sampling velocities from the Maxwell
@@ -48,7 +50,12 @@ relevant simulations.
 metallic atoms and the solvent, the effective viscosity is a
 essentially a free parameter that must be tuned to give experimentally
 relevant simulations.
-
+\begin{figure}[htbp]
+\centering
+\includegraphics[width=\linewidth]{images/stochbound.pdf}
+\caption{Methodology for nanoparticle cooling. Equations of motion for metal atoms contained in the outer 4 {\AA} were determined by Langevins' Equations of motion. Metal atoms outside this region were allowed to evolve under Newtonian dynamics.}
+\label{fig:langevinSketch}
+\end{figure}
 The viscosity ($\eta$) can be tuned by comparing the cooling rate that
 a set of nanoparticles experience with the known cooling rates for
 those particles obtained via the laser heating experiments.
@@ -98,7 +105,8 @@ $(\mathrm{Wm^{-2}K^{-1}})$.\cite{XXX}
 Values for the interfacial conductance have been determined by a
 number of groups for similar nanoparticles and range from a low
 $87.5\times 10^{6}$ $(\mathrm{Wm^{-2}K^{-1}})$ to $120\times 10^{6}$
-$(\mathrm{Wm^{-2}K^{-1}})$.\cite{XXX}
+$(\mathrm{Wm^{-2}K^{-1}})$.\cite{hartlandPrv2007} Plech {\it et al.} reported a value for the interfacial conductance of $G=105\pm 15$ $(\mathrm{Wm^{-2}K^{-1}})$ and
+$G=130\pm 15$ $(\mathrm{Wm^{-2}K^{-1}})$ for Pt nanoparticles.\cite{plech:195423,PhysRevB.66.224301}
 
 We conducted our simulations at both ends of the range of
 experimentally-determined values for the interfacial conductance.
@@ -110,9 +118,9 @@ closer to the faster regime: $117\times 10^{6}$
 $(\mathrm{Wm^{-2}K^{-1}})$ was used.  Based on calculations we have
 done using raw data from the Hartland group's thermal half-time
 experiments on Au nanospheres, we believe that the true G values are
-closer to the faster regime: $117\times 10^{6}$
-$(\mathrm{Wm^{-2}K^{-1}})$.
+closer to the faster regime: $117\times 10^{6}$ $(\mathrm{Wm^{-2}K^{-1}})$.
 
+
 The rate of cooling for the nanoparticles in a molecular dynamics
 simulation can then be tuned by changing the effective solvent
 viscosity ($\eta$) until the nanoparticle cooling rate matches the
@@ -127,7 +135,7 @@ profiles from Eq. \ref{eq:laplacetransform} exactly.
 Fig. \ref{fig:images_cooling_plot}. It should be noted that the
 Langevin thermostat produces cooling curves that are consistent with
 Newtonian (single-exponential) cooling, which cannot match the cooling
-profiles from Eq. \ref{eq:laplacetransform} exactly.
+profiles from Eq. \ref{eq:laplacetransform} exactly. Fitting the Langevin cooling profiles to a single-exponential produces $\tau=25.576$ ps, $\tau=43.786$ ps, and $\tau=56.621$ ps for the 20, 30 and 40 {\AA} nanoparticles and a G of $87.5\times 10^{6}$ $(\mathrm{Wm^{-2}K^{-1}})$. The faster cooling G of $117\times 10^{6}$ $(\mathrm{Wm^{-2}K^{-1}})$ produced a $\tau=13.391$ ps, $\tau=30.426$ ps, $\tau=43.857$ ps for the 20, 30 and 40 {\AA} nanoparticles.
 
 \begin{figure}[htbp]
 \centering
@@ -191,7 +199,17 @@ metals.
 $2\alpha_{ij}$ and include up to the sixth coordination shell in FCC
 metals.
 
-\subsection{Sampling single-temperature configurations from a cooling
-trajectory} 
+%\subsection{Sampling single-temperature configurations from a cooling
+%trajectory} 
 
-ffdsafjdksalfdsa
+To better understand the structural changes occurring in the nanoparticles throughout the cooling trajectory, configurations were sampled at temperatures throughout the cooling trajectory. These configurations were then allowed to evolve under NVE dynamics to sample from the proper distribution in phase space. Figure \ref{fig:images_cooling_time_traces} illustrates this sampling.
+
+
+\begin{figure}[htbp]
+	\centering
+		\includegraphics[height=3in]{images/cooling_time_traces.pdf}
+	\caption{Illustrative cooling profile for the 40 {\AA} nanoparticle evolving under stochastic boundary conditions corresponding to $G=$$87.5\times 10^{6}$ $(\mathrm{Wm^{-2}K^{-1}})$. At temperatures along the cooling trajectory, configurations were sampled and allowed to evolve in the NVE ensemble. These subsequent trajectories were analyzed for structural features associated with bulk glass formation.}
+	\label{fig:images_cooling_time_traces}
+\end{figure}
+
+