tecktom.com

この章のまとめ

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第1章 原子2個:結合軌道と反結合軌道

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第2部 第1章 §2.1.7 ★★★ 2つの軌道を残すと重なりつきの 2×22\times2 行列になり、結合軌道と反結合軌道の分裂がホッピングの大きさを与える。これを NN 個に広げると帯になる

H2+\mathrm{H}_2^+ の電子の状態を、2つの陽子のまわりの 1s1s 軌道が張る部分空間に限ると、Schrödinger 方程式は一般化固有値問題 Hc=ESc\mathsf H\bm{c}=E\mathsf S\bm{c} になる。反転対称性で固有ベクトルは ϕa±ϕb\phi_a\pm\phi_b に決まり、準位は結合軌道 E+=(Haa−t)/(1+s)E_+=(H_{aa}-t)/(1+s) と反結合軌道 E−=(Haa+t)/(1−s)E_-=(H_{aa}+t)/(1-s) に分かれる。結合軌道は陽子の間に電子を集めてエネルギーを下げ、Re=2.49a0R_{\mathrm e}=2.49a_0 に谷を作る。LCAO の値は Rayleigh–Ritz の上界で、厳密解(Re=2.00a0R_{\mathrm e}=2.00a_0、De=2.79 eVD_{\mathrm e}=2.79\,\mathrm{eV})より結合を弱く見積もる。

重なりを直交化で取り込むと、行列は第1部第1章の tight-binding の形 ε′∓t′\varepsilon '\mp t' に戻る。サイトのエネルギーは重なりの分 st′st' だけ上がるので、閉殻どうしの原子は近づくと反発し、He2\mathrm{He}_2 は結合しない。電子を独立に詰めたときの結合の次数 (n+−n−)/2(n_+-n_-)/2 が結合の有無の目安になる。2つの軌道だけを残す近似は、井戸が深く離れた極限で正しくなり、そのときホッピングは井戸の間のトンネルの振幅 e−κ0d\mathrm{e}^{-\kappa_0d} に比例する。

次の章(第2部第2章)では、原子が3個以上つながった分子の π 電子に同じ方法を当て、Hückel 法で N×NN\times N の行列を対角化する。第2部第3章で原子を輪にして NN を大きくすると、結合軌道と反結合軌道の2準位が幅 4t′4t' の帯に詰まる。

この章の結果

  • 2つの原子軌道で Hc=ESc\mathsf H\bm{c}=E\mathsf S\bm{c}。同じ原子なら E±=(Haa∓t)/(1±s)E_\pm=(H_{aa}\mp t)/(1\pm s)、Hab=−tH_{ab}=-t。
  • H2+\mathrm{H}_2^+:s=e−R(1+R+R2/3)s=\mathrm{e}^{-R}(1+R+R^2/3)、C=1/R−e−2R(1+1/R)C=1/R-\mathrm{e}^{-2R}(1+1/R)、X=e−R(1+R)X=\mathrm{e}^{-R}(1+R)、E±(R)=εa+1/R−(C±X)/(1±s)E_\pm(R)=\varepsilon _{\mathrm a}+1/R-(C\pm X)/(1\pm s)。
  • 直交化すると ε′=(Haa+st)/(1−s2)\varepsilon '=(H_{aa}+st)/(1-s^2)、t′=(t+sHaa)/(1−s2)t'=(t+sH_{aa})/(1-s^2)。tight-binding の tt と Hückel の −β-\beta はこの t′t'。
  • 閉殻どうしの反発は 4st′>04st'>0。結合の次数 b=(n+−n−)/2b=(n_+-n_-)/2。
  • LCAO の固有値は厳密な値の上界(A 型)。深い井戸の極限で分裂は e−κ0d\mathrm{e}^{-\kappa_0d} に比例し、指数は正しいが前の係数は原子軌道のままでは合わないことがある(B 型)。

総合問題 水素分子イオンの結合の長さ・結合エネルギー・振動の量子 ★★★

