トポロジー解析の理論:QTAIM と ELF

トポロジー解析の理論:QTAIM と ELF#

電荷分割解析の章 は密度を原子間で分ける方法の族を導出しました — そのどれもが分割の選択でした。本章は代わりに電子密度(および密接に関連するスカラー場)を直接、 それ自身の臨界点・盆地・局在構造を持つ真の数学的対象として研究します — \(\rho(\mathbf r)\) 自体の トポロジーから化学的意味(結合・孤立電子対・環・ケージ)を抽出し、分割の決定はまったく行いません。 これが Bader の**分子内原子の量子論(QTAIM)と Becke-Edgecombe の電子局在化関数(ELF)**であり、 crates/qc-prop/src/{qtaim,basin,realspace}.rs に直接基づきます。

QTAIM:トポロジー的対象としての密度#

電子密度は、化学的に決定的な単一の構造的事実を持つ、なめらかで至るところ定義されたスカラー場です: それはほぼ厳密に各原子核で局所極大を持ち(核カスプ)、それらの極大の対を結ぶ鞍点が、構造化学者が 手で描くまさに化学結合パスと一致することが判明します。QTAIM はこれを、\(\nabla\rho\) の臨界点 — \(\nabla\rho=\mathbf 0\) となる点 — を研究し、そこでの密度 Hessian の固有値符号によって各点を 分類することで厳密にします。

曲率符号による臨界点の分類#

臨界点で、\(\rho\) の(実対称)Hessian は 3 つの実固有値を持ちます;化学的に意味のある臨界点はすべて 非退化です(3 つの固有値すべてが非ゼロ — 真の局所極値または鞍点であり、平坦/退化した方向では ない);その種類は、3 つの固有値のうちいくつが負であるかだけで決まります:

符号 \((3,\sigma)\)

負の固有値数

種類

化学的意味

\((3,-3)\)

3

原子核(NCP)

密度極大 — (ほぼ厳密に)原子核上にある

\((3,-1)\)

2

結合(BCP)

結合した 2 原子間の鞍点 — 結合に垂直な 2 方向で密度が減少し、結合に沿って増加する

\((3,+1)\)

1

環(RCP)

結合の環内での密度極小

\((3,+3)\)

0

ケージ(CCP)

環に囲まれた真の局所密度極小

\(\sigma\)(トレースの符号数え:正の固有値ごとに \(+1\)、負の固有値ごとに \(-1\))とランク(これらの非 退化型では常に 3)を合わせたものが伝統的な \((3,\sigma)\) 記法です;結合 CP の定義的性質 — 結合パスに 垂直な 2 つの負の曲率とそれに沿った 1 つの正の曲率 — はまさに「これは 2 つの峰の間の峠だ」という 数学的言明であり、この枠組みにおける真の化学結合のトポロジー的シグネチャそのものです、単なる示唆的 アナロジーではありません。

Poincaré-Hopf 不変量:無料の分子全体整合性検査#

注目すべきトポロジー的事実が、どの分子であるかによらず、完全な孤立分子系についてすべての臨界点数 を結びつけます:

\[ n_{\text{NCP}} - n_{\text{BCP}} + n_{\text{RCP}} - n_{\text{CCP}} = 1. \]

これは密度の勾配ベクトル場に適用された Poincaré-Hopf 指数定理であり、任意の正しく完全に特徴付け られた分子トポロジーで成り立ちます — 水の 3 原子核 + 2 結合 + 0 環 + 0 ケージの臨界点は \(3-2+0-0=1\);ベンゼンの 12 原子核 + 12 結合 + 1 環 + 0 ケージは \(12-12+1-0=1\);四面体 \(P_4\) ケージの 4 原子核 + 6 結合 + 4 環 + 1 ケージは \(4-6+4-1=1\) を与えます。他のどんな値も、臨界点が見逃された (検索が見つけられなかった期待される結合点や環点)か、偽のものが見つかったことを意味します — この 恒等式は、比較すべき参照計算をまったく必要とせず内部一貫性だけで済む、数値臨界点検索に対する真に 無料の正しさ検査です。

