-
-
Notifications
You must be signed in to change notification settings - Fork 1
Tips for force field selection in 3D optimization
When you run Optimize 3D, you choose a force field to calculate the molecule's energy. This page compares MMFF and UFF so you can pick the right one for your molecule.
| Feature | MMFF (Merck Molecular Force Field) | UFF (Universal Force Field) |
|---|---|---|
| Primary Goal | High-accuracy structure/energy for organic & medicinal molecules. | General-purpose estimation for any molecule. |
| Supported Elements | Common organic elements (C, H, N, O, F, Si, P, S, Cl, Br, I) and some ions (Na+, K+, Zn2+, etc.). Most metals are not supported. | Nearly the entire periodic table (H to Lr, elements 1-103). |
| Accuracy | Very high accuracy for molecules within its parameter set (i.e., organic molecules). | General-purpose. Can handle organometallics but is less accurate than MMFF for organics. |
-
Description:
- Developed by the Merck pharmaceutical company specifically for drug design.
- Parameters are rigorously tuned using a vast set of organic molecules to accurately reproduce high-level quantum chemistry calculations and experimental data.
-
Pros:
- Provides extremely high-accuracy structures and relative energies for organic and medicinal molecules composed of its supported elements (esp. C, H, N, O, P, S, and halogens).
-
Cons:
- Supported elements are limited. It cannot be used for molecules containing most metals or other unsupported elements.
MMFF94 and MMFF94s use the same core MMFF parameter set but have a subtle difference:
-
MMFF94:
- The full, validated version. It is designed to accurately evaluate molecular energies and conformations.
- In some cases (e.g., the -NH2 group in aniline), it may correctly calculate a slightly more stable "pyramidal" (non-planar) structure.
-
MMFF94s:
- The "s" often stands for "static."
- This is a modified version designed to be more robust (i.e., less likely to fail) during minimizations, especially for tasks like docking.
- It simplifies or constrains certain structures. For example, it forces nitrogen atoms in amide bonds or attached to aromatic rings to be treated as "planar."
- For general-purpose structure optimization or cleanup,
MMFF94sis often preferred as it tends to produce chemically intuitive structures (like planar amides).
-
Description:
- As its "Universal" name implies, its greatest strength is its coverage of nearly the entire periodic table.
- Parameters are generated automatically based on general properties for each element (atomic radius, ionization potential, hybridization, etc.).
-
Pros:
- Often the only choice for optimizing structures that MMFF does not support, such as metal complexes, inorganic compounds, and organometallics.
-
Cons:
- Because the parameters are general, its accuracy for common organic molecules (especially conformational energy differences) is generally lower than that of the highly-tuned MMFF.
The force field/backend selectors live in Settings > 3D Optimization Settings (sets the default for the Optimize 3D button and Ctrl+L), and in 3D Edit > Constrained Optimization... (per-run "Force Field" dropdown, defaulting to the Settings value). Right-click Optimize 3D for a one-off run with a different method without changing your default. See 3D structure optimization in the User Manual.
| Backend | Available Force Fields |
|---|---|
| RDKit (always available) | MMFF94s (default), MMFF94, UFF |
| Open Babel (optional dependency) | MMFF94s, MMFF94, UFF, GAFF (General Amber Force Field), Ghemical |
| Plugins | Some installed plugins register their own optimizers, which then appear alongside the built-in force fields in both menus above. |
Linux note: the
moleditpy-linuxpackage ships with Open Babel disabled (see Installation for Linux), so on Linux only the RDKit force fields (MMFF94s/MMFF94/UFF) — and any plugin-provided optimizers — are available; GAFF and Ghemical require Open Babel and won't appear.
If MMFF94/MMFF94s setup fails because the molecule contains an element the parameterization doesn't cover (most transition metals), MoleditPy automatically falls back to UFF rather than crashing — you'll typically want to just switch to UFF yourself for such molecules going forward.
- User Manual — 3D structure optimization in context
- Tips for 3D editing
- Tips for 2D to 3D conversion
- Installation for Linux — why GAFF/Ghemical aren't available there