量子化学
原田○ 講義概要
第1回
第2回
第3回
第4回
第5回
第6回
第7回

概論、量子化学の基礎演習1
分子の電子状態の計算法（Hückel法）
演習2
近似を高めた理論化学計算法
演習3
試験

準教科書・参考書
準教科書
・ 「入門 分子軌道法」藤永 茂 著（講談社サイエンティフィク、1990）
参考書
・「三訂 量子化学入門（上）」米澤･永田･加藤･今村･諸熊（化学同人、1983）
・「分子軌道法」藤永 茂（岩波、1980）
・「はじめて学ぶ量子化学」阿部正紀（培風館、1996）
・「量子化学演習」宮崎智雄（朝倉書店）
・「量子化学問題の解き方 付大学院入試問題」犬塚（東京化学同人、1981）
・「アトキンス物理化学（上）（下）」（第6版） P. W. Atkins著、千原・中村 訳、（東京化学同人、2001）
・M. J. Winter, “Chemical Bonding”, (Oxford, 1994) ISBN0-19-8556942

便利な概念が真実とは限らない！

例

作り物と本物（人工の概念と自然の真理）の区別が付かなくなる
O3 分子

(I)

自然なLewis構造

(II)

(III)

共鳴構造
（自然界には無い）

(IV)

実際の構造

はじめに
近代化学結合論にみる自然観
1．原子は、正電荷を持つ原子核（1中心）のまわりに負電荷を持った電子が群がったもの．分子は、正電荷を持つ原子核（多中心）と負電荷を持つ電子が集まったもの．
2. 自然界には４つの力．分子で重要なのはクーロン力のみ．
3. 電子は波であり粒子．
4. 原子核同士が電子を介してどのように集まれば安定になるのか、安定な集団（分子）はどのような性質を持つのか、どのようにすればうまく組み替えができるのか、即ち「構造・機能・物性・反応」が問題．

○ 量子力学/量子化学・統計力学・熱力学の相互関連性

熱力学

統計力学
ÍÎ 分配関数 ÍÎ 量子力学/量子化学
原子・分子分光学

量子力学/量子化学で計算（または分光的に測定）
Î 分子の並進・回転・振動・電子エネルギー準位
Î 分配関数（状態和：ある温度で分子が取る平均的な状態の数）
Î 全ての熱力学的関数（U, H, S, A, G, Kp, ・・・）
○理論化学の二大潮流
化学反応論（動的）と化学構造論（静的）
化学結合論は後者．伴によって立つのは量子化学．

【１】 量子化学の基礎◎ 到達目標：量子化学の概要および量子化学計算の本質を把握する．
◎ 歴史的外観： 量子力学から量子化学へ
Kekulé van't Hoff Rutherford Bohr Mendeleev Planck
Heisenberg de Broglie Schrödinger Pauli Heitler London
Pauling Mulliken Hückel Hartree Fock Slater

展開： スペクトルの説明、物理的測定の解釈、有機反応論、分子レベルの薬理・生体反応…

◎ 量子力学と量子化学
・量子力学
波動関数（＝状態関数, state function）Ψ
物理量を示す演算子αと物理量の観測値aの関係
波動方程式 αΨ＝aΨ系のエネルギーＥに関心が有れば αとしてＨを使うＨΨ＝ＥΨ （Schrödinger の波動方程式）
h2 2
H =−
∇ + V ( x, y , z )2m
主要な目的： 注目する系について方程式を解いて、Ｅを求め、Ψの形を知る

幾つかの約束事
１）粒子の存在確率：２）規格化：

Ψ＊Ψｄτ

Ψ＊Ψｄτ= 1３）Ψは、一価・連続・有限４）異なる固有値に対応する固有関数の直交性
５）物理量の期待値（平均値）： f = ∫ Ψ fΨdτ*６）重ね合わせの原理：Ψ = C A ΨA + C B ΨBcf. 直交関数系展開 i.e., Euler展開、Maclaurin
展開、Fourier変換