結合臨界点の energetics:局所ビリアル定理#

結合臨界点の位置と密度値だけで既に定性的な情報を運びますが、Bader の局所ビリアル定理はその 単一の点で真のエネルギー密度を抽出し、純粋にトポロジー的な分類をエネルギー的な分類に変えます。 この定理は密度の Laplacian を局所運動エネルギー密度とポテンシャルエネルギー密度に関係づけます:

\[ \tfrac14\nabla^2\rho = 2G+V, \qquad G(\mathbf r) = \tfrac12\sum_{\mu\nu}D_{\mu\nu}\nabla\varphi_\mu\cdot \nabla\varphi_\nu\ (\ge0), \]

\(G\) は(常に正の)Lagrangian 運動エネルギー密度です。この 1 つの関係は、臨界点で既に分かっている \(\rho\)、\(\nabla^2\rho\)、\(G\) と合わせて、ポテンシャルエネルギー密度と全局所エネルギー密度を直接 解くのに十分です:

\[ V = \tfrac14\nabla^2\rho - 2G\ (\text{結合領域では}\le0), \qquad H = G+V = \tfrac14\nabla^2\rho - G. \]

これらから構築される 2 つの無次元比が、Hessian 符号だけからではなく局所 energetics から直接結合の 性格を分類します:結合性格比 \(|V|/G\)(経験的に、\(>2\) は共有(電子共有殻)相互作用、\(<1\) は 閉殻/イオン性あるいは van der Waals 的なもの、その間は部分的共有性)と結合次数 \(H/\rho\)(共有 結合では負 — 全局所エネルギー密度が安定化的 — 閉殻相互作用では正)です。

検証済み例 — 水の完全な QTAIM トポロジー、自動的に計算・分類:

import qc
water = "O 0 0 0.117; H 0 0.757 -0.469; H 0 -0.757 -0.469"
m = qc.chk.new(atom=water, ao="def2-svp", unit="angstrom").scf(ref="r").run()

topo = m.prop.qtaim.topology()
topo["counts"]          # {'nuclear': 3, 'bond': 2, 'ring': 0, 'cage': 0}
topo["poincare_hopf"]    # {'sum': 1, 'holds': True}

bcp = topo["critical_points"][3]   # 1 つの O-H 結合臨界点
bcp["rho"], bcp["laplacian"]        # (0.3659441002679062, -2.5397077300947384)
bcp["v_over_g"], bcp["h_over_rho"]  # (11.152930010504011, -1.9245989766373959)

両方の O–H 結合 CP(分子の \(C_{2v}\) 対称性により等価)は \(|V|/G\approx11.2\gg2\) と \(H/\rho\approx-1.92<0\) を与えます — 曖昧さのない共有結合分類であり、O–H 単結合について化学的直感 が要求するまさにその通りですが、今度はどんな分割スキームからでもなく、単一の数学的に well-defined な点での厳密な局所エネルギー密度計算から導出されています。

結合楕円率#

結合 CP での 2 つの負の固有値(\(\lambda_1,\lambda_2\)、慣例により \(|\lambda_1|\ge|\lambda_2|\))は 大きさが等しい必要はありません — 実際の結合は完全に円柱対称ではありません — その比が結合楕円率 を定義します:

\[ \varepsilon = \frac{\lambda_1}{\lambda_2} - 1, \]

これは密度が一方の垂直方向で他方よりどれだけ多く枯渇しているかを測ります。完全に円柱対称な結合 (真の単結合、あるいは完全な軸対称を持つ三重結合)は \(\varepsilon\approx0\) を持ちます;優先的な \(\pi\) 蓄積方向を持つ二重結合は測定可能なほど大きな楕円率を示します。上記の水の O–H 結合 CP は \(\varepsilon\approx0.024\) を持ちます — 小さく、特別な \(\pi\) 性格のない通常のほぼ円柱対称な単結合と 整合します。

盆地積分:密度自身の分水嶺からの Bader 電荷#

