Methods#
The plug-in exposes four Hamiltonians / force fields from the
xTB family. The xTB method parameter is shared by all
sub-steps.
Choosing among the GFN family#
GFN2-xTB (default)#
The recommended general-purpose method. Self-consistent with multipole electrostatics (up to quadrupoles on each atom) and a density-dependent dispersion correction (D4 family). Parametrized for Z = 1-86.
GFN2-xTB is the right starting point for almost everything: organic, inorganic, organometallic, and main-group systems. It is the most accurate of the GFN methods for structures, conformer energies, and non-covalent interactions, and is well-tested across a broad range of benchmark sets.
Citation: Bannwarth, Ehlert, Grimme, J. Chem. Theory Comput. 2019, 15, 1652.
GFN1-xTB#
The earlier self-consistent xTB Hamiltonian. Uses only monopole-monopole electrostatics and a less elaborate dispersion correction. Parametrized for Z = 1-86.
GFN1-xTB is typically less accurate than GFN2 for thermochemistry but is sometimes more robust for difficult electronic-structure cases (heavily strained systems, unusual oxidation states). Useful when GFN2 fails to converge or gives implausible geometries.
Citation: Grimme, Bannwarth, Shushkov, J. Chem. Theory Comput. 2017, 13, 1989.
GFN0-xTB#
A non-self-consistent (single-shot) method. Faster than GFN1/GFN2, and very robust as a starting-point method or for very large systems. Parametrized for Z = 1-86.
Best uses: rough screening of large libraries; getting a reasonable starting geometry for a higher-level method; cases where SCC convergence is unreliable.
Note that GFN0-xTB is not parametrized for ALPB or GBSA implicit solvation. If you need implicit solvation with GFN0, use CPCM-X. The plug-in will pass the combination to xtb regardless; xtb itself will emit a warning or refuse to run.
Citation: Pracht, Caldeweyher, Ehlert, Grimme, ChemRxiv 2019, DOI 10.26434/chemrxiv.8326202.v1. (As of the plug-in v1 release this paper has only a preprint; a peer-reviewed version may appear later.)
GFN-FF#
A generic, non-quantum force field automatically parametrized by xTB from the structure. No SCC, no Hessian build-up of polarization response.
GFN-FF is by far the fastest of the four and is intended for very large systems (thousands of atoms) where even GFN0 is too slow, and for MD or conformational sampling where a force field is sufficient. Accuracy is markedly lower than the SCC GFN methods – do not use GFN-FF for energy differences where you care about precision under a kcal/mol.
Citation: Spicher, Grimme, Angew. Chem. Int. Ed. 2020, 59, 15665.
Quick comparison#
Method |
Type |
Speed |
Best for |
|---|---|---|---|
GFN-FF |
Force field |
Fastest |
Very large systems, MD |
GFN0-xTB |
Non-SCC |
Fast |
Screening, starting point |
GFN1-xTB |
SCC |
Moderate |
Backup when GFN2 struggles |
GFN2-xTB |
SCC + D4 |
Moderate |
General-purpose default |
Accuracy parameter#
The Accuracy field maps to xTB’s --acc flag. Lower values
tighten the integral cutoffs and SCC convergence thresholds; higher
values loosen them. The default of 1.0 is appropriate for almost
all calculations. A typical usable range is roughly 0.0001 (very
tight) to 1000 (very loose, only for crude screening).
You generally do not need to change this. The main case where
tightening it helps is when you are about to do a
Frequencies calculation and want a tightly
converged starting geometry; the tight/vtight/extreme
optimization levels already arrange for this.