・Schrödinger の方程式
時間に依存しない方程式

Î

原子・分子の構造、物性

Î反応、ダイナミクス、光との相互作用HΦ = EΦ
時間に依存する方程式

h ∂Φ
= HΦ
−i ∂tcf.

Φ ∝ exp{− i(E / h )t}・量子化学ＨΨ＝ＥΨ を解く

（定常状態を扱う）

Ｈ は項数が膨大で、かつ、パラメータも含むが“既知”
一方、ΨとＥとは分からない．そこで、(1) 既知の関数,,･･･を重ね合わせる．
すなわち、全3次元空間に渡ってΨの形を仮定する．
ψ (1,2 ) = C1 χ a (1)χ b (2 ) + C 2 χ a (2 )χ b (1) +・・・
ex.
(2) エネルギー（固有値）を最適化（最小化）する．（摂動法や変分法）(3) 固有ベクトルΨの形を求める（C1,C2,・・・を決める）(4) 満足する．or より良い結果を得るために修正を続ける(繰り返し)(5) ΨとＥを使って、機能・物性・反応についての考察を進める注）χa(1), χb(2) ,･･･には、例えば水素類似原子の電子波動関数的なものを使う． （だから、水素のs, p, d,･･･軌道が重要）．

◎ 量子化学・理論化学の現状
・物質の系： Ｈ （既知）は項数が膨大でパラメータを含む複雑さ(Complexity in Chemistry)の中に、新機能が発現
cf. 物理の系： モデルポテンシャルで表される無限に広がる系、（ÍÎ 孤立系）
・Quantum Chemistry has come a long way!
化学の問題に量子力学が応用可能に！ （1970～．本当はここ10年弱）
最大の問題は計算時間！
原子数Ｎに対しての計算時間は、一般にＮ3～Ｎ4
Ｎ2.2（平尾）．Ｎ1（世界的目標． cf. 青木 elongation法）．
SCF法で100個、もっと精度を上げると数10個が限界（1999年12月）
励起状態、遷移状態などの計算はもっと時間が
しかし・・・、タンパク質などの巨大分子、溶媒中の溶質分子など、計算されてないものは無い．ただし、精度についてはある程度我慢
コンピュータ性能
CPU･MPUの高性能化、大型・並列化で処理速度は1.5～2倍/年で向上（10年で1000倍）
あと20年は続くだろう．2050年には脳の性能に達する？ でも、・・・・・理論化学者へ与えられたノーベル賞
Linus Carl Pauling (1954) Nature of the chemical bond ･･･
Robert S. Mulliken (1966) ･･･ electronic structure of molecules bond ･･･
Lars Onsager (1968) ･･･ reciprocal relations ･･･ thermodynamics ･･･
Ilya Prigogine (1977) Non-equilibrium thermodynamics ･･･
Kenichi Fukui and Roald Hoffmann (1981) ･･･ chemical reactions
Rudolph A. Marcus (1992) ･･･ electron transfer reactions ･･･
Walter Kohn and John A. Pople (1998) ･･･ density-functional theory（密度汎関数理論: ‘90～急増．今半分がこれ）･･･ computational methods in quantum chemistry(Ab-Initio, Gaussian70 ～)・ブレークスルーが求められる課題
メゾスコピック系（ナノ系）： N～100 Î N>1000．超分子、タンパク質
インターフェース系： 界面、溶媒
超高状態密度系： 多重励起
多原子系の量子効果： 量子ダイナミクス
・・・

h2 2
H =−
∇ + V ( x, y , z )2m◎ 水素類似原子： 1中心1電子系（静止核１個と電子１個）=== 軌道が重要 ===2()
Vr=−e/ (4πε 0 r )・特徴： 厳密解がある．