臨界点だけでも結合を分類しますが、盆地積分はさらに進んで、(電荷分割解析の章 のどのスキームとも異なり)外部参照密度をまったく使わずにトポロジーを使って空間全体を分割します — 各原子核が \(\nabla\rho\) の勾配流のアトラクタであり、その盆地 \(\Omega_A\) は \(\rho\) における 最急上昇パスがそこで終わるすべての点の集合です。盆地はゼロフラックス面(面上のすべての点 \(\mathbf r\) について \(\nabla\rho(\mathbf r)\cdot\hat{\mathbf n}(\mathbf r)=0\))で境界づけられます — 勾配場が決して横切らない面であり、これがまさに、結果の原子占有 \(N_A=\int_{\Omega_A}\rho\,d\mathbf r\) を、なめらかな重み関数もどこにも恣意的な参照密度もない、真の曖昧さのない分割にするものです。 数値的には、立方実空間グリッド上で、qc-rs はグリッド上最急上昇(Henkelman 式アルゴリズム)に よってすべてのグリッド点を盆地に割り当てます:各点から、距離重み付きの密度増加する隣接点へ、局所 極大(原子核)に達するまで歩き、パス圧縮によって近くの点からの繰り返しの歩行を安価にします。

検証済み例 — 水の Bader(盆地積分)原子占有と電荷:

b = m.prop.qtaim.basin_integrate()
b["atoms"]
# O: 9.2826916251811 electrons -> q = 8 - 9.283 = -1.283
# H: 0.37835478988956356       -> q = 1 - 0.378 = +0.622
# H: 0.3626332357145994        -> q = 1 - 0.363 = +0.637

ここでの Bader 電荷(\(q_O\approx-1.28\))は、電荷分割解析の章 で 導出したどのスキームよりも顕著に絶対値が大きくなります — NPA/MBIS/Hirshfeld-I の \(\sim-0.87\)〜 \(-0.92\) よりも大きいです。これはバグでも外れ値でもありません:QTAIM のゼロフラックス分割は、 どんなストックホルダー重みや AO インデックススキームとも真に異なる、トポロジー駆動の基準で空間を 分けます、そして(その章自身の締めくくりが強調したように)収束すべき単一の「真の」原子電荷は 存在しません — QTAIM は単にもう 1 つの、独立に十分な根拠を持つ定義であり、他のものよりも「より 正しい」版ではありません。

ELF と LOL:臨界点検索なしの局在#

構造的に異なる実空間ツールが、「\(\rho\) 自体のトポロジー的特徴はどこか」ではなく「各点で電子対密度 はどれだけ空間的に局在しているか」を問います — Becke と Edgecombe による電子局在化関数(ELF) です。その構成はPauli 運動エネルギー密度から始まります、真の(正定値の、Kohn-Sham 軌道の) 運動エネルギー密度 \(\tau\) が、同じ密度の単一軌道(ボゾン的、完全に非局在化)記述が持つはずの von Weizsäcker 運動エネルギー \(\tau_W\) を超える超過分です:

\[ D(\mathbf r) = \tau(\mathbf r) - \tau_W(\mathbf r), \qquad \tau_W = \frac{|\nabla\rho|^2}{8\rho}. \]

\(D\) は、局所密度が真に単一軌道で記述可能などこでも(まったく Pauli 反発ペナルティがない — 理論的な 床)厳密にゼロになります、これが高い局在(共有結合、孤立電子対、内殻)の領域が小さな \(D\) を示す 理由です。これを有界で化学的に解釈可能な数値に変えるために、\(D\) は局所密度での一様電子ガスの 同じ量(Thomas-Fermi 運動エネルギー密度 \(D_h=\tfrac{3}{10}(3\pi^2)^{2/3}\rho^{5/3}\)、\(\rho\) だけの 純粋に局所的な関数)と比較されます:

\[ \chi(\mathbf r) = \frac{D(\mathbf r)}{D_h(\mathbf r)}, \qquad \text{ELF}(\mathbf r) = \frac{1}{1+\chi^2} \in(0,1]. \]

