454 |
|
density corrected SSD models can be found in reference |
455 |
|
\ref{Gezelter04}. |
456 |
|
|
457 |
< |
!!!Place a {\sc BASS} scheme showing SSD parameters around here!!! |
457 |
> |
\begin{lstlisting}[caption={[A simulation of {\sc ssd} water]An example file showing a simulation including {\sc ssd} water.},label={sch:ssd}] |
458 |
> |
|
459 |
> |
#include "water.mdl" |
460 |
|
|
461 |
+ |
nComponents = 1; |
462 |
+ |
component{ |
463 |
+ |
type = "SSD_water"; |
464 |
+ |
nMol = 864; |
465 |
+ |
} |
466 |
+ |
|
467 |
+ |
initialConfig = "liquidWater.init"; |
468 |
+ |
|
469 |
+ |
forceField = "DUFF"; |
470 |
+ |
|
471 |
+ |
/* |
472 |
+ |
* The reactionField flag toggles reaction |
473 |
+ |
* field corrections. |
474 |
+ |
*/ |
475 |
+ |
|
476 |
+ |
reactionField = false; // defaults to false |
477 |
+ |
dielectric = 80.0; // dielectric for reaction field |
478 |
+ |
|
479 |
+ |
/* |
480 |
+ |
* The following two flags set the cutoff |
481 |
+ |
* radius for the electrostatic forces |
482 |
+ |
* as well as the skin thickness of the switching |
483 |
+ |
* function. |
484 |
+ |
*/ |
485 |
+ |
|
486 |
+ |
electrostaticCutoffRadius = 9.2; |
487 |
+ |
electrostaticSkinThickness = 1.38; |
488 |
+ |
|
489 |
+ |
\end{lstlisting} |
490 |
+ |
|
491 |
+ |
|
492 |
|
\subsection{\label{sec:eam}Embedded Atom Method} |
493 |
|
|
494 |
|
Several other molecular dynamics packages\cite{dynamo86} exist which have the |