芳香族性指標

芳香族性指標#

芳香族性 — ベンゼンのような共役環の特別な安定性と非局在化 — は化学で最も有用な概念の 1 つで、原子電荷と同様、 単一の定義がないものです。幾何（等しくなった結合長）にも、電子の非局在化にも、磁気応答にも現れるので、指標 ごとに異なる側面を探ります。qc-rs はそれらを qc.prop.arom にまとめています;本章は、参照芳香族としてベンゼンを 使い、計算と比較の方法を示します。

理論:芳香族性への 3 つの窓#

芳香環は非芳香環と測定可能な形で異なり、指標のファミリごとに異なる窓から見ます:

  • 幾何的 — 芳香環は結合長が等しく（局在した単/二重の交替がない）なります。HOMA（芳香族性の調和振動子 モデル）がこれを採点します:完全に芳香族な環で HOMA ≈ 1、結合が局在するにつれ 0（や負）へ下がります。 BIRD は結合次数ベースの関連指標です。

  • 電子的 / 非局在化 — 芳香族性は電子が環の周りで共有されることを意味します。FLU（芳香族揺らぎ指数、 芳香族で ≈ 0）、PDI（パラ非局在化指数）、MCI（多中心指数）、I_ring が、密度行列からこの非局在化 を測ります。

  • 情報理論的 / その他 — Shannon エントロピーベースおよび関連記述子（shannon, info_theoretic, lolipop, …）。

異なるものを測るので、単一の数を信じるのではなく、複数を報告して一致を見ます。

使い方#

指標は環 — 環を構成する原子番号のリスト、順序どおり — を取ります。ベンゼン（下の幾何で炭素が番号 0, 2, 4, 6, 8, 10）では:

import qc, math
# D6h ベンゼンを作る
R, Rh = 1.39, 1.39 + 1.08
atoms = []
for i in range(6):
    a = math.radians(60 * i)
    atoms.append(f"C {R*math.cos(a):.4f} {R*math.sin(a):.4f} 0")
    atoms.append(f"H {Rh*math.cos(a):.4f} {Rh*math.sin(a):.4f} 0")
benzene = ";".join(atoms)

m = qc.chk.new(atom=benzene, ao="6-31g", unit="angstrom").scf(ref="r").run()   # E = -230.624263
ring = [0, 2, 4, 6, 8, 10]   # 6 個の炭素

qc.prop.arom.homa(m, ring=ring)      # 0.999    幾何的 — ほぼ完璧
qc.prop.arom.bird(m, ring=ring)      # 99.99    結合次数の幾何的指標（~100 = 芳香族）
qc.prop.arom.indices(m, ring=ring)   # {'pdi': 0.093, 'flu': 0.001, 'iring': 0.034, 'mci': 0.051, ...}

ベンゼンはすべての指標で教科書的な芳香族と採点されます:HOMA 0.999（≈ 1）、BIRD 99.99（≈ 100）、そして ほぼゼロの FLU 0.0006 — 3 つの窓すべてが環が完全に芳香族だと一致します。indices leaf は非局在化尺度（PDI, FLU, I_ring, MCI、および大きな環用の AV1245/AVmin 記述子）を 1 回の呼び出しでまとめます。

Tip

環を正しく指定する ring= リストは環の原子を、環に沿った接続順（恣意的でなく）で与えねばなりません。順序を誤ると、環を一周 することに依存する多中心指標（MCI, I_ring）が乱れます。幾何的な HOMA はより寛容です。

環を比較する#

芳香族性指標は相対的に最も有用です — この環はあの環より芳香族か？各環で同じ指標を計算して比較します:

# 概略:ベンゼン環を非芳香環と比較
homa_benzene = qc.prop.arom.homa(m, ring=ring)      # ~0.999（芳香族）
# シクロヘキサンの HOMA は 1 をはるかに下回る（局在した単結合）

HOMA が 1 に近く FLU が 0 に近ければ芳香環;HOMA が 1 をはるかに下回る（や負）で FLU が大きければ、局在した非芳香 （や反芳香）環を示します。

総合例#

import qc, math
R, Rh = 1.39, 2.47
atoms = []
for i in range(6):
    a = math.radians(60 * i)
    atoms += [f"C {R*math.cos(a):.4f} {R*math.sin(a):.4f} 0",
              f"H {Rh*math.cos(a):.4f} {Rh*math.sin(a):.4f} 0"]
m = qc.chk.new(atom=";".join(atoms), ao="6-31g", unit="angstrom").scf(ref="r").run()
ring = [0, 2, 4, 6, 8, 10]

print("HOMA :", round(qc.prop.arom.homa(m, ring=ring), 3))          # 0.999
print("BIRD :", round(qc.prop.arom.bird(m, ring=ring), 2))          # 99.99
print("FLU  :", round(qc.prop.arom.indices(m, ring=ring)["flu"], 4)) # 0.0006

練習 13

  1. ある環が HOMA = 0.30、FLU = 0.05 を与えます。芳香族か否か。2 つの指標は一致しますか、そして低い HOMA は物理的に 何を意味しますか。

  2. 1 つだけでなく、幾何的指標（HOMA）と非局在化指標（FLU/MCI）の両方を報告するのが良い作法なのはなぜですか。

  3. 芳香族と分かっている環で arom.indices の MCI 値がおかしく見えます。最も可能性の高い入力ミスは何ですか。

最後の物性章、スペクトルと DOS は、結合解析から軌道エネルギーのスペクトル — 状態密度と HOMO–LUMO ギャップ — へ移ります。