- 1.7
- 1.6
- 1.5
- Beta
- Examples (for 1.6a)
- Misc
Table of Contents
- Overview
- Supplied Forcefields
- Canongia-Lopes & Padua Ionic Liquids (cldp-il.ff)
- Youngs, Kohanoff, & Del Pópolo [dmim]Cl (dmimcl-fm.ff)
- Youngs & Hardacre [dmim]Cl (dmimcl-fm2.ff)
- Jorgensen at al. OPLS-AA (oplsaa.ff)
- Berensen et al.s Simple Point Charge Water (spc.ff)
- Berensen et al.s Extended Simple Point Charge Water (spce.ff)
- Rappe et al.s Universal Forcefield (testing/uff.ff)
- Mayo, Olafson & Goddard II's Generic Forcefield (testing/dreiding.ff)
- General Amber Forcefield (testing/gaff.ff)
- Liu, Wu & Wang's United-Atom Ionic Liquids Forcefield (lww-il.ff)
- Keyword Reference
- Rule-Based Forcefields
- Typing
- NETA Reference
- VDW Functional Forms
- Bond Functional Forms
- Angle Functional Forms
- Torsion Functional Forms
If you're doing anything interesting with a model or a molecule, a suitable forcefield description of the system is a must. A forcefield contains lists of parameters that describe the interactions between atoms, for example bonds, angles, and van der Waals interactions. More specifically, a forcefield contains parameters to describe many such interactions in many different types of molecule or chemical environment. An 'ex<x>pression', referred to throughout the manual, should be thought of as the subset of terms from a given forcefield necessary to describe all the interactions within a model.
Aten has its own free format for forcefield files, described in the following sections. Once loaded in, the energy and forces in Models can then be calculated, and allows energy minimisation etc. More so, once a set of forcefield parameters has been read in and used to describe a model, this expression can be written out using a custom format ready for input into something else.
The basic forcefield file format is designed to be as readable as possible by both machine and user. Lines are free format, meaning that any number of tabs, commas, and spaces may separate items. Text items should be enclosed with either double or single quotes if they themselves contain these delimiters (in particular, this applies to NETA descriptions). There is no terminating symbol at the end of a line, c.f. the command language where a ';' typically ends every command.
The majority of data is contained within blocks in the file. Blocks begin with a specific keyword (e.g. inter), contain one or more lines of definitions, and are terminated by an 'end' keyword. Forcefield files may contain many blocks of the same type, permitting terms of the same type but with differing functional forms to be defined easily.
More advanced forcefields may contain generator sections and functions that enable them to either generate all their parameters on the fly from a basic set of data, or fill in missing terms that are not defined in any blocks.
The format is keyword-based and as simple as possible. Most input is enclosed within blocks, beginning with a keyword and terminated after several lines of data with and 'end' statement. As an example of the overall structure of a forcefield, consider the simple point charge (SPC) forcefield for water as is provided with Aten:
name "SPC Water" units kj types 1 HW H "nbonds=1" 2 OW O "-H,-H" end inter lj 1 HW 0.41 0.0 0.0 2 OW -0.82 0.650 3.166 end bonds constraint HW OW 4184.0 1.000 end angles bondconstraint HW OW HW 4184.0 1.62398 end
After giving the forcefield a name and defining the energy units used for the parameters, the two types of atom described by the forcefield (OW and HW) are listed. A unique id, type name, base element, and type description are provided, providing Aten with all it needs to be able to recognise and type these atoms within a model. For each of these types the van der Waals data are provided in the Lennard-Jones style ('inter lj') - again, the type id and type name are specified, followed by the atomic charge and the epsilon and sigma values. Note here that there are default combination rules set for each functional form - see VDW functional forms for a list. Finally, the single bond and angle within the molecule are defined. The type names involved in the interactions are given, followed by the necessary parameters for the functional form specified ('constraint' for the bond, and 'bondconstraint' for the angle). And that's it. Detailed explanations of each section follow.