ϕ n ,l ,m (r ,θ , φ ) ∝= Rn ,l (r )Θ l ,m (θ ) exp(imφ )ｎｌｍ[ｓに対して、(E n = −me 4 / 8ε 0 h 2 n 22

主量子数（principal quantum number） 1,2,3,･･･
方位量子数 （azimuthal quantum number） 0,1,2,･･･,(n-1)
磁気量子数 （magnetic quantum number） 0,±1,±2,･･･,±l
スピン量子数（spin quantum number）]
cf. E1 = -13.6 (eV))

水素原子×1012
n=1
l=0
1s6×1012

n=2
l=02s2×1012xz

R(r)

R(r)

(b) l=143020.8

10.40

1 2 3
r[Å]

0400.40.2

r[Å]

0z
Pxz
Py

ψ8z

ϕxzdz

Pz

yx(c) l=2ϕyx4ϕyl=23d0.40.28ϕ

l=13p0l=12psy0.51ϕn=3l=03s1.05R(r)z(a) l=02ϕd x 2− y 212r[Å]Rπ .l (r)yyx

動径波動関数xz

位相（+,-）があることに注意

ϕzd xyyx

ϕzd zxϕd yzyxyx

波動関数（軌道）の角度依存h2 2H =−∇ + V ( x, y , z )2m◎ ヘリウム原子：===

1中心2電子系

三体問題：電子間反発が問題

===

⎛ h2
Ze 2 ⎞ ⎛ h 2
Ze 2 ⎞e2⎟⎟ + ⎜⎜ −⎟⎟ +H = ⎜⎜ −Δ1 −Δ2 −r1 ⎠ ⎝ 2mr2 ⎠ 4πε 0 r12⎝ 2m
注1) ハミルトニアンは式で表せる．解けないことが問題．
注2) 以降、簡素化のために原子単位系（=1, m =1, e =1）を使う．
なお、単位系の換算に頭を悩ます必要は全く無い．

⎛ 1Z⎞ ⎛ 1Z⎞ 11
H = ⎜⎜ − Δ1 − ⎟⎟ + ⎜⎜ − Δ 2 − ⎟⎟ += h1 + h2 +r1 ⎠ ⎝ 2r2 ⎠ r12r12⎝ 2
電子間反発項
水素原子１水素原子2「厳密には解けない」ことが問題！ Î 近似解法

ε=*Ψ
∫ HΨdτ*Ψ∫ Ψdτ

・摂動法
元の状態にわずかな変化が加わったときの結果を考える
cf. 振動電場中にある分子のハミルトニアン・変分法
変分原理： （真のE） < （近似波動関数で求めたE）
Î 安心して使える波動関数 をパラメータを含む適当な形に仮定．例えば、適切な関数の一次結合（linear combination）

Ψ = c1 χ 1 + c 2 χ 2 + ⋅ ⋅ ⋅ + c n χ nとして、エネルギーを最小にするようにパラメータを決めると、変分原理に基づき、最も確からしい解が得られる．
（行列の固有値・固有関数問題に帰着．後述）

He2+ - 2e
0-10エネルギー(eV)-20-3054.40eV-40
He ＋ e-50-60-70-80

24.58eV

20.40eV23.05eV

23.58eV

He-74.80eV
78.98eV-77.45eV -77.98eV

実測

変分法

摂動法

Hatree法

図 He原子のエネルギー

適切な関数 の一次結合（linear combination）

Ψ = c1 χ 1 + c 2 χ 2 + ⋅ ⋅ ⋅ + c n χ n

χi として代表的なのは、スレーター型軌道(STO)関数
原子波動関数に類似．
AOの計算では優れる．
MOの計算では長時間．
ガウス型軌道(GTO)関数
2電子積分が容易．
MOの計算で短時間．
例．STO-3G
3-21G図 STOとGTOの比較（各軌道は規格化されている）．

・Hartree法
１電子波動関数近似
波動関数 を１電子軌道 の積で表す．

