Awesome
OpenFF Force Fields
This repository contains force fields released by the Open Force Field Initiative.
These force fields use the SMIRKS Native Open Force Field (SMIRNOFF) format.
By convention these files use the .offxml
file extension.
The SMIRNOFF format has a specification and is discussed in a JCTC publication and associated pre-print.
The OpenFF Toolkit provides a reference implementation of the SMIRNOFF format. In particular, the ForceField
class is used to load SMIRNOFF-format force fields and the create_openmm_system
method enables the parametrization of small molecules into OpenMM objects.
Detailed usage examples can be found in the OpenFF Toolkit repository.
Each mainline force field is currently available in two forms -- both with and without bond constraints to hydrogen. The default version of each force field (i.e. openff-2.0.0.offxml
) is suitable for typical molecular dynamics simulations with constrained bonds to hydrogen. The "unconstrained" version of each force field (i.e. openff_unconstrained-2.0.0.offxml
) should be used when single-point energies are a major concern (e.g. geometry optimizations) and when comparing the force field to QM data.
This repository may also contain other useful force fields, data, or utilities, separate from mainline OpenFF force fields. More information on those can be found in the docs/
directory in this repository.
How to cite
To cite the Parsley line of force fields (openff-1.Y.Z
) please use this citation.
To cite the Sage line of force fields (openff-2.Y.Z
) please use this citation.
Optionally, consider also citing the version-specific Zenodo DOI for the particular force field from the table below.
Details for each force field in this repository can be found in the following table:
Installation
conda install -c conda-forge openff-forcefields
Use
Installing this package exposes an entry point that makes the openforcefields/offxml/
directory easily accessible by other packages in the same Python installation. If the OpenFF Toolkit is installed, it will automatically detect and use this entry point when loading OFFXML files:
from openff.toolkit.typing.engines.smirnoff import ForceField
ff = ForceField('openff-2.0.0.offxml')
Otherwise, the entry point can be accessed by querying the openforcefield.smirnoff_forcefield_directory
entry point group.
from importlib.metadata import entry_points
for entry_point in iter_entry_points(group='openforcefield.smirnoff_forcefield_directory'):
print(entry_point.load()())
For convenience, the OpenFF Toolkit also has a standalone function get_available_force_fields
that lists the force fields discered by this entry point group:
>>> from openff.toolkit.typing.engines.smirnoff.forcefield import get_available_force_fields
>>> get_available_force_fields()
['smirnoff99Frosst-1.0.2.offxml', 'smirnoff99Frosst-1.0.0.offxml', 'smirnoff99Frosst-1.1.0.offxml', 'smirnoff99Frosst-1.0.4.offxml', 'smirnoff99Frosst-1.0.8.offxml', 'smirnoff99Frosst-1.0.6.offxml', 'smirnoff99Frosst-1.0.3.offxml', 'smirnoff99Frosst-1.0.1.offxml', 'smirnoff99Frosst-1.0.5.offxml', 'smirnoff99Frosst-1.0.9.offxml', 'smirnoff99Frosst-1.0.7.offxml', 'ff14sb_off_impropers_0.0.2.offxml', 'ff14sb_0.0.1.offxml', 'ff14sb_0.0.3.offxml', 'ff14sb_off_impropers_0.0.1.offxml', 'ff14sb_off_impropers_0.0.3.offxml', 'ff14sb_0.0.2.offxml', 'openff-1.0.1.offxml', 'openff-1.1.1.offxml', 'openff-1.0.0-RC1.offxml', 'openff-1.2.0.offxml', 'openff-1.3.0.offxml', 'openff_unconstrained-2.0.0-rc.1.offxml', 'openff_unconstrained-1.3.1.offxml', 'openff_unconstrained-1.2.1.offxml', 'openff-2.0.0-rc.2.offxml', 'openff_unconstrained-1.0.0-RC2.offxml', 'openff_unconstrained-1.1.0.offxml', 'openff_unconstrained-1.0.0.offxml', 'openff-2.0.0.offxml', 'openff_unconstrained-2.0.0.offxml', 'openff_unconstrained-2.0.0-rc.2.offxml', 'openff-1.1.0.offxml', 'openff-1.0.0.offxml', 'openff-1.0.0-RC2.offxml', 'openff-1.3.1.offxml', 'openff-1.2.1.offxml', 'openff-1.3.1-alpha.1.offxml', 'openff_unconstrained-1.0.0-RC1.offxml', 'openff_unconstrained-1.2.0.offxml', 'openff_unconstrained-1.3.0.offxml', 'openff-2.0.0-rc.1.offxml', 'openff_unconstrained-1.0.1.offxml', 'openff_unconstrained-1.1.1.offxml', 'openff_unconstrained-1.3.1-alpha.1.offxml']
History
Force fields in the Parsley and Sage lines are descended from the SMIRNOFF99Frosst line of force fields.
The first official release (openff-1.0.0.offxml
, code name "Parsley") was made in September 2019 as a result of the Open Force Field Initiative's refitting efforts.
General versioning guidelines
Applicable in general to SMIRNOFF-format FFs produced by the Open Force Field Consortium
Force fields moving forward will be called name-X.Y.Z
X
denotes some major change in functional form or fitting strategy.Y
is the parameterization epoch / generation, or a minor change that can affect energy.Z
is a bugfix version -- e.g. something we've caught and corrected.
Versions
-
v1.0.0 Parsley
: First major forcefield release. -
v1.0.1 Parsley
: This bugfix release contains following changes: (1) Addition of monatomic ionLibraryCharges
. -
v1.1.0 Parsley
: This minor release contains following changes: (1) Addition of new proper torsions and improper torsions for tetrazole; (2) Corrections to N-N bond rotation periodicity; (3) Removal of redundant periodicity component int19
; (4) Addition of three new bond and angle terms,a22a
,b14a
andb36a
. -
v1.1.1 Parsley
: This bugfix release contains following changes: (1) Addition of monatomic ionLibraryCharges
. -
v1.2.0 Parsley
: This minor release contains following changes: (1) New, carefully designed quantum chemical dataset was utilized in training valence parameters in the force field and (2) Removal of redundancy int108
SMIRKS pattern -
v1.2.1 Parsley
: This bugfix release manually changes two bond force constants to resolve an issue seen in propyne substituents when using hydrogen mass repartitioning with a 4fs timestep. Full details are available at https://github.com/openforcefield/openforcefields/issues/19. -
v1.3.0 Parsley
: This minor release contains a fix of amide-related issues; (1) a poor performance of v1.2 in reproducing amide torsional energy profiles and (2) absence of appropriate torsion parameters for dialkyl amides. -
v1.3.1-alpha.1 Parsley
: The new force field files in this release are adapted from the openff-1.3.0 release, except that two angle parameters have been reverted to their original values from smirnoff99Frosst 1.1.0, in an attempt to fix sulphonamide geometries. -
v1.3.1 Parsley
: This release is identical tov1.3.1-alpha.1 Parsley
. -
v2.0.0-rc.1 Sage
: This major release candidate contains both refit valence and vdW terms. Full details are available at https://github.com/openforcefield/openff-sage/releases/tag/2.0.0-rc.1 -
v2.0.0-rc.2 Sage
: This major release candidate is identical tov2.0.0-rc.1 Sage
except that theangle
value fora16
has been changed to180.0 * degree
, as the previous value of183... * degree
causes geometry optimizers to fail to converge. -
v2.0.0 Sage
: This major release contains the same physical parameters asv2.0.0-rc.2 Sage
, but has the parameter ids changed. For more information see the openff-sage repository. -
2023.02.1
: This release switches to a calendar-versioning system and addstip3p-1.0.0.offxml
andtip3p-fb-1.0.0.offxml
, as well astip3p.offxml
andtip3p-fb.offxml
, which will always point to the latest files in their respective lines. -
2023.04.1
: This release addsopenff-2.1.0-rc.1.offxml
Sage 2.1.0 Release Candidate 1, with several small parameter fixes and improvements. -
2023.05.0
: This release adds theopc-1.0.0.offxml
, our initial port of the OPC water model. -
2023.05.1
: This release addsopenff-2.1.0.offxml
, Sage 2.1.0, with identical physics parameters toopenff-2.1.0-rc.1.offxml
. -
2023.06.0
: This release addsopc-1.0.1.offxml
, standardizing the use of the string "kilocalorie_per_mole". It also addstip3p-fb-1.1.0.offxml
, which fixes a unit error in the geometric constraints and adds additional ion parameters. Finally, it addsopc3-1.0.0.offxml
andtip4p_fb-1.0.0.offxml
. -
2024.01.0
: This release addsopenff-2.1.1.offxml
andopenff_unconstrained-2.1.1.offxml
, Sage 2.1.1, which has identical parameters to Sage 2.1.0 (originally released in version2023.05.1
of theopenff-forcefields
package) but adds Xe van der Waals parameters from Tang, K.T., Toennies, J.P. New combining rules for well parameters and shapes of the van der Waals potential of mixed rare gas systems. Z Phys D - Atoms, Molecules and Clusters 1, 91–101 (1986). -
2024.03.0
: This release addsopenff-2.2.0-rc1.offxml
andopenff_unconstrained-2.2.0-rc1.offxml
, Sage 2.2.0 Release Candidate 1, which modifies some small ring internal angles and remedies issues with sulfamide geometries. -
2024.04.0
: This release addsopenff-2.2.0.offxml
andopenff_unconstrained-2.2.0.offxml
, Sage 2.2.0, with identical parameters toopenff-2.2.0-rc1.offxml
. -
2024.07.0
: This release addsopenff-2.2.1-rc1.offxml
andopenff_unconstrained-2.2.1-rc1.offxml
, which contains re-fit valence parameters where linear angles have been pinned to linear equilibrium values. -
2024.09.0
This release adds the Sage 2.2.1 force field withopenff-2.2.1.offxml
andopenff_unconstrained-2.2.1.offxml
, with identifical parameters to 2.2.1-rc1.
Acknowledgements
Project based on the Computational Molecular Science Python Cookiecutter version 1.0.