ハミルトニアンの選び方で固体を考える 第2部 第1章 §2.1.4 ★★★ 電子どうしの斥力を平均に入れて独立に詰めると、結合の次数が結合の有無を決める。閉殻どうしは重なりのために反発する
独立な電子として詰める
電子が2個以上ある分子では、電子どうしの斥力があるので 2.1.3 節のようには解けない。第1部第1章 の1電子近似を使い、電子どうしの斥力を平均のポテンシャルに入れて、電子が互いに独立に結合軌道と反結合軌道に入るとみなす。Pauli の原理で各軌道にはスピンの向きの違う電子が2個まで入るので、電子の数 N N N に応じてH 2 + \mathrm{H}_2^+ H 2 + (N = 1 N=1 N = 1 )は結合軌道に1個、H 2 \mathrm{H}_2 H 2 (N = 2 N=2 N = 2 )は結合軌道に2個、H e 2 + \mathrm{He}_2^+ He 2 + (N = 3 N=3 N = 3 )はそれに反結合軌道の1個、H e 2 \mathrm{He}_2 He 2 (N = 4 N=4 N = 4 )は両方に2個ずつ入る。結合の次数 を
b = n + − n − 2 (2.1.15) b=\frac{n_+-n_-}{2}\tag{2.1.15} b = 2 n + − n − ( 2.1.15 )
で定義する。n + n_+ n + 、n − n_- n − は結合軌道と反結合軌道の電子の数である。H 2 + \mathrm{H}_2^+ H 2 + は b = 1 2 b=\frac12 b = 2 1 、H 2 \mathrm{H}_2 H 2 は b = 1 b=1 b = 1 、H e 2 + \mathrm{He}_2^+ He 2 + は b = 1 2 b=\frac12 b = 2 1 、H e 2 \mathrm{He}_2 He 2 は b = 0 b=0 b = 0 になる。
閉殻どうしは反発する
H e 2 \mathrm{He}_2 He 2 のように結合軌道と反結合軌道が2個ずつ埋まると、重なりが無ければ上がる分と下がる分が打ち消し合うはずである。ところが重なりがあると打ち消し合わない。原子が離れているときのエネルギー 4 H a a 4H_{aa} 4 H aa から測った電子のエネルギーの変化は、式 (2.1.11) の ε ′ \varepsilon ' ε ′ を使って
2 E + + 2 E − − 4 H a a = 4 ( ε ′ − H a a ) = 4 s ( s H a a + t ) 1 − s 2 = 4 s t ′ (2.1.16) 2E_++2E_--4H_{aa}=4(\varepsilon '-H_{aa})=\frac{4s(sH_{aa}+t)}{1-s^2}=4st'\tag{2.1.16} 2 E + + 2 E − − 4 H aa = 4 ( ε ′ − H aa ) = 1 − s 2 4 s ( s H aa + t ) = 4 s t ′ ( 2.1.16 )
と書ける。s > 0 s>0 s > 0 で、結合軌道が反結合軌道より低い限り t ′ > 0 t'>0 t ′ > 0 なので、これは正である。閉殻どうしの原子は近づくとエネルギーが上がる。2.1.2 節の直交化の言葉では、重なりのある軌道を直交化すると、各軌道に相手の軌道を引いた成分が混ざって運動エネルギーが上がり、サイトのエネルギーが ε ′ − H a a = s t ′ \varepsilon '-H_{aa}=st' ε ′ − H aa = s t ′ だけ上がる。2つの軌道を4個の電子で全部埋めると、結合と反結合の分裂 ± t ′ \pm t' ± t ′ は打ち消し合い、この上がった分だけが残る。希ガスの原子が分子を作らないこと、閉殻のイオンや分子が近づくと強く反発すること(固体ではイオン結晶の短距離の斥力)は、この式で説明できる。
H 2 \mathrm{H}_2 H 2 の分子軌道が削ったもの
H 2 \mathrm{H}_2 H 2 の2個の電子を両方とも結合軌道に入れた状態は、空間部分が
ψ + ( r 1 ) ψ + ( r 2 ) ∝ ϕ a ( r 1 ) ϕ b ( r 2 ) + ϕ b ( r 1 ) ϕ a ( r 2 ) ⏟ 各原子に1個ずつ + ϕ a ( r 1 ) ϕ a ( r 2 ) + ϕ b ( r 1 ) ϕ b ( r 2 ) ⏟ 片方の原子に2個 \psi_+(\bm{r}_1)\psi_+(\bm{r}_2)\propto\underbrace{\phi_a(\bm{r}_1)\phi_b(\bm{r}_2)+\phi_b(\bm{r}_1)\phi_a(\bm{r}_2)}_{\text{各原子に1個ずつ}}+\underbrace{\phi_a(\bm{r}_1)\phi_a(\bm{r}_2)+\phi_b(\bm{r}_1)\phi_b(\bm{r}_2)}_{\text{片方の原子に2個}} ψ + ( r 1 ) ψ + ( r 2 ) ∝ 各原子に 1 個ずつ ϕ a ( r 1 ) ϕ b ( r 2 ) + ϕ b ( r 1 ) ϕ a ( r 2 ) + 片方の原子に 2 個 ϕ a ( r 1 ) ϕ a ( r 2 ) + ϕ b ( r 1 ) ϕ b ( r 2 )
である(スピンは1重項)。重なりが小さければ、各原子に1個ずつの配置と片方に2個の配置が同じ重みで混ざる。R → ∞ R\to\infty R → ∞ でもこの重みは変わらないので、この状態は分子が解離した後も半分の確率で H + \mathrm{H}^+ H + と H − \mathrm{H}^- H − に分かれている。実際には2個の電子が同じ原子に入ると大きな斥力 U U U がかかるので、解離の極限は2つの中性の水素原子でなければならない。1電子近似は、この U U U を平均に入れたために、離れた原子の上で電子が互いを避けることを表せない。同じサイトの斥力 U U U を戻した2サイトの Hubbard 模型(第1部第1章 の総合問題)では、U / t U/t U / t が大きくなると片方に2個の配置の重みが 0 0 0 に近づき、各原子に1個ずつの配置(化学の Heitler–London の波動関数)に移る。
例題 2.1.4 電子を1個から4個まで詰める ★★
オリジナル 同じ原子2個からなる分子で、H a a = − 10 e V H_{aa}=-10\,\mathrm{eV} H aa = − 10 eV 、H a b = − t = − 4 e V H_{ab}=-t=-4\,\mathrm{eV} H ab = − t = − 4 eV 、s = 0.3 s=0.3 s = 0.3 とする。
(1) 結合軌道と反結合軌道のエネルギー E ± E_\pm E ± と、直交化した ε ′ \varepsilon ' ε ′ 、t ′ t' t ′ を求めよ。
(2) 電子を N = 1 , 2 , 3 , 4 N=1,2,3,4 N = 1 , 2 , 3 , 4 個詰めたとき、N H a a NH_{aa} N H aa から測った電子のエネルギーの変化を求め、結合の次数と比べよ。電子どうしの斥力の変化は無視する。
(3) N = 4 N=4 N = 4 の値が式 (2.1.16) の 4 s t ′ 4st' 4 s t ′ と一致することを確かめよ。
解答 (1) E + = ( − 10 − 4 ) / 1.3 = − 10.769 e V E_+=(-10-4)/1.3=-10.769\,\mathrm{eV} E + = ( − 10 − 4 ) /1.3 = − 10.769 eV 、E − = ( − 10 + 4 ) / 0.7 = − 8.571 e V E_-=(-10+4)/0.7=-8.571\,\mathrm{eV} E − = ( − 10 + 4 ) /0.7 = − 8.571 eV である。ε ′ = ( H a a + s t ) / ( 1 − s 2 ) = ( − 10 + 1.2 ) / 0.91 = − 9.670 e V \varepsilon '=(H_{aa}+st)/(1-s^2)=(-10+1.2)/0.91=-9.670\,\mathrm{eV} ε ′ = ( H aa + s t ) / ( 1 − s 2 ) = ( − 10 + 1.2 ) /0.91 = − 9.670 eV 、t ′ = ( t + s H a a ) / ( 1 − s 2 ) = ( 4 − 3 ) / 0.91 = 1.099 e V t'=(t+sH_{aa})/(1-s^2)=(4-3)/0.91=1.099\,\mathrm{eV} t ′ = ( t + s H aa ) / ( 1 − s 2 ) = ( 4 − 3 ) /0.91 = 1.099 eV である。H a a H_{aa} H aa から測ると結合軌道は 0.769 e V 0.769\,\mathrm{eV} 0.769 eV 下がり、反結合軌道は 1.429 e V 1.429\,\mathrm{eV} 1.429 eV 上がる。
(2) N = 1 N=1 N = 1 :− 0.769 e V -0.769\,\mathrm{eV} − 0.769 eV (b = 1 2 b=\frac12 b = 2 1 )、N = 2 N=2 N = 2 :− 1.538 e V -1.538\,\mathrm{eV} − 1.538 eV (b = 1 b=1 b = 1 )、N = 3 N=3 N = 3 :− 1.538 + 1.429 = − 0.110 e V -1.538+1.429=-0.110\,\mathrm{eV} − 1.538 + 1.429 = − 0.110 eV (b = 1 2 b=\frac12 b = 2 1 )、N = 4 N=4 N = 4 :− 1.538 + 2.857 = + 1.319 e V -1.538+2.857=+1.319\,\mathrm{eV} − 1.538 + 2.857 = + 1.319 eV (b = 0 b=0 b = 0 )である。結合の次数が同じ 1 2 \frac12 2 1 でも、N = 3 N=3 N = 3 は N = 1 N=1 N = 1 よりずっと弱く結合する。反結合軌道が上がる分の方が大きいためで、b = 0 b=0 b = 0 の N = 4 N=4 N = 4 では結合しないどころか反発する。
(3) 4 s t ′ = 4 × 0.3 × 1.099 = 1.319 e V 4st'=4\times0.3\times1.099=1.319\,\mathrm{eV} 4 s t ′ = 4 × 0.3 × 1.099 = 1.319 eV で、(2) の N = 4 N=4 N = 4 と一致する。
ほかの巻では
近似の数理 第5章 :Hartree–Fock が解離で誤る理由を相関エネルギーとして定義し、その大きさと限界を述べた所が読める
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com