The general form of the force field implemented in the program is:
where
is the sum of Lennard-Jones interactions taken over non-bonded atom
pairs. A pair of atoms is considered as non-bonded if they are on
different molecules or if they are on the same molecule but separated by more
than two covalent bonds. A case when a pair of atoms is separated by exactly
3 bonds (the so-called 1-4 neighbors) is treated separately, see below.
is the distance between atoms
and
. Parameters
and
for atom pair (
) are determined
by default from the individual
and
parameters
for atoms using combination rules; they can be also specified separately.
is the sum of electrostatic interactions
The remaining terms in (1) are inramolecular interactions due to covalent bonds, covalent angles and torsion angles. They are described below.
For each molecule type used in the simulation a file describing molecular
structure and parameters of the force field is needed. The
file must have extension *.mmol . Some examples of .mmol files are given
in moldb
directory.
.mmol
files consist of several parts which are described below.
Lines, beginning with "#" are commentaries and they are ignored by the
program. Other lines contain parameters in free format. The delimiter
between keywords and parameters is one or more spaces. Do not use Tab as
a delimiter, it may result in unpredictable behaviour.