ガウシアンと群論 ( 名古屋工業大学 ) 川崎晋司
ガウシアンの特徴非経験的分子軌道計算 分子のシュレディンガー方程式をどう解くか HΨ = EΨ 電子だけでなく原子核も入る もちろん複数 一電子波動関数の形にして解こう = 分子軌道法 例えばハートリー法では多電子波動関数 Ψを一電子波動関数 φの積で近似 Ψ r 1, r, = ϕ r 1 ϕ r しかし この近似ではパウリの原理 ( 電子の入れ替えに反対称 ) を満足しない スレーター行列を使う (HF 法 ) さらに電子相関も考慮して (MP 法 CI 法 )
入力 出力ファイル 入力 *.gif ファイル *.chk ファイル 出力 *.log ファイル バイナリ テキスト
入力ファイル なし : 一点計算 opt: 構造最適化 Freq: 振動計算 %chk=c: Work Kawasaki CH4.chk # opt mp/6-31g geom=connectivity pop=full 分子軌道の成分などを出力するときに必要 (full にしておけばすべての情報が出て良いと思う ) GaussView からは additional Keywords 欄に入れる Title Card Required 0 1 C 0.00000000 0.00000000 0.00000000 H 0.647544 0.647544 0.647544 H -0.647544-0.647544 0.647544 H -0.647544 0.647544-0.647544 H 0.647544-0.647544-0.647544 計算方法 基底関数 1 1.0 3 1.0 4 1.0 5 1.0 3 4 5
Method 電子密度汎関数 DFT とそれ以外に大別 コーンシャム方程式 ハートリーフォック方程式 1 Z A r A + ρ r dr + VXC(r) ψ r = εψ(r) r r 1 Z A r A + J K ψ r = εψ(r) レベルは HF << MP < CISD MP4 < CCSD DFT
DFT ( 分子軌道計算の場合 ) B3LYP 交換汎関数 相関汎関数 普通はこれで良いようだ 励起状態を扱う場合は TDDFT が良いらしい 通常利用される汎関数では離れた電子間の相関をうまく取り込めないらしい ( 分散力の評価がうまくない ) van der Waals 力の評価に問題あり
Basis Set 多電子原子の電子の軌道 : 球面調和関数は水素と同様 動径関数の部分は近似関数としてスレータ型原子軌道 (STO) がある R r = c ξ a0 n 1 r n 1 e ξ a 0 r この STO は分子軌道計算を行う上で時間を食うのでガウシアンではこの STO をいくつかのガウス関数 (GTO) で近似する http://www.chem.ous.ac.jp/~waka/compchem /general_atom/ga-7.html 1 つの STO を複数の GTO で表したものを contracted GTO= CGTO という
DZ (Double-Zeta) 1 つの AO に つの CGTO をあてたもの 3-1G 内殻軌道を 3 つのガウス関数で表す 価電子は つの CGTO を使う (DZ) 内訳は つの G でできた CGTO と 1 つの G でできた CGTO TZ (Triple-Zeta) 6-311G 内殻軌道 :6 つのガウス関数価電子軌道 :3 つの CGTO 内訳は 3 つの G でできた CGTO と 1 つの G でできた CGTO が つ
分極関数と分散関数 (Polarization) (Diffuse) 大きいアニオンのように広がったもの + をつけて示す 例 ) 6-31G 6-31+G 球対称からのずれを軌道を足して回避 DZ DZP 例 ) 6-31G 6-31G(d p)=6-31g** 水素以外は d 軌道を足す 水素には p 軌道を足す
MP DFT-B3LYP なら DZP 以上 MP4 CCSD なら TZP 以上 らしい 遷移元素 cc-pvdz cc-pvtz VZP TZP
ガウシアンで何を求める? 1. ( 熱力学 ) エネルギー. 分子軌道 ( 各原子からの寄与 ) 3. 振動 (IR Raman)
振動計算 Freq=Raman で IR と Raman 両方計算 (additional keywords の欄に入力 ) つの目的 (1) 構造の正しさ評価 求めた固有振動数に負のものがあればダメ ( ただし 励起状態では 1 つだけ負 ) () 振動スペクトル求める (IR, Raman)
調和振動子近似 振動が激しくなっても平均位置は変わらない 非調和項が入ると
構造最適化はポテンシャルミニマムを探すべく 微分係数がゼロとなるところを探すのですが このケースも微係数ゼロ V=1/(-k)x 平方根が振動数 : 虚数
振動計算 Result Vibration 赤外スペクトル ラマンスペクトル 振動数選ぶとグラフにマークでる 選択した振動数のモードが動く
IR スペクトル Raman スペクトル Epsilon Obs. Calc. (0.96*x) 500 1000 1500 000 Frequency / cm 1 Intensity / arb. units obs. calc. 1000 1500 000 Raman Shift / cm 1
pop=full (additional keywords 欄に入力 ) 分子軌道の成り立ちを調べる O: 占有軌道 V: 非占有軌道
構造最適化と対称性
混成軌道? sp, sp 3 と言いますが sp pz sp3 1s 1s 果たして等価な 3 つの sp や等価な 4 つの sp3 軌道は見えるのか?
8=HOMO 7 6 4 5 3
三中心四電子結合 He H +
1 =HOMO 結合性軌道 非結合性軌道 3=LUMO 反結合性軌道
中松博英 向山毅 XANES と分子軌道法 放射光 7 89 (1994)
群論 目標 : マリケンの記号の意味を少しだけ理解する
注意! 本稿を通じて例として取り上げている C3v について : y 左図のように x 軸上に H 原子がのっているように x, y 軸を取りました しかし どうも一般的には y 軸上に H 原子をのせるようです x 軸の取り方によってもちろん結論は変わりませんが行列 ( 表現 ) は変わってきます 他書との比較ではその点を気をつけてください
NH 3 y N を含む z 軸が主軸 (3 回軸 ) x C 3v σ v '' σ v ' 3 枚の主軸を含む鏡面 σ v
対称操作 ( 群の要素 ) C 3v には以下の要素がある 1 恒等操作 E( すべての群にある ) 何もしないという操作 3 回軸 10 回す C 3 と 40 回す C 3 の つ 3 鏡面 ( 主軸を含む ) 前頁のように σ v, σ v ', σ v '' E, C 3, C 3, σ v, σ v ', σ v '' の 6 つ
対称操作 ( 要素 ) の行列表現 対称操作による (x,y,z) 座標の変換を例として考えると 恒等操作 E は 表現 基底 同様に 10 回す C3 は以下のように表せる
C 3v の表現 1 0 0 0 1 0 0 0 1 1 3 0 3 1 0 0 0 1 E C 3 C 3 1 3 0 3 1 0 0 0 1 σ v 1 0 0 0 1 0 0 0 1 σ v ' 1 3 0 3 1 0 0 0 1 σ v '' 1 3 0 3 1 0 0 0 1
要素の掛け算 ( 対称操作 回 ) 10 40 10 + C 3 = C 3 C 3
掛け算表 ( 積表 ) E C 3 C 3 σ v ' σ v E E C 3 C 3 σ v ' σ v '' σ v '' σ v '' C 3 C 3 C 3 E σ v σ ' v σ v 先 C 3 C 3 E C 3 σ v ' '' σ v σ v σ v σ ' '' v σ v σ v E C 3 C 3 ' ' σ v σ v σ v '' '' σ v σ v C 3 E C 3 '' σ v σ v ' σ v C 3 C 3 E 後
表現の掛け算 C 3 = C 3 C 3 を表現行列でも確認
大きな行列表現 z 3 z 4 y x y z x y x z y x 各原子の座標変位ベクトルを基底とする表現を考える 1
大きな行列表現 座標変位ベクトルに対する C3 表現 x1 y1 z1 x y z x3 y3 z3 x4 y4 z4 = 0 0 0 1 0 0 0 3 1 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 3 1 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 3 1 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 3 1 3 0 0 x1 y1 z1 x y z x3 y3 z3 x4 y4 z4
相似変換 行列 A に対してある行列 X とその逆行列 X -1 をつぎのように作用させることを相似変換という A = X 1 A X もし A, B, C. がある群の要素である時 A, B, C. もまた その群の表現となる
表現行列の区画因数化 大きな行列においてうまい相似変換により次のような表現を見つけることができることがある 0 可約表現 もとの大きい行列 0 規約表現
指標 行列の対角成分を足し算したものを指標という E 1 0 0 0 1 0 0 0 1 この行列では 1+1+1 = 3 が指標 C 3 1 3 0 3 1 0 0 0 1 この行列では -1/-1/+1 = 0 が指標
指標表 C 3v の指標表 類 ( 要素 ) E C 3 3σ v A 1 1 1 1 z x +y, z A 1 1-1 R z E -1 0 (x,y) (R x, R y ) (x -y, xy) (xz, yz) 基底の例 既約表現の組 指標
既約表現の記号 一次元の表現は A か B 二次元は E, 三次元は T 1. 主軸の C n 軸の回転に対して対称的 A 反対称的 B. 添え字の 1, は主軸に垂直な C 軸 ( または σd) に対して対称的 1 反対称的 3. プライム 二重プライムは σh に対する対称性 4. 対称心がある群では反転に対して対称 g 反対称 u
指標表に関わるルール 1. 既約表現の数は類 ( 要素の固まりの数 ) の数と同じ. 表現の次元の二乗和が位数 h( 要素の数 ) に等しい 3. 既約表現の指標の二乗和は位数 h( 要素の数 ) に等しい 4. つの規約表現の指標のベクトルは直交
指標表の決め方 C 3v : E, C 3, 3σ v 類の数 3 規約表現は 3 組 規約表現は 3 組 表現の次元を l 1, l, l 3 とする 位数は 6 l 1 + l + l 3 = h = 6 l 1 = l = 1 l 3 = (1 次元 つと 次元 1つ ) Γ1 Γ Γ3 E C 3 3σ v 1 1 1 a 1 b c d 次元と同じになる どの群にも恒等表現がある 1 + a + 3 c =6 + c + 3 d =6 連立して解く (1,1,1) (1,a,b) = 0, (1,1,1) (,c,d) = 0, (1,a,b) (,c,d) = 0 前頁の指標表
指標表の使い方 1 - 振動モードの解析 - IR やラマン活性モードを調べる 4 z y x z 各原子の座標変位ベクトルを基底とする 3 z y x z y x y x 1
指標を数える E 1 0 0 1 1 1 なので対角成分に 1 個 1 が並ぶ指標 χ(e) = 1 C 3 σ v については つの原子 (1,4) のみが対角成分に関わる 1-1 1 0 0 0 0 0 0 1-1 1 x1 y1 z1 x y z x3 y3 z3 x4 y4 z4 0 0 0 1 0 0 0 3 1 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 3 1 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 3 1 3 指標 χ(c3) = 0 とわかる 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 指標 χ(σ v ) = とわかる 3 1 3 0 0 x1 y1 z1 x y z x3 y3 z3 x4 y4 z4
NH3 E C 3 3σ v A 1 1 1 1 z x +y, z A 1 1-1 R z E -1 0 (x,y) (R x, R y ) (x -y, xy) (xz, yz) Γ 振動 1 0 この Γ 振動の中にどんな規約表現があるかを数えます これには便利な規則があって次のようにして数えます 位数 a i = 1 χ(r)χ h i (R) 指標表の指標 a 1 = 1 6 a = 1 6 a 3 = 1 6 1 1 1 + 1 0 + 3(1)() = 3 可約表現の指標 1 1 1 + 1 0 + 3( 1)() = 1 1 1 1 + 1 0 + 3(0)() = 4 Γ 振動 =3A 1 +A +4E が得られた
NH3 Γ 振動 =3A 1 +A +4E とわかりましたが これは多すぎます なぜなら非直線型の 4 原子分子の基準振動の数は自由度 3 4 原子 6 = 6 個だからです 6 を引いていますが これは分子全体の並進と回転です すべての原子が x y z の同じ方向に動くときなので 指標表で x y z が基底となるものを探す A1 + E 同じ理屈で 指標表で Rx Ry Rz が基底となるものを探す A + E これを引くと Γ 振動 =A 1 +E
赤外活性とラマン活性 例外もあるので注意! 赤外活性 : 励起される振動がデカルト座標と同じ表現に属する 指標表で基底が x, y, z のものを探す ラマン活性 : 振動が分極率テンソル成分と同じ表現に属する x, y, z, xy, yz, zx のものを探す
NH3 基準振動は Γ 振動 =A 1 +E とわかった E C 3 3σ v A 1 1 1 1 z x +y, z A 1 1-1 R z E -1 0 (x,y) (R x, R y ) (x -y, xy) (xz, yz) 前頁のルールから 赤外活性 :A1, E ラマン活性 :A1, E
CO 3 - D 3h コットンの 群論の化学への応用 の例そのまま E C 3 3C σ h S 3 3σ v A 1 1 1 1 1 1 1 x +y, z A 1 1-1 1 1-1 Rz E -1 0-1 0 (x,y) (x -y, xy ) A 1 1 1 1-1 -1-1 A 1 1-1 -1-1 1 z E -1 0-1 0 (Rx,Ry) (xz,yz) Γ 振動 1 0-4 - Γ 振動 = A 1 + A + 3E + A + E 並進と回転をとると Γ 振動 = A 1 + E + A 赤外活性 : E + A ラマン活性 :A 1 + E
BF 3 (D 3h ) 分子の振動計算 ( ガウシアン ) A 1 A ラマン活性 赤外活性 E 赤外活性 ラマン活性 E 赤外活性 ラマン活性
BF3 分子の振動計算 ( ガウシアン ) E A E
指標表の使い方 - 錯体の d 電子の分裂 - O h T d 八面体配位 四面体配位
C 3 C 3 4 1 5 1 3 3 1 4 5 5 6 6 4 6 指標はゼロ指標は 1 4 5 3 3 4 1 5 6 6
O h 中心金属イオンの電子について考える E 8C 3 6C 6C 4 3C i 6S 4 8S 6 3σ h 6σ d A 1g 1 1 1 1 1 1 1 1 1 1 x +y +z A g 1 1-1 -1 1 1-1 1 1-1 E g -1 0 0 0-1 0 (z -x -y, x - y ) T 1g 3 0-1 1-1 3-1 0-1 -1 (Rx,Ry,Rz) T g 3 0 1-1 -1 3-1 0-1 1 (xz,yz,xy) A 1u 1 1 1 1 1-1 -1-1 -1-1 A u 1 1-1 -1 1-1 1-1 -1 1 E u -1 0 0-0 1-0 T 1u 3 0-1 1-1 -3-1 0 1 1 (x,y,z) T u 3 0 1-1 -1-3 1 0 1-1 指標表から s : A 1g p : T 1u dxy, dxz, dyz: T g dz, dx -y : E g
T d E 8C 3 3C 6S 4 6σ d A 1 1 1 1 1 1 x +y +z A 1 1 1-1 -1 E -1 0 0 (z -x -y, x - y ) T 1 3 0-1 1-1 (Rx,Ry,Rz) T 3 0-1 -1 1 (x,y,z) (xz,yz,xy) s : A 1 p : T dxy, dxz, dyz: T dz, dx -y : E
配位子の電子 ( ただし s と p だけ ) について考える z y x 配位子の s 電子と px 電子は σ 結合を形成する π 結合について考える
配位子の s を基底として指標をとる (px を基底としてとっても同じ ) 対称操作によって原子の位置が動いてしまうようなものは表現の対角成分が 0 になることに注意すれば意外に簡単に指標を取れる 例えば反転対称があるもの (i, S 4, S 6 ) や原子を通らない回転軸 (8C 3, 6C ( 八面体の辺の中心を通る )) の場合は直ちに指標がゼロとわかる Γ s(px) 6 0 0 0 0 0 4
E 8C 3 6C 6C 4 3C i 6S 4 8S 6 3σ h 6σ d A 1g 1 1 1 1 1 1 1 1 1 1 x +y +z A g 1 1-1 -1 1 1-1 1 1-1 E g -1 0 0 0-1 0 (z -x -y, x - y ) T 1g 3 0-1 1-1 3-1 0-1 -1 (Rx,Ry,Rz) T g 3 0 1-1 -1 3-1 0-1 1 (xz,yz,xy) A 1u 1 1 1 1 1-1 -1-1 -1-1 A u 1 1-1 -1 1-1 1-1 -1 1 E u -1 0 0-0 1-0 T 1u 3 0-1 1-1 -3-1 0 1 1 (x,y,z) T u 3 0 1-1 -1-3 1 0 1-1 Γ s(px) 6 0 0 0 0 0 4 a i = 1 h χ(r)χ i (R) でひたすら数え上げると a 1g = 1 48 e g = 1 48 t 1u = 1 48 1 1 6 + 6 1 + 3(1)() + 3 1 4 + 6(1)() = 1 1 6 + 6 0 + 3()() + 3 4 + 6(0)() = 1 1 3 6 + 6 1 + 3( 1)() + 3 1 4 + 6(1)() = 1 他はすべてゼロで Γ s = A 1g + E g + T 1u とわかる
配位子の py, pz を基底として指標をとる 振動座標のときと同様 対称操作で動かない原子 にだけ注目すればよい E 8C 3 6C 6C 4 3C i 6S 4 8S 6 3σ h 6σ d A 1g 1 1 1 1 1 1 1 1 1 1 x +y +z A g 1 1-1 -1 1 1-1 1 1-1 E g -1 0 0 0-1 0 (z -x -y, x - y ) T 1g 3 0-1 1-1 3-1 0-1 -1 (Rx,Ry,Rz) T g 3 0 1-1 -1 3-1 0-1 1 (xz,yz,xy) A 1u 1 1 1 1 1-1 -1-1 -1-1 A u 1 1-1 -1 1-1 1-1 -1 1 E u -1 0 0-0 1-0 T 1u 3 0-1 1-1 -3-1 0 1 1 (x,y,z) T u 3 0 1-1 -1-3 1 0 1-1 Γ py,pz 1 0 0 0-4 0 0 0 0 0 Γ py,pz = T 1g + T g + T 1u + T u 中心金属イオンの T g (dxy dxz dyz) と π 結合形成 T 1u の p 電子は σ 結合に使用される
C Z y z z y y Y z
C Z y z z y y Y z