Ψ = ϕ (r1 )ϕ (r2 )Schrödinger 方程式はZ⎡ 1⎤()
rϕ (r ) = εϕ (r ) ただし−Δ−+v
⎢ 2
⎥r
⎣
⎦

v(r ) = ∫ dv'

ϕ (r ')r − r'

SCF法（Self-consistent Field法 cf. 繰り返し計算）
適切な ϕ = ϕ(0) で ν = ν(0) を計算． ν(0) を使って方程式を解きϕ(1)を求める
引き続きを ϕ(1) 用いて ν(1) を計算． ν(1) を使ってを ϕ(2)求める
ϕ(n-1) ＝ ϕ(n) となるまで繰り返す

このHartree-Fock方程式の拡張がHartree-Fock-Roothaan方程式（後述）で、軌道計算の基礎式になる

22

◎ 原子の電子構造：

1中心N電子系===多体問題・ Pauliの原理： N電子系の波動関数は任意の2つの電子の座標交換に対して反対称（量子数の同じ状態に2個の電子が入ることはあり得ない．
かつ、電子には個別性がない）Ψ (ξ1 , ξ 2, ⋅ ⋅⋅, ξ i ,⋅ ⋅ ⋅, ξ k ,⋅ ⋅ ⋅, ξ N ) = − Ψ (ξ1 , ξ 2, ⋅ ⋅⋅, ξ k ,⋅ ⋅ ⋅, ξ i ,⋅ ⋅ ⋅, ξ N )Î

同じ空間軌道には、スピンの異なる2つの電子しか収容できない（Pauliの禁制）

１電子波動関数近似でこの制約をうまく満たす表現が、Slater行列式===

Slater行列式
１全電子波動関数を、直交化した分子軌道、ϕ1, ϕ2, ･･･の積（Hartree積）で表す： ψ1･ψ2･ψ3･ψ4･････ψN
ただし、ψ1=ϕ1(1)α(1), ψ2=ϕ1(2)β(2), ψ3=ϕ2(3)α(3), ････
α, β はスピン関数

電子の入れ換えに関してN!通り．２電子の交換に対して“反対称性”を示す表現は、

Ψ (1,2,⋅ ⋅ ⋅, N ) =略記法

ψ 1 (1) ψ 2 (1)1 ψ 1 (2) ψ 2 (2 )N!

MM

ψ 3 (1) L ψ N (1)ψ 3 (2) L ψ N (2 )MM

ψ 1 (N ) ψ 2 (N ) ψ 3 (N ) L ψ N (N )

Ψ (1,2,⋅ ⋅ ⋅, N ) = ψ 1 (1)ψ 2 (2) ⋅ ⋅ ⋅ ⋅β (n )“行列式は列の入れ替えで符号が変わる“ Î “反対称”

Slater行列式で表現された全電子波動関数を用いてHartree-FockRoothaan方程式を解く
Î お馴染み原子の電子配置を再現

・電子配置：
n=1,2,3,4, 5, 6, に K殻, L殻, M殻, N殻, O殻, P殻が対応電子の割り当て法： Pauli禁制則、Hundの規則、エネルギー
cf. 閉殻、最外殻・状態の表し方： ・状態は n 2S+1LJ で表す例えば、3P1/2 など．
・合成軌道角運動量Ｌ=0,1,2,3,… に、S, P, D, F, … が対応．
・合成スピン角運動量Sとして、スピン多重度2S+1
S=0, 1, 2, … に、1重項（singlet）, 2重項（doublet）が対応
・Jは軌道とスピンの合成角運動の量子数（ = L+S, L+S-1,･･･, L − S ）◎ 水素分子イオンの電子構造：===

ここから分子２中心1電子系

===
e・

1
1 1 1
H =− Δ− − +rA rB R2

分子軌道

HA－HB＋○ 厳密解（省略．詳細は参考書2に記載）○ 分子軌道(MO, molecular orbital)
・分子全体に広がった波動関数
・原子軌道(AO, atomic orbital)が無数にあったように、１つの分子についての分子軌道も無数にあり、異なる量子状態を持つそのうちの一つを図示したのが右図．

○ LCAO MO近似
分子軌道 ϕMO を原子軌道 ϕA と ϕB との重ね合わせで近似

ϕ MO = aϕ A + bϕ B
注）この場合水素の原子軌道．
一般論としては（n, l, m）で指定される軌道のどれでも良い
共鳴
H2+状態

HB状態

HA状態
erA

HA+HB+分子軌道

ϕ MOe-e-rB
HB+HA+

原子軌道

ϕ A ; n ,l , ml

HA－

HB+原子軌道ϕ B ; n ,l , ml

○ LCAO MO近似
分子軌道 ϕMO を原子軌道 ϕA と ϕB との重ね合わせで近似ϕ MO = aϕ A + bϕ B
変分原理に基づいてa, b を求めると

ϕ± =ϕA ±ϕB

2 (1 ± S )E± = En +ここで、
Q, クーロン積分（クーロンポテンシャルの期待値）**
Q = ∫ ϕ A V Aϕ A dτ = ∫ ϕ B V Bϕ B dτϕK, 共鳴積分（K < 0）
K = ∫ ϕ B V Aϕ A dτ = ∫ ϕ A V Bϕ B dτ**S, 重なり積分（S << 1）
S = ∫ ϕ A ϕ B dτ*Q±K
≅ En + Q ± K
1± S

SABを与える原子軌道の重なりϕBA

HA＋HB＋

○水素分子イオンの特徴+++

6E+

E-E |eV|420

ϕ_∞ϕA + ϕB
ϕBϕA

ϕ ∞ϕA − ϕB

ϕAHB+－200.100.20
R |nm|0.30

HA+0.40

HA+

ϕB(a )

水素分子イオン（1s状態）のポテ
ンシャル曲線HB+

結合性分子軌道ϕ 1( b)反結合性分子軌道ϕ

水素分子イオンの分子軌道（実線）と原子軌道（波線）・結合性電子と反結合性電子．結合性分子軌道と反結合性分子軌道
・結合エネルギーと結合距離・光励起に伴う解離機構などが、説明できる．

◎ 水素分子の化学結合：

２中心２電子系-er121r a2ra2-er b21111111 1−−−++
H = − Δ1 − Δ 2 −22ra1 ra 2 rb1 rb 2 r12 R1rb1+e+eabR・ LCAO MO近似＋１電子波動関数近似
これまでの手法の総合：
Slater行列式の軌道波動関数ϕi を適切な軌道の組{χμ }の線形結合

e-eϕAϕB
HBHAHA+HBmϕ i = ∑ χ μ C μiμ =1で表し、Hartree-Fock-Roothaan
方程式を解くH2e-ϕi (MO)1sσ *E－K1s1sQElEl－K1sσEHAH2

2つの水素原子が結合して水素分子を形成したときの電子エネルギー準位とその占有状態

HB・残る問題
1) 現実のH2の中では、2つの電子は互いにうまく避け合うように運動するはず．ここまでのMO近似にはこの効果（電子間相関, electroncorrelation）が含まれない．これは、平均場を使っているから．
2) H2の波動関数が核間距離Rが大きいところで、H+Hにならない．
これらは、配置間相互作用（configuration interaction, CI）を取り入れると改善される．簡単には、別の電子配置の軌道も考慮に入れて精度を上げる（準テキスト186～192頁参照）．