オリジナルH2+\mathrm{H}_2^+ の全エネルギーを、LCAO の式 (2.1.14) の E+(R)E_+(R) で近似する。比べる厳密な値として、同じハミルトニアンの数値解のRe=1.997a0R_{\mathrm e}=1.997a_0、De=0.1026EhD_{\mathrm e}=0.1026E_{\mathrm h}、谷の曲率 d2E/dR2=0.1030Eh/a02\mathrm{d}^2E/\mathrm{d} R^2=0.1030E_{\mathrm h}/a_0^2 を使う。陽子と電子の質量比は Mp/m=1836.15M_{\mathrm p}/m=1836.15 とする。

  1. (1) (使う道具:2.1.1節・2.1.3節)E+(R)E_+(R) を数値で最小にして、LCAO の ReR_{\mathrm e} と De=−12−E+(Re)D_{\mathrm e}=-\frac12-E_+(R_{\mathrm e}) を求めよ。
  2. (2) (使う道具:2.1.2節・2.1.5節)(1) の DeD_{\mathrm e} が厳密な値より小さいことを定理 2.1.1 から説明せよ。また、ζ\zeta を加えて部分空間を広げると DeD_{\mathrm e} がどう変わるはずかを述べよ。
  3. (3) (使う道具:第1部第1章の Born–Oppenheimer 近似と、ばね定数を E0(R)E_0(R) の2階微分とみる見方)E+(R)E_+(R) の谷の曲率 K=d2E+/dR2K=\mathrm{d}^2E_+/\mathrm{d} R^2 を求め、換算質量 μ=Mp/2\mu=M_{\mathrm p}/2 を使って振動の量子 ℏω=ℏK/μ\hbar\omega=\hbar\sqrt{K/\mu} を eV で求めよ。厳密な曲率からも ℏω\hbar\omega を求めて比べよ。
  4. (4) (使う道具:第1部第1章のエネルギーの階層)厳密な ℏω\hbar\omega を使い、零点振動のエネルギー ℏω/2\hbar\omega/2 と DeD_{\mathrm e} の比、300 K300\,\mathrm K で振動の第1励起状態の Boltzmann 重み e−ℏω/(kBT)\mathrm{e}^{-\hbar\omega/( k_{\mathrm B} T)} を求め、室温で振動の自由度を凍らせてよいかを述べよ。また ℏω/Eh\hbar\omega/E_{\mathrm h} を κ2=(m/Mp)1/2\kappa^2=(m/M_{\mathrm p})^{1/2} と比べよ。
解答

(1) 式 (2.1.14) の E+(R)=−12+1R−C+X1+sE_+(R)=-\frac12+\frac1R-\frac{C+X}{1+s} を数値で最小にすると、Re=2.493a0=1.319 A˚R_{\mathrm e}=2.493a_0=1.319\,\text{\AA}、E+(Re)=−0.5648EhE_+(R_{\mathrm e})=-0.5648E_{\mathrm h} である。De=0.0648Eh=1.764 eVD_{\mathrm e}=0.0648E_{\mathrm h}=1.764\,\mathrm{eV} になる。

(2) 各 RR で E+(R)E_+(R) は厳密な偶の基底エネルギー E(R)E(R) 以上である(定理 2.1.1)。解離の極限 −12-\frac12 はどちらも同じなので、DeLCAO=−12−min⁡RE+(R)≤−12−min⁡RE(R)=DeD_{\mathrm e}^{\mathrm{LCAO}}=-\frac12-\min_RE_+(R)\le-\frac12-\min_RE(R)=D_{\mathrm e} が成り立つ。ζ\zeta を変分のパラメータに加えると部分空間が広がり、各 RR で上界は下がるか変わらないので、DeD_{\mathrm e} は大きくなるか変わらない。実際に最適化すると ζ=1.24\zeta=1.24 で De=2.35 eVD_{\mathrm e}=2.35\,\mathrm{eV} になる(2.1.5節)。

(3) ReR_{\mathrm e} のまわりで差分を取ると K=0.0628Eh/a02K=0.0628E_{\mathrm h}/a_0^2 である。μ=918.08m\mu=918.08m なので、原子単位で ω=0.0628/918.08=8.27×10−3\omega=\sqrt{0.0628/918.08}=8.27\times10^{-3}、ℏω=8.27×10−3Eh=0.225 eV\hbar\omega=8.27\times10^{-3}E_{\mathrm h}=0.225\,\mathrm{eV} になる。厳密な曲率 0.10300.1030 からは ω=0.1030/918.08=1.059×10−2\omega=\sqrt{0.1030/918.08}=1.059\times10^{-2}、ℏω=0.288 eV\hbar\omega=0.288\,\mathrm{eV} である。LCAO は谷を浅く、広く見積もるので、振動の量子を 22%22\% 小さく出す。曲率は SI で 0.1030Eh/a02=160 N m−10.1030E_{\mathrm h}/a_0^2=160\,\mathrm{N\,m^{-1}} で、これが陽子2個をつなぐばねのばね定数である。

(4) ℏω/2=0.144 eV\hbar\omega/2=0.144\,\mathrm{eV} は De=2.79 eVD_{\mathrm e}=2.79\,\mathrm{eV} の 5.2%5.2\% である。300 K300\,\mathrm K の kBT=0.0259 eVk_{\mathrm B} T=0.0259\,\mathrm{eV} に対して ℏω/(kBT)=11.1\hbar\omega/( k_{\mathrm B} T)=11.1、e−ℏω/(kBT)=1.4×10−5\mathrm{e}^{-\hbar\omega/( k_{\mathrm B} T)}=1.4\times10^{-5} で、室温では振動はほぼ全部が基底状態にあり、振動の自由度は凍らせてよい。κ2=(1/1836.15)1/2=0.0233\kappa^2=(1/1836.15)^{1/2}=0.0233 に対して ℏω/Eh=0.0106\hbar\omega/E_{\mathrm h}=0.0106 で、同じ桁に入る。第1部第1章で見積もった「振動のエネルギーは電子のエネルギーの κ2\kappa^2 倍の桁」がここでも成り立っている。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