background image
c CCLRC
Section 4.3
4.3.3
Adding Solvent to a Structure
The utility wateradd adds water from an equilibrated configuration of 256 SPC water molecules
at 300 K to fill out the MD cell. The utility solvadd fills out the MD box with single-site solvent
molecules from a fcc lattice. The FIELD files will then need to be edited to account for the solvent
molecules added to the file.
Hint: to save yourself some work in entering the non-bonded interactions variables involving solvent
sites to the FIELD file put two bogus atoms of each solvent type at the end of the CONNECT DAT
file (for AMBER force-fields) the utility ambforce will then evaluate all the non-bonded variables
required by DL POLY 3 . Remember to delete the bogus entries from the CONFIG file before
running DL POLY 3 .
4.3.4
Analysing Results
DL POLY 3 is not designed to calculate every conceivable property you might wish from a sim-
ulation. Apart from some obvious thermodynamic quantities and radial distribution functions,
it does not calculate anything beyond the atomic trajectories. You must therefore be prepared
to post-process the HISTORY file if you want other information. There are some utilities in the
DL POLY 3 package to help with this, but the list is far from exhaustive. In time, we hope to have
many more. Our users are invited to submit code to the DL POLY 3 public library to help with
this.
The utilities available are described in the DL POLY 2 Reference Manual. Users should also be
aware that many of these utilities are incorporated into the DL POLY Graphical User Interface
[
15
].
4.3.5
Choosing Ewald Sum Variables
4.3.5.1
Ewald sum and SPME
This section outlines how to optimise the accuracy of the Smoothed Particle Mesh Ewald sum
parameters for a given simulation..
As a guide to beginners DL POLY 3 will calculate reasonable parameters if the ewald precision
directive is used in the CONTROL file (see Section
5.1.1
). A relative error (see below) of 10
-6
is
normally sufficient so the directive
ewald precision 1d-6
will make DL POLY 3 evaluate its best guess at the Ewald parameters , kmax1, kmax2 and kmax3.
(The user should note that this represents an estimate, and there are sometimes circumstances
where the estimate can be improved upon. This is especially the case when the system contains a
strong directional anisotropy, such as a surface.) These four parameters may also be set explicitly
by the ewald sum directive in the CONTROL file. For example the directive
ewald sum 0.35 6 6 8
would set = 0.35 °
A
-1
, kmax1 = 6, kmax2 = 6 and kmax3 = 8. The quickest check on the accuracy
of the Ewald sum is to compare the coulombic energy (U ) and virial (W) in a short simulation.
70