正規化は、単一軌道で真に記述される領域が厳密に \(D=0\Rightarrow\chi=0\Rightarrow\text{ELF}=1\) を 与える(単なる極限近似ではなく、qc-rs のテストスイートが直接検査する厳密な恒等式)ように意図的に 選ばれており、厳密に一様電子ガスの運動エネルギー密度を持つ領域は厳密に \(D=D_h\Rightarrow\chi=1\Rightarrow\text{ELF}=\tfrac12\) を与えます — この関数にその解釈尺度を与える 正準的な参照点です(ELF が \(\tfrac12\) を超えることは「この密度での一様電子ガスより局在している」を 意味し、\(\tfrac12\) を下回ることはそれより少ないことを意味します)。実際には、分子の外側境界を 超えると \(D\) は \(D_h\) より速く減衰することがあり、ELF を真空の裾で偽って 1 へ駆り立てます (「完全に局在した真空」、明らかに非物理的)。qc-rs は、化学的に興味深い結合/孤立電子対領域を 乱すことなくこの裾を手なずけるために特に \(D\) に小さな床を加える Multiwfn の規約に従います。

LOL(局在軌道位置指標、Schmider-Becke)は同じ2 つの材料、\(\tau\) と \(D_h\) から構築されますが、 異なる方法で組み合わされます — 正規化された二乗形式ではなく比です:

\[ t(\mathbf r) = \frac{D_h(\mathbf r)}{\tau(\mathbf r)}, \qquad \text{LOL}(\mathbf r) = \frac{t}{1+t} \in[0,1), \]

ELF と同じ定性的局在マッピングの目的を持つが数値尺度が異なり(実際には)内殻構造への視覚的強調が いくらか異なる有界関数を与えます — 2 つは競合する理論ではなく、密接に関連する背後の物理の相補的な 表現です。

関連する量、電子局在可能性指標(ELI-D)、\(\Upsilon_D=\rho_\sigma[12/g^\sigma]^{3/8}\) (\(g^\sigma=\rho_\sigma D_\sigma\)、ELF が使うのと同じ Pauli 運動エネルギー密度の分子)は、単一行列式 波動関数について ELF と厳密に同じトポロジーを共有します — 同じ位置で同じ殻・結合・孤立電子対を 示します — が、まったく一様電子ガス参照を必要とせず、その値は(恣意的な正規化点ではなく普遍 定数として一様ガスが \(\approx1.10\) を与える)絶対的な物理的意味を持ちます、ELF の相対的な、ガス 正規化された尺度とは異なります。

練習 17

  1. 完全な分子トポロジーに対する Poincaré-Hopf 和は厳密に \(+1\) でなければなりません。適度に大きな、 環を含む分子の数値臨界点検索が和 \(0\) を報告した場合、比較すべき参照計算を必要とせずに、ほぼ 確実に何が間違っているかを教えてくれますか。

  2. 水の酸素上の Bader 電荷(\(\approx-1.28\))は、その NPA 電荷(前章から \(\approx-0.91\))より顕著に絶対値が大きいです。この違いがどちらかの方法が単に間違っていることを 意味しない理由を一文で説明しなさい。

  3. ELF と ELI-D は単一行列式波動関数についてまったく同じトポロジー(同じ臨界点、同じ定性的な結合/ 孤立電子対の描像)を共有しますが、ELI-D は一様電子ガス参照をまったく必要とせず、ELF の構成全体は 1 つとの比較を中心に構築されています。これは、一様電子ガス比較が本質的な物理なのか、それとも 同じ背後の局在情報を捉えるための便利な正規化の選択にすぎないのか、について何を教えてくれ ますか。

トポロジーと局在は、収束済み密度から化学構造を直接読み取るための本マニュアルの実空間ツールキット を完結させます。反応性理論の章 は構造から反応性へ向かいます — 概念的 DFT 記述子(電気陰性度、硬度、Fukui 関数)と、分子が現在どう配置されているかだけでなく、 それがどうふるまうかを予測する芳香族性指標です。