・別アプローチ： Heitler-London の取り扱い （基底状態の水素分子）分子の状態とは、２つの原子に属する２つの価電子が、原子に属している状態を保ちながら相互作用し、交換し合っているとする
cf. 核 a, b； 電子 1, 2１）水素原子の1s状態： χa(1), χb(2)ψ (1,2) = C1 χ a (1)χ b (2) + C 2 χ a (2)χ b (1)から、一次結合を作る (共有結合)
C1 = ±C 2
等核だから、Î ψ 1 = N1 {χ a (1)χ b (2) + χ a (2)χ b (1)} 、 ψ 3 = N 3 {χ a (1)χ b (2) − χ a (2)χ b (1)}
規格化すると、 N1 = 1 / 2 + 2S 2
N 3 = 1 / 2 − 2S 2
（Ｓ： 重なり積分overlap integral）
注）２つの状態を組合わせて別の２つの状態を作る．
（状態の数は変わらない）
２）スピンを考えた全波動関数にPauliの原理（電子の交換に対して、反対称＝符号が変わる）を適用．

３）エネルギーを求めると
E 1 = 2EH +J+K1+ S 2

E 3 = 2EH +J −K
1− S 2

Ｊ： クーロン積分Coulomb integral
Ｋ：交換積分exchange integra

2 + 2S2{χ a (1)χ a (2) + χ b (1)χ b (2)}

考えて、さらに、
ψ = a cψ cov + aiψ ion

21Re 10-12

3.14eV-2-3に変分法を適用.
エネルギーを求める．

1− S

谷の深さ

ψ ion =1エネルギー(eV)４）イオン結合状態を表す関数J −K2

核間距離(Å)

J +K1+ S2

◎ 要するに、・H の形が分かっても全波動関数は複雑で分からない（エネルギーも計算できない）
・そこで、最もその状態を反映しそうな関数を幾つか選び出し、状態の線形結合を作って（重ね合わせて）、変分法でエネルギーの最小値を求める
・最小エネルギーに対応する波動関数は、真に近いとして満足する
・その関数を使って分子の性質を予想し実験値と比較、良く一致するように改良する
なお、励起状態を考える際には、励起状態を表すのに相応しい関数を用いる

○ 原子価結合法と分子軌道法===

関数近似法：まとめ===・原子価結合法（Valence bond 法、VB法． cf. 前述の水素分子の取扱法）
Heitler, London 共有結合の理論的解明、交換エネルギー、結合の飽和性(1927)
Pauling混成軌道(1931)、共鳴理論(1933)
・分子軌道法（MO法）
Mulliken(1927), Hund(1928)
原子価結合法 ψ (1,2) = C1 χ a (1)χ b (2) + C 2 χ a (2)χ b (1)①ψ = C a χ a + Cb χ b
分子軌道法②
原子価結合法： 化学者に馴染み易い．理解し易い．直感に訴えやすい．
長時間
分子軌道法： 計算しないと理解し難いが、多くの理解が得られる．
短時間
大雑把な話：

◎ ２原子分子の電子状態

孤立電子○ σ 結合とπ 結合
σ 結合： 軌道の重なり//結合軸
sdz 2π 結合： 軌道の重なり⊥結合軸
pxdzr=== 2核多電子系（復習+α） ===

分子

融合原子

S=0

○ 結合の強さ：
重なり積分Sが大きいほ
ど強く、S = 0なら結合は生じない．

sspz

pxpxdz2sdz2pz

dxz

pz

pypxdz2pxpydz2px

dx2-y2pxspz

px

dx2-y2S=0結合軸z軸dxzdxzspy○ 等核2原子分子○ 異核2原子分子AB

2pxσ *

ϕA
2p2 p yπ * 2 p zπ *ϕ_ = cAϕA − cBϕB

2p2pxσ2 p yπ2s1sϕB

2 p zπ2sσ *2sσ

2s1sσ *
1sσ

1s(a) B2,C2,N2ϕl = cAϕA + cBϕBAABB◎混成軌道=== 多核多電子系（復習+α） ===◎Cの混成軌道sp3

2p2p混成励起

基底状態(2価)C2H4sp2(昇位)2sCH4(4価)2pspC2H2