Cerius²·Forcefield Engines



C       File Formats

This section includes information on:

Files documented at our website

Forcefield parameter files

Pf forcefield parameter file converter

Trajectory files


Files documented at our website

The types of trajectory files that originated with Insight® II (.his and .arc), as well as the forcefield files for forcefields that were originally supplied with Discover®, are documented at our website:

http://www.msi.com/doc/


Forcefield parameter files

Important

Only the forcefield files for forcefields that were originally supplied with Cerius2 are documented in this section.  

This section shows the format of the forcefield parameter files used by Cerius2. They are ASCII files divided into 16 sections, each starting with a left-justified keyword and ending with the left-justified keyword END. This makes the parameter files easy to read; new sections can also be added easily.

The VERSION section must come first and the ATOMTYPES section must come before any atom type label is used. Comments begin with a hash mark (#) at the beginning of the line and can be entered within any of the sections after the file header.

Not all parameters listed below are used by each forcefield. For information about a particular forcefield, see its parameter file. The parameter files are located in the Cerius2-Resources/FORCE_FIELD directory and its subdirectories.

Version number and product name

Keyword


VERSION

The product name and the version number appear here. For example: CERIUS2 2

Header information

Keyword


HEADER

Information for the user should be included here (such as forcefield and references). This is sent to standard output unchanged.

Preferences

Keyword


PREFERENCES

Variable Format
Terms

BONDS   (1X,1A,T40,L1)  
ANGLES   (1X,1A,T40,L1)  
COULOMB   (1X,1A,T40,L1)  
INVERSIONS   (1X,1A,T40,L1)  
TORSIONS   (1X,1A,T40,L1)  
UREY_BRADLEY   (1X,1A,T40,L1)  
HYDROGEN_BONDS   (1X,1A,T40,L1)  
HARMONIC_CONSTRAINTS   (1X,1A,T40,L1)  
DIAGONAL_VAN_DER_WAALS   (1X,1A,T40,L1)  
OFF_DIAGONAL_VAN_DER_WAALS   (1X,1A,T40,L1)  
GENERATE_UNDEFINED_TERMS   (1X,1A,T40,L1)  
Bonds

SHRINK_CH_BONDS   (1X,1A,T40,L1)  
SHRINK_CH_H_ATOM   (1X,1A,T40,1A)  
SHRINK_CH_FACTOR   (1X,1A,T40,F10.5)  
Torsions

SINGLE_TORSION   (1X,1A,T40,L1)  
SCALE_TORSIONS_ABOUT_COMMON_BOND   (1X,1A,T40,L1)  
EXOCYCLIC_TORSIONS_SCALE_FACTOR   (1X,1A,T40,F10.5)  
Inversions

SINGLE_INVERSION   (1X,1A,T40,L1)  
Hydrogen bonds

H-BOND-METHOD   (1X,1A,T40,1A)  
H-BOND-LIST   (1X,1A,T40,L1)  
H-BOND_DIRECT_RCUT   (1X,1A,T40,F10.5)  
H-BOND_SPLINE_DISTANCE_ON   (1X,1A,T40,F10.5)  
H-BOND_SPLINE_DISTANCE_OFF   (1X,1A,T40,F10.5)  
H-BOND_SPLINE_ANGLE_ON   (1X,1A,T40,F10.5)  
H-BOND_SPLINE_ANGLE_OFF   (1X,1A,T40,F10.5)  
Coulomb

DISTANCE_DEPENDENT_DIELECTRIC_ CONSTANT   (1X,1A,T40,L1)  
COU_DIELECTRIC_CONSTANT   (1X,1A,T40,F10.5)  
COU_INTER_CUT_OFF   (1X,1A,T40,F10.5)  
COU_SPLINE_OFF   (1X,1A,T40,F10.5)  
COU_SPLINE_ON   (1X,1A,T40,F10.5)  
EWALD_SUM_COU_ACCURACY   (1X,1A,T40,F10.5)  
EWALD_SUM_COU_ETA   (1X,1A,T40,F10.5)  
EWALD_SUM_COU_KCUT   (1X,1A,T40,F10.5)  
EWALD_SUM_COU_RCUT   (1X,1A,T40,F10.5)  
EWALD_SUM_COU_AUTO_OPT   (1X,1A,T40,L1)  
COU_EXCLUDE_1-2   (1X,1A,T40,L1)  
COU_EXCLUDE_1-3   (1X,1A,T40,L1)  
COU_EXCLUDE_1-4   (1X,1A,T40,L1)  
COU_1-4_SCALE_FACTOR   (1X,1A,T40,F10.5)  
COU_METHOD   (1X,1A,T40,1A)  
COU_DIRECT_CUT-OFF   (1X,1A,T40,F10.5)  
van der Waals

VDW_RADII_COMBINATION_RULE   (1X,1A,T40,1A)  
VDW_INTER_CUT_OFF   (1X,1A,T40,F10.5)  
VDW_EXCLUDE_1-2   (1X,1A,T40,L1)  
VDW_EXCLUDE_1-3   (1X,1A,T40,L1)  
VDW_EXCLUDE_1-4   (1X,1A,T40,L1)  
VDW_1-4_SCALE_FACTOR   (1X,1A,T40,F10.5)  
VDW_METHOD   (1X,1A,T40,1A)  
VDW_SPLINE_ON   (1X,1A,T40,F10.5)  
VDW_SPLINE_OFF   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_AUTO_OPT   (1X,1A,T40,L1)  
EWALD_SUM_VDW_ACCURACY   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_ETA   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_KCUT   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_RCUT   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_REP_CUT   (1X,1A,T40,F10.5)  
EWALD_SUM_VDW_RATIO   (1X,1A,T40,F10.5)  
General nonbond

MINIMUM_IMAGE   (1X,1A,T40,L1)  
NON_BOND_LIST   (1X,1A,T40,L1)  
NON_BOND_LIST_CUT_OFF   (1X,1A,T40,F10.5)  
Automatic setup

ASSIGN_MASS   (1X,1A,T40,L1)  
ASSIGN_CHARGE   (1X,1A,T40,L1)  
ASSIGN_HYBRIDIZATION   (1X,1A,T40,L1)  
ATOM_TYPE   (1X,1A,T40,L1)  
CALCULATE_BOND_ORDER   (1X,1A,T40,L1)  
IGNORE_UNDEFINED_TERMS   (1X,1A,T40,L1)  

Atom types

Keyword


ATOMTYPES

Variable Format
atomtype,element,mass,charge,hyb,implhyds,lonepairs   (1X,A8,3X,A2,3X,F10.5,F8.4,3I4)  

Diagonal van der Waals

Keyword


DIAGONAL_VDW 

Variable Format
ffi,'LJ_6_12 ',r0,d0   (1X,A5,3X,4X,A10,1X,1F10.4,E12.4)  
ffi,'EXPO_6 ',r0,d0,y   (1X,A5,3X,4X,A10,1X,1F10.4, E12.4,1F10.4)  
ffi,'VDW_MORSE ',r0,d0,y   (1X,A5,3X,4X,A10,1X,1F10.4, E12.4,1F10.4)  
ffi,'R6_ATTRACT',r0,d0   (1X,A5,3X,4X,A10,1X,1F10.4,E12.4)  
ffi,'REP_EXPO ',r0,d0,y   (1X,A5,3X,4X,A10,1X,1F10.4,E12.4,1F10.4)  
ffi,'IGNORE ',d0   (1X,A5,3X,4X,A10,1X,E12.4)  

Atom typing rules

Keyword


ATOM_TYPING_RULES

Variable Format
atomtype,element,hyb,ring,shell,factor   (1X,A5,3X,8X,A2,4I12)  

Off-diagonal van der Waals

Keyword


OFF_DIAGONAL_VDW

Variable Format
ffi,ffj,'LJ_6_12 ',r0,d0   (1X,A5,4X,A5,3X,4X,A10,1X,1F10.4,E12.4)  
ffi,ffj,'EXPO_6 ',r0,d0,y   (1X,A5,4X,A5,3X,4X,A10,1X,1F10.4,E12.4,1F10.4)  
ffi,ffj,'VDW_MORSE ',r0,d0,y   (1X,A5,4X,A5,3X,4X,A10,1X,1F10.4,E12.4,1F10.4)  
ffi,ffj,'R6_ATTRACT',r0,d0   (1X,A5,4X,A5,3X,4X,A10,1X,1F10.4,E12.4)  
ffi,ffj,'REP_EXPO ',r0,d0,y   (1X,A5,4X,A5,3X,4X,A10,1X,1F10.4,E12.4,1F10.4)  
ffi,ffj,'IGNORE '   (1X,A5,4X,A5,3X,4X,A10)  

Bond stretch

Keyword


BOND_STRETCH

Variable Format
ffi,ffj,'HARMONIC ',k0,r0   (1X,A5,4X,A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,'MORSE ',k0,r0,alpha,shift   (1X,A5,4X,A5,3X,4X,A10,1X,4F10.4)  
ffi,ffj,'CUBIC ',k0,r0,d   (1X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,'IGNORE '   (1X,A5,4X,A5,3X,4X,A10)  

Angle bend

Keyword


ANGLE_BEND

Variable Format
ffi,ffj,ffk,'COS_HARMON',k0,t0   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,ffk,'COS_PERIOD',k0,b,p   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,'FOURIER3 ',k0,t0   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,ffk,'SCHLEYER ',k0,t0,d   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,'QUARTIC ',k0,t0,d,e,f   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,5F10.4)  
ffi,ffj,ffk,'CFF ',k0,t0,d   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,'THETA_HARM',k0,t0   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,ffk,'MM2 ',k0,t0,k01h,t1h,t2h,d   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,6F10.4)  
ffi,ffj,ffk,'IGNORE '   (1X,A5,4X,A5,4X,A5,3X,4X,A10)  

Torsions

Keyword


TORSIONS

Variable Format
ffi,ffj,ffk,ffl,'DIHEDRAL ',kw,n,w0   (1X,3(A5,4X),A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,ffl,'IGNORE '   (1X,3(A5,4X),A5,3X,4X,A10)  

Inversions

Keyword


INVERSIONS

Variable Format
ffi,ffj,ffk,ffl,'UMBRELLA ',k0,t0   (1X,3(A5,4X),A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,ffk,ffl,'IT_JIKL ',k0,t0,n   (1X,3(A5,4X),A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,ffl,'IT_IJKL ',k0,t0   (1X,3(A5,4X),A5,3X,4X,A10,1X,2F10.4)  
ffi,ffj,ffk,ffl,'IGNORE '   (1X,3(A5,4X),A5,3X,4X,A10)  

Urey-Bradley

Keyword


UREY_BRADLEY

Variable Format
ffi,ffj,ffk,'HARMONIC ',kub,r0,d   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  
ffi,ffj,ffk,'IGNORE '   (1X,A5,4X,A5,4X,A5,3X,4X,A10)  

Hydrogen bonds

Keyword


HYDROGEN_BONDS

Variable Format
ffi,ffj,'MORSE ',de,re,y   (1X,A5,4X,A5,3X,4X,A10,1X,E12.4,2F10.4)  
ffi,ffj,'LJ_12_10 ',de,re   (1X,A5,4X,A5,3X,4X,A10,1X,E12.4,1F10.4)  
ffi,ffj,'IGNORE '   (1X,A5,4X,A5,3X,4X,A10)  

Universal forcefield parameter generator

Keyword


GENERATOR 

Variable Format
ffi, aradius, aangle, azstar, azeta, auenerg, auang, iaprd, iacis, atorbar, ElectNeg   (1X,A5,3X,2X,6F10.4,2I10,F10.4,F10.5)  

Stretch-bend-stretch cross terms

Keyword


STRETCH_BEND_STRETCH

Variable Format
ffi,ffj,ffk,'MM2 ',r01,r02,t0,kr01,kr02,k1h,c1h,c2h   (1X,A5,3X,1X,A5,3X,1X,A5,3X,4X,A10,1X,8F10.4)  
ffi,ffj,ffk,'COS_STRCH ',r01,r02,t0,kr01,kr02   (1X,A5,3X,2X,6F10.4,2I10,F10.4,F10.5)  
ffi,ffj,ffk,'ANG_STRCH ',r01,r02,t0,kr01,kr02   (1X,A5,3X,2X,6F10.4,2I10,F10.4,F10.5)  

Stretch-stretch cross terms

Keyword


STRETCH_STRETCH

Variable Format
ffi,ffj,ffk,'SIMPLE ',r01,r02,k12   (1X,A5,4X,A5,4X,A5,3X,4X,A10,1X,3F10.4)  

Default preferences for Dreiding and Universal forcefields

The default preferences (listed in the PREFERENCES section) for the Cerius2 general-purpose forcefields, Dreiding and Universal, are given in the table below.

Variable Dreiding Universal
BONDS   T   T  
ANGLES   T   T  
COULOMB   T   T  
INVERSIONS   T   T  
TORSIONS   T   T  
UREY_BRADLEY   F   F  
HYDROGEN_BONDS   T   F  
DIAGONAL_VAN_DER_WAALS   T   T  
OFF_DIAGONAL_VAN_DER_WAALS   T   T  
GENERATE_UNDEFINED_TERMS   F   T  
IGNORE_UNDEFINED_TERMS   T   F  
SHRINK_CH_BONDS   F   F  
SINGLE_TORSION   F   F  
SCALE_TORSIONS_ABOUT_COMMON_BOND   T   T  
EXOCYCLIC_TORSIONS_SCALE_FACTOR   0.40000   1.0000  
SINGLE_INVERSION   F   F  
H-BOND_METHOD   SPLINE   SPLINE  
H-BOND_LIST   T   T  
NON_BOND_LIST   T   T  
NON_BOND_LIST_CUT_OFF   ---   15.00000  
DISTANCE_DEPENDENT_DIELETRIC_CONSTANT   T   F  
COU_DIELETRIC_CONSTANT   1.00000   1.00000  
COU_INTER_CUT_OFF   ---   14.00000  
COU_SPLINE_OFF   ---   14.00000  
COU_SPLINE_ON   ---   11.00000  
COU_EXCLUDE_1-2   T   T  
COU_EXCLUDE_1-3   T   T  
COU_EXCLUDE_1-4   F   F  
COU_1-4_SCALE_FACTOR   1.00000   1.00000  
COU_METHOD   SPLINE   SPLINE  
COU_DIRECT_CUT_OFF   ---   11.00000  
VDW_RADII_COMBINATION_RULE   ARITHMETIC   GEOMETRIC  
VDW_INTER_CUT_OFF   ---   14.00000  
VDW_EXCLUDE_1-2   T   T  
VDW_EXCLUDE_1-3   T   T  
VDW_EXCLUDE_1-4   F   F  
VDW_1-4_SCALE_FACTOR   1.00000   1.00000  
VDW_METHOD   SPLINE   SPLINE  
VDW_SPLINE_ON   ---   11.00000  
VDW_SPLINE_OFF   ---   14.00000  
MINIMUM_IMAGE   F   F  
ASSIGN_MASS   T   T  
ASSIGN_CHARGE   F   F  
ASSIGN_HYBRIDIZATION   F   T  
ATOM_TYPE   T   T  
CALCULATE_BOND_ORDER   F   T  


Pf forcefield parameter file converter

Because forcefield parameter files produced by CERIUS or POLYGRAF cannot be read by Cerius2, all CERIUS forcefield parameter files and the Dreiding321 of POLYGRAF have been converted into the new Cerius2 format. These forcefields are available within the Cerius2-Resources/FORCE-FIELD directory and its subdirectories (see 1, "Using the Open Force Field").

If you have your own forcefield(s) created for use in CERIUS or POLYGRAF, you need to convert them to the new format before they can be loaded into Cerius2. To facilitate this a standalone file format converter program has been provided with this release of Cerius2.

The pf program allows you to convert previously created CERIUS or POLYGRAF forcefields into Cerius2 forcefields.

The program is not accessible through the GUI, but can be run from any directory by entering the following command at the system prompt:


>	pf [-CI or -G] IN_ff_parameter_file OUT_ff_parameter_file 

The optional arguments specify if the origin of the input file is CERIUS (-CI) or POLYGRAF (-G) and are a safeguard in case the converter cannot decipher the source of the parameter file.

For example, to convert a CERIUS forcefield file named myoldforcefield.par to a Cerius2 forcefield file named mynewforce-field1.01, type:


>	pf -CI myoldforcefield.par mynewforcefield1.01

If the energies obtained with the converted forcefield are different from those obtained with the original forcefield, then carefully examine the cross terms and options described below.

Options

Make sure that the following three options are chosen in both the old and the newly converted forcefield file.

POLYGRAF file conversion warnings

During conversion of a POLYGRAF forcefield parameter file you may notice a number of warning messages reported to the text window. This is because some options are not yet supported, or their underlying actions may be treated in a different way by Cerius2 than by POLYGRAF. When such an option is found by the converter, the parameter file section for that option is skipped.

Several such POLYGRAF options are listed in the following table:

Keyword Underlying action
DEL RE FF or DEL RE ATOMS   Del Re charges.  
GASTEIGER   Gasteiger charges.  
ANGX 2 K   Cosine angle-stretch cross term.  
ANGANGINV   Angle-angle inversion.  
BNDXANG   Bond cross angle terms.  
ANGXANG   Angle cross angle terms.  
ANGANGTOR   Coupling of 1-2-3 and 2-3-4 angles of torsions.  
BNDBNDTOR   Coupling of the 1-2 and 3-4 bonds of torsions.  
PITWIST   Cross term -- twist terms.  
CROSS   Cross bond-angle terms.  
USER   User-supplied parameters.  
POLYENE   Whether the bond and angle force constants depend on bond order.  
MASSZER   Zero-mass option.  
Q EQUIL   This action is dealt with by Cerius2 separately in the Charges module.  

Atom mass warning

During file conversion the following warning may appear in the text window:


WARNING: Some atom types masses read in differ from their expected mass. The expected mass is given by the element mass + Number of implicit hydrogens. Therefore some masses have been adjusted.

This happens when an atom type with implicit hydrogens is found during conversion. In this case, the converter first determines whether the masses of the implicit hydrogens have already been added to the atom's mass and, if not, it adds them.


Trajectory files

Documentation is provided at MSI's website.




Last updated July 09, 1998 at 07:49PM PDT.
Copyright © 1997, 1998 Molecular Simulations Inc. All rights reserved.