Semiempirical Tight-Binding

xTB

多数の分子・配座を量子化学的な情報を含む低コストモデルで評価し、より高精度な計算へ渡す候補を整えます。

位置づけ: 半経験的tight-binding
主な入力: XYZ/SDF、電荷、不対電子数
主な出力: エネルギー、構造、電荷、Hessian

1. 概要

xTBはextended tight-bindingに基づく半経験的電子状態計算プログラムです。ab initio法のようにすべての積分を同じ形で評価せず、parameterizationと近似によって大幅に高速化します。

RDKit力場より電子構造を反映した評価を低コストで行えるため、配座の再最適化、大規模系の初期構造修正、DFT前のscreeningに適します。ただし小さな選択性差や反応障壁の最終値には、対象に応じたbenchmarkが必要です。

2. GFN family

GFNはgeometry、frequency、noncovalent interactionを主目的としたparameterizationです。

方法	特徴	使いどころ
GFN0-xTB	非自己無撞着で高速	難しい初期構造の修正、大規模screening
GFN1-xTB	速度と電子的記述の中間	構造・振動数・一般的物性
GFN2-xTB	多極子静電相互作用等を含むSCC法	通常の高品質な構造最適化候補
GFN-FF	一般化力場	非常に大きな系の構造前処理

方法間のenergyを直接混ぜず、同一方法・同一条件内で比較します。元素、電荷状態、結合様式がparameterization範囲に適しているかも確認します。

3. Single-point計算

XYZ形式は元素と座標しか保持しないため、全電荷と不対電子数を明示します。--uhfは多重度ではなく不対電子数です。

Shell

xtb acetophenone.xyz --gfn 2 --chrg 0 --uhf 0

Energy同一Hamiltonianと条件での相対比較に用います。
Gradient初期構造の歪みや最適化の必要性を判断します。
Chargespopulation model依存のため、絶対的な原子電荷とはみなしません。
Orbital information近似準位であることを明記し、高精度法の軌道エネルギーと混同しません。

4. 構造最適化

--optでエネルギー勾配に基づく最適化を行います。品質指定は収束閾値と計算時間の両方に影響します。

Shell

xtb acetophenone.xyz --gfn 2 --opt tight --chrg 0 --uhf 0

# 主な生成物
# xtbopt.xyz : optimized geometry
# xtbopt.log : optimization trajectory
# xtbrestart : restart information

終了メッセージ、gradient norm、最終構造、意図しない解離・結合形成を確認します。SCCが難しい場合にGFN0-xTBで構造を修正してからGFN2-xTBへ進む方法は、公式資料でも案内されています。

5. Hessian、振動数、熱化学量

最適化構造に対してHessianを評価し、停留点の性質を確認します。

Shell

xtb xtbopt.xyz --gfn 2 --hess --chrg 0 --uhf 0

極小構造では大きな虚振動がないことを確認します。非常に低い振動数は柔軟なtorsionや数値誤差に敏感です。調和振動子・rigid rotor近似に基づく熱化学補正を、溶液中の厳密な自由エネルギーとみなしません。

6. ALPB・GBSA溶媒モデル

implicit solventは連続誘電体として溶媒効果を近似します。ALPBを使う例は次の通りです。

Shell

xtb acetophenone.xyz --gfn 2 --opt tight --alpb chloroform

気相と溶媒条件を混ぜない

配座順位、イオン対、分子間距離は溶媒モデルで変化します。比較する全構造に同じmodel・溶媒・温度条件を適用し、明示的溶媒が必要な現象かを検討します。

7. 実行と検証

RDKitから候補を出力
配座ID、formal charge、spinをmetadataに保持します。
xTBで再最適化
各jobを独立directoryで実行し、出力の上書きを防ぎます。
重複排除
最適化後にRMSDとenergyで同一配座へ収束した候補を統合します。
高精度法へ接続
energy window内の代表構造をGaussian等で再最適化します。

XYZ例をダウンロード

8. 参考文献・公式資料

xTB project, User Guide to Semiempirical Tight Binding.
C. Bannwarth, S. Ehlert, S. Grimme, “GFN2-xTB,” J. Chem. Theory Comput. 2019. DOI
S. Grimme, C. Bannwarth, P. Shushkov, “GFN-xTB,” J. Chem. Theory Comput. 2017. DOI

次章: Gaussian