Overview
This section covers the TIP5P force field for water as it is implemented into the towhee_ff_TIP5P file in the ForceFields directory.
All of the Towhee atom types for this force field are listed, along with a short description of their meanings.
TIP5P is a rigid water model that has a Lennard-Jones site on the Oxygen and bare charge sites on the Hydrogens and on the lone electron pairs.
TIP5P is a branched, rigid molecule (when you consider the lone pairs) and therefore configurational-bias
is currently unable to generate structures so you will need to provide an initial structure and also use the rotational-bias
moves instead of the configurational-bias moves for molecule transfers between simulation boxes.
The Jorgensen group has also parameterized two other water models that utilize different numbers of interaction sites.
Please see the TIP3Pand TIP4P web pages for more information
about those water models. This water model is suitable for use with any Lennard-Jones forcefield, but is
especially well suited for the OPLS family of potentials. Please see the OPLS-ua,
OPLS-aa, and OPLS-2001 web pages for more information.
Any discrepencies (especially typos) from the published force field values are the sole responsibility of Marcus Martin,
and I welcome feedback on how this implementation compares with other programs.
References for TIP5P
TIP5P was originally published in the following paper
TIP5P in Towhee
The official force field name for TIP5P in Towhee is 'TIP5P'.
This list contains the atom names for use in the towhee_input file, along with a brief description taken from the literature.
Please note that the capitalization and spacing pattern is important and must be followed exactly as listed here.
- 'H' : hydrogen in water
- 'O' : oxygen in water
- 'L' : charge site on lone pair positions in an sp3 hybridization of oxygen
Coulombic interactions
TIP5P uses atom-centered point charges to represent the electrostatic interactions.
Here are the suggested charge assignments for water.
- 'H': 0.241
- 'O': 0.0
- 'L': -0.241
Improper torsions
There are no improper torsions for this force field.
Return to the Towhee Capabilities web page
|