| 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 |