88 |
|
is often unclear what shape of gradient should be imposed at the |
89 |
|
boundary between materials. |
90 |
|
|
91 |
< |
% \begin{figure} |
92 |
< |
% \includegraphics[width=\linewidth]{figures/VSS} |
93 |
< |
% \caption{Schematics of periodic (left) and non-periodic (right) |
94 |
< |
% Velocity Shearing and Scaling RNEMD. A kinetic energy or momentum |
95 |
< |
% flux is applied from region B to region A. Thermal gradients are |
96 |
< |
% depicted by a color gradient. Linear or angular velocity gradients |
97 |
< |
% are shown as arrows.} |
98 |
< |
% \label{fig:VSS} |
99 |
< |
% \end{figure} |
91 |
> |
\begin{figure} |
92 |
> |
\includegraphics[width=\linewidth]{figures/npVSS} |
93 |
> |
\caption{Schematics of periodic (left) and non-periodic (right) |
94 |
> |
Velocity Shearing and Scaling RNEMD. A kinetic energy or momentum |
95 |
> |
flux is applied from region B to region A. Thermal gradients are |
96 |
> |
depicted by a color gradient. Linear or angular velocity gradients |
97 |
> |
are shown as arrows.} |
98 |
> |
\label{fig:VSS} |
99 |
> |
\end{figure} |
100 |
|
|
101 |
|
Reverse Non-Equilibrium Molecular Dynamics (RNEMD) methods impose an |
102 |
|
unphysical {\it flux} between different regions or ``slabs'' of the |
132 |
|
We have extended the VSS method for use in {\it non-periodic} |
133 |
|
simulations, in which the ``slabs'' have been generalized to two |
134 |
|
separated regions of space. These regions could be defined as |
135 |
< |
concentric spheres (as in figure \ref{fig:VSS}), or one of the regions |
135 |
> |
concentric spheres (as in figure \ref{fig:npVSS}), or one of the regions |
136 |
|
can be defined in terms of a dynamically changing ``hull'' comprising |
137 |
|
the surface atoms of the cluster. This latter definition is identical |
138 |
|
to the hull used in the Langevin Hull algorithm. |
426 |
|
Table \ref{table:couple} shows the calculated rotational friction coefficients $\Xi^{rr}$ for spherical gold nanoparticles and a prolate ellipsoidal gold nanorod in TraPPE-UA hexane. An angular momentum flux was applied between the A and B regions defined as the gold structure and hexane molecules beyond a certain radius, respectively. The resulting angular velocity gradient causes the gold structure to rotate about the prescribed axis. |
427 |
|
|
428 |
|
\begin{longtable}{lcccc} |
429 |
< |
\caption{Comparison of rotational friction coefficients under ideal ``stick'' conditions ($\Xi^{rr}_{stick}$) calculated via Stokes' and Perrin's laws and effective rotational friction coefficients ($\Xi^{rr}_{\mathit{eff}}$) of gold nanostructures solvated in TraPPE-UA hexane at 230 K. The ellipsoid is oriented with the long axis along the $z$ direction.} |
429 |
> |
\caption{Comparison of rotational friction coefficients under ideal ``stick'' conditions ($\Xi^{rr}_{\mathit{stick}}$) calculated via Stokes' and Perrin's laws and effective rotational friction coefficients ($\Xi^{rr}_{\mathit{eff}}$) of gold nanostructures solvated in TraPPE-UA hexane at 230 K. The ellipsoid is oriented with the long axis along the $z$ direction.} |
430 |
|
\\ \hline \hline |
431 |
< |
{Structure} & {Axis of Rotation} & {$\Xi^{rr}_{stick}$} & {$\Xi^{rr}_{\mathit{eff}}$} & {$\Xi^{rr}_{\mathit{eff}}$ / $\Xi^{rr}_{stick}$}\\ |
431 |
> |
{Structure} & {Axis of Rotation} & {$\Xi^{rr}_{\mathit{stick}}$} & {$\Xi^{rr}_{\mathit{eff}}$} & {$\Xi^{rr}_{\mathit{eff}}$ / $\Xi^{rr}_{\mathit{stick}}$}\\ |
432 |
|
{} & {} & {\small(amu A$^2$ fs$^{-1}$)} & {\small(amu A$^2$ fs$^{-1}$)} & \\ \hline |
433 |
|
Sphere (r = 20 \AA) & {$x = y = z$} & {3314} & {2386} & {0.720}\\ |
434 |
|
Sphere (r = 30 \AA) & {$x = y = z$} & {11749} & {8415} & {0.716}\\ |
439 |
|
\label{table:couple} |
440 |
|
\end{longtable} |
441 |
|
|
442 |
+ |
|
443 |
+ |
|
444 |
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |
445 |
|
% **DISCUSSION** |
446 |
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |