概念 DFT と反応性の理論

概念 DFT と反応性の理論#

これまでの Part II の各章はすべて「固定された電子数の固定された分子について、エネルギーは 何で、密度はどう構造化されているか」を問うてきました。本章は異なる種類の問いを発します:電子数 自体が変化したらエネルギーはどう変わるか、そして分子中のどの原子が求核剤・求電子剤・ラジカルに 対して最も反応性があるか。概念 DFTは、\(N\) に関する密度汎関数摂動論から直接導出される一握りの 大域的スカラー記述子で最初の問いに答えます;Fukui 関数は、密度自体がその同じ摂動に局所的に どう応答するかを問うことで 2 番目に答えます。本章はまた、電子非局在化に基づく芳香族性指標族 — PDI、FLU、\(I_{\text{ring}}\)、MCI — を導出します。これは、ガイドの章で既に導入 された純粋に幾何学的な HOMA 指標を、真に異なる、波動関数に基づいて補完するものです。 crates/qc-prop/src/{cdft,fukui,aromaticity,homa,bird}.rs に直接基づきます。

概念 DFT:\(\partial E/\partial N\) からの反応性#

Kohn-Sham DFT の基底状態エネルギーは、形式的には外部ポテンシャルだけでなく電子数 \(N\) の汎関数でも あります — そして固定幾何での \(N\) に関するその微分は、まさに化学的に馴染みのある電気陰性度と硬度の 概念であり、単に定性的なものではなく数学的に精密になったものです。ここですべての土台となる 2 つの エネルギーは垂直イオン化ポテンシャル \(I\) と垂直電子親和力 \(A\) です — カチオンとアニオンの エネルギーを中性分子自身の幾何で評価したもの(それら自身の緩和した幾何ではない)であり、2 通りで 得られます:Koopmans の定理(\(I\approx-\varepsilon_{\text{HOMO}}\)、 \(A\approx-\varepsilon_{\text{LUMO}}\)、中性 SCF 自身の軌道エネルギーから無料で、追加計算不要)または \(\Delta\)SCF(\(I=E(N-1)-E(N)\)、\(A=E(N)-E(N+1)\)、2 回の追加の一点エネルギー、より正確だが実際の 電子数摂動が必要)。

\((I,A)\) だけから、\(\partial E/\partial N\) と \(\partial^2E/\partial N^2\) への有限差分近似が化学 ポテンシャルと硬度を直接与えます:

\[ \mu = -\frac{I+A}{2}, \qquad \chi = -\mu = \frac{I+A}{2}\ \text{(Mulliken 電気陰性度)}, \qquad \eta = I-A\ \Bigl(=\frac{\partial^2E}{\partial N^2}\Bigr), \qquad S=\frac1\eta. \]

\(\mu\) は文字通り、中性分子での \(E(N)\) の(有限差分)傾きです — 電子を 1 つ加えるか除くかにかかる エネルギーコストであり、平衡において分子全体で均等化されています、これはまさに「電気陰性度」が 捉えるべき化学的意味です:非常に負の \(\mu\)(大きな \(\chi\))を持つ物質は電子密度を自分の方へ引き 寄せます。曲率 \(\eta\) はその電子移動プロセスへの抵抗を測ります — 大きな \(\eta\)(「硬い」種)は 電子密度を得ることも失うことも抵抗し、小さな \(\eta\)(「柔らかい」種、大きな \(S=1/\eta\))は分極 しやすく反応性があります。

\(\mu\) と \(\eta\) を組み合わせると、電気陰性度/硬度単独よりもさらに具体的な反応性の問いに答える いくつかの高次記述子が得られます:

\[ \omega = \frac{\mu^2}{2\eta}\ \text{(求電子性指標、Parr-Szentpaly-Liu)}, \qquad \Delta N_{\max} = -\frac{\mu}{\eta}\ \text{(完全な電子貯蔵庫への最大電子流)}, \]
\[ \omega^- = \frac{(3I+A)^2}{16(I-A)}\ \text{(電子供与力)}, \qquad \omega^+ = \frac{(I+3A)^2}{16(I-A)}\ \text{(電子受容力、Gazquez-Cedillo-Vela)}, \]
\[ \Delta_{\text{nu}} = \frac{(\mu+\eta)^2}{2\eta}\ \text{(求核脱離能)}, \qquad \Delta_{\text{el}} = \frac{(\mu-\eta)^2}{2\eta}\ \text{(求電子脱離能、Ayers-Anderson-Bartolotti)}, \qquad \Delta\omega^\pm = \omega^++\omega^-. \]

求電子性指標 \(\omega\) はこれらの中で最も広く使われるものです — それは、完全な電子貯蔵庫が供給 するであろう最大の電荷流を受け入れることでその種が安定化するエネルギーを測るので、大きな \(\omega\) は、その挙動を \(I\) か \(A\) のどちらが個別に駆動しているかによらず良い求電子剤を識別します。 求核脱離能と求電子脱離能は関連するが異なる問い — 「この分子全体はどれだけ反応性があるか」ではなく 「この種は、他所で結合した後、アニオンとして(求核脱離能)あるいはカチオンとして(求電子脱離能) どれだけ容易に脱離するか」— を定量化します、攻撃する求核剤/求電子剤の問いではなく脱離基の問いです。

検証済み例 — 水の大域反応性記述子(\(\Delta\)SCF 品質の \(I,A\) から、eV 単位):

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()
r = m.prop.cdft.reactivity()
r["ip"], r["ea"]                    # (13.554, -4.795)  eV
r["electronegativity"], r["hardness"]  # (4.379, 18.349)  eV
r["electrophilicity"]               # 0.523  eV -- 小さい:水は貧弱な求電子剤

水の大きな硬度(\(\eta\approx18.3\ \text{eV}\))と小さな求電子性(\(\omega\approx0.52\ \text{eV}\))は、 大きな HOMO-LUMO ギャップを持つ小さく閉殻で運動学的に不活性な分子について化学的直感が予測する まさにその通りです — それは電子密度を得ることも失うことも強く抵抗し、対応して貯蔵庫からの電荷を 受け入れることで安定化する傾向がほとんどありません。

Fukui 関数:大域反応性が局所的にどこに現れるか#

大域記述子は「この分子全体はどれだけ反応性があるか」に答えますが、「どの原子が最初に反応するか」 には答えません。Fukui 関数 \(f(\mathbf r)=\partial\rho(\mathbf r)/\partial N\) はまさにそれに 答えます — 特定の点での密度が電子の追加や除去にどう応答するか — そしてその原子への凝縮形は、 分子の異なる部位にわたって直接比較可能な原子ごとの反応性マップを与えます:

\[ f_A^+ = q_A(N) - q_A(N+1)\ \text{(求核攻撃部位 — 電子密度を得ることを好む)}, \qquad f_A^- = q_A(N-1) - q_A(N)\ \text{(求電子攻撃部位)}, \]
\[ f_A^0 = \tfrac12(f_A^++f_A^-)\ \text{(ラジカル攻撃部位)}, \qquad \Delta f_A = f_A^+-f_A^-\ \text{(二重記述子 — 符号が支配的な性格を識別する)}. \]

qc-rs はこれらの**凍結分子軌道(FMO)**近似を実装します — パラメータフリーで本質的に無料で計算 できる近道であり、中性分子の収束済み波動関数だけが必要で、追加の \(N\pm1\) 一点計算はまったく 不要です。この近似は、真の密度差 Fukui 関数をフロンティア軌道密度で直接置き換えます、 \(f^+(\mathbf r)\approx|\varphi_{\text{LUMO}}(\mathbf r)|^2\) および \(f^-(\mathbf r)\approx|\varphi_{\text{HOMO}}(\mathbf r)|^2\) — 物理的には、「次の電子はどこへ行く か」が「LUMO が既に振幅を持つ場所」で近似され、除去については対称的に HOMO で近似されます。 Mulliken 式軌道組成を介して原子へ凝縮すると、これはまさにその組成行列のフロンティア行が そのものになります:

\[ f_A^+ = \text{comp}_{\text{LUMO},A}, \qquad f_A^- = \text{comp}_{\text{HOMO},A}, \]

そして Mulliken 軌道組成行列のすべての行が構成上厳密に 1 に合計するので、和則 \(\sum_Af_A^+=\sum_Af_A^-=1\) と \(\sum_A\Delta f_A=0\) は近似的にではなく厳密に成り立ちます — QTAIM の Poincaré-Hopf 恒等式 と同じ精神での、参照計算を 必要としない真に無料の正しさ検査です。(より正確だがより高価な有限差分 Fukui 変種 — 凍結 フロンティア軌道で近似するのではなく実際の \(N,N\pm1\) SCF 原子電荷を差分する — は自然な将来の 追加です;FMO 形が qc-rs が今日実装するものです。)

検証済み例 — 水の凝縮 Fukui 関数:

f = m.prop.cdft.fukui()
f["fukui_minus"]      # [0.665, 0.167, 0.167]  -- HOMO 組成:O が支配的(求電子攻撃部位)
f["fukui_plus"]       # [0.237, 0.381, 0.381]  -- LUMO 組成:H が支配的(求核攻撃部位)
f["dual_descriptor"]  # [-0.428, 0.214, 0.214] -- O で負(求電子性)、H で正(求核性)

水の酸素は HOMO を支配します(その孤立電子対は最高エネルギーの占有密度であり、まさに求電子剤が 攻撃するであろう場所です)一方、水素は LUMO を支配します(求核剤の電子対が流れ込むであろう \(\sigma^*\) 性格)— O(負、求電子性)と H(正、求核性)の間の二重記述子の符号反転が、この定性的な 描像を原子ごとの単一の曖昧さのない数値にします。

幾何を超えた芳香族性:電子非局在化指標#

ガイドの芳香族性の章 は既にHOMA(調和振動子芳香族性 モデル)を導入しました、各環結合の長さが表になった「理想的な芳香族」参照長にどれだけ近いかから 構築される純粋に幾何学的な指標です。HOMA の構造的限界はまさにそれが純粋に幾何学的であることです — それは、幾何がその代理であるはずの実際の電子非局在化について何も述べません。電子非局在化 芳香族性族は代わりに波動関数から直接非局在化を測ります、同じ非局在化指標 \(\delta(A,B)\) (原子 \(A,B\) 間の共有電子対の尺度、Mayer 結合次数 と精神的に 密接に関連)と占有 MO 基底での原子重なり行列 \(S(A)\) に基づきます:

\[ \text{PDI} = \text{6 員環の 3 つの}\textit{para}\text{ }\delta\text{の平均}\ \text{(高いほど芳香族性)}, \qquad \text{FLU} = \Bigl\langle\Bigl(\frac{\delta_i-\delta_i^{\text{ref}}}{\delta_i^{\text{ref}}}\Bigr)^2 \Bigr\rangle\ \text{(0 = 芳香族)}, \]
\[ I_{\text{ring}} = 2^n\operatorname{Tr}\bigl[S(A_1)S(A_2)\cdots S(A_n)\bigr]\ \text{(1 つの環パスに 沿った真の}n\text{中心非局在化)}, \qquad \text{MCI} = \text{すべての環順列にわたる}I_{\text{ring}}\text{の平均}. \]

\(I_{\text{ring}}\) はここで概念的に中心となる対象です — 2 中心非局在化指標の真の \(n\) 中心一般化 (\(I_{\text{ring}}\) の 2 中心の場合は \(\delta(A,B)\) 自体に厳密に帰着します)であり、隣接原子間で ペアごとにだけではなく環パス全体にわたって一度にどれだけの電子密度が非局在化しているかを測り ます。\(I_{\text{ring}}\) の値は環をどちらの方向にたどるか、どこから始めるかに依存するので、 MCI(多中心指標)は環原子のあらゆる異なる順列にわたってそれを平均し、単一の順列不変な数値を 与えます — その不変性の代償は、より大きな環についてより組み合わせ的に多くの環パス評価が必要な ことです。PDI と FLU はより安価で、より的を絞った量であり、完全な多中心トレースではなくペアごとの \(\delta\) 値から直接構築されます:PDI は特に 6 員環に固有のパラ関係(環をまたぐ反対側の原子)を 利用し、FLU は同じ理論レベルの参照芳香族系に対して非局在化がどれだけ均一に環に広がっているかを 測ります(FLU は呼び出し側が供給する外部参照 \(\delta^{\text{ref}}\) を必要とします、「どれだけの 非局在化が正常か」は本質的に比較であり絶対尺度ではないからです)。

検証済み例 — ベンゼン、HOMA/Bird(幾何学的)と電子非局在化指標:

benzene = """
C  0.0000  1.3970  0.0000
C  1.2098  0.6985  0.0000
C  1.2098 -0.6985  0.0000
C  0.0000 -1.3970  0.0000
C -1.2098 -0.6985  0.0000
C -1.2098  0.6985  0.0000
H  0.0000  2.4810  0.0000
H  2.1486  1.2405  0.0000
H  2.1486 -1.2405  0.0000
H  0.0000 -2.4810  0.0000
H -2.1486 -1.2405  0.0000
H -2.1486  1.2405  0.0000
"""
mb = qc.chk.new(atom=benzene, ao="sto-3g", unit="angstrom").scf(ref="r").run()
ring = [0, 1, 2, 3, 4, 5]
mb.prop.arom.homa(ring=ring)      # 0.9792  -- 1 に近い:結合が理想芳香族長に近い
mb.prop.arom.bird(ring=ring)      # 99.987  -- 100 に近い:結合次数が本質的に一様
mb.prop.arom.indices(ring=ring)
# {'pdi': 0.1029, 'flu': 0.000259, 'iring': 0.0434, 'mci': 0.0662, ...}

すべての指標が独立にベンゼンの教科書的な芳香族性を確認します — HOMA と Bird はそれぞれの理想値 (1 と 100)に近く、PDI は非局在化指標の文献がベンゼンについて具体的に報告する参照「芳香族」値 \(\approx0.1\) に近く、FLU は本質的にゼロです(一様な非局在化、有意な結合交替駆動のずれなし)。 幾何学的(HOMA/Bird)族と電子非局在化(PDI/FLU/\(I_{\text{ring}}\)/MCI)族は関連するが真に異なる 問いに答えています — 一方は結合長/結合次数の一様性について、他方は波動関数自身の多中心電子共有 について — そしてここでのそれらの同時の一致は、同じ計算の冗長な繰り返しではなく、強力で複数手法 にわたる確認です。

練習 18

  1. 水の求電子性指標 \(\omega\approx0.52\ \text{eV}\) は小さく、その硬度 \(\eta\approx18.3\ \text{eV}\) は大きいです。\(\omega=\mu^2/(2\eta)\) の分母にある大きな \(\eta\) が、\(\mu\) が何であれ、硬く 反応性のない閉殻分子について小さな求電子性を生むとまさに期待される理由を一文で説明しなさい。

  2. Fukui 関数への FMO 近似は、より正確な有限差分変種とは異なり、まったく追加の \(N\pm1\) 計算を 必要としません。真の \(N\pm1\) 計算なら必要としないであろう、FMO が行う 1 つの具体的な物理的 仮定は何ですか。そしてどんな状況で(考えてみてください:ほぼ縮退したフロンティア軌道、あるいは イオン化時に強い軌道緩和を持つ系)FMO が著しく信頼できなくなると期待されますか。

  3. PDI は(パラ関係を介して)6 員環についてのみ定義されますが、\(I_{\text{ring}}\)/MCI にはそう した制限がありません。なぜパラの概念が具体的に 6 員環を要求するのか、そしてこれが \(I_{\text{ring}}\)/MCI がこの指標族のより一般的に適用可能なメンバーである理由について何を 教えてくれますか。

本章は、収束済みの密度と波動関数が何を教えてくれるかについての Part II の巡回を完結させます: エネルギーと構造(Part I-II の基礎の章)、力と振動(解析的微分、 Hessian)、環境(ECP、溶媒和、 分散)、そして今や電荷分布・トポロジー・反応性です。 マルチリファレンス展望 は、qc-rs がまだ実装していないもの — CASSCF ベースの相関法 — とその理由を見ることで Part II を締めくくります。