三重大学工学部

Transcription

1 反応理論化学 ( その5 6 ポテンシャルエネルギー面と反応経路最も簡単な反応 X + Y X + Y 反応物 ( 生成物 (P X 結合が切断反応系全体のエネルギーは X と Y の Y 結合が形成原子間距離によって変化 r(x と r( Y に対してエネルギーを等高線で表す赤矢印 P:X 結合の切断と Y 結合の形成が同時進行青矢印 P: まず X 結合が切断し次いで Y 結合が形成 谷 X + + Y r(-y 峠 山 谷 P エネルギー X Y P r(x- 反応経路 P の経路は無数にあるが P を通る経路が最も容易である : 反応経路に沿ってエネルギーの最も高い状態 = 遷移状態 (ranston State 反応は遷移状態に到るエネルギーを必要とする活性化エネルギー : E( E ( E P E 反応熱 : ( ( 協奏反応は遷移状態が つ存在する 段階反応例えば置換反応 H HO - + H 5 H Br HO H H 5 H Br - H HO + Br - H 5 H 遷移状態 Dels-lder 反応 H H H H + H H H H H H H H H H H H H H 遷移状態

2 反応経路の途中にエネルギーの窪みがあれば反応中間体が存在する多段階反応例えば反応物 遷移状態 中間体 遷移状態 生成物 ( (S (I (S (P S S エネルギー I P 反応経路 ポテンシャルエネルギー面の立体図と等高線図 : 分子理論の展開永瀬茂 平尾公彦岩波書店 (00 山 谷 谷 安定構造 反応経路 : つの谷底である反応物から出発して別の谷底である生成物へ行くのに最もエネルギーの低い道筋 安定構造 遷移状態 極小 極大 鞍点 ( 遷移状態 : 反応経路に沿ってエネルギーは極大他の方向に関してエネルギーは極小 反応経路 反応熱 : 反応物 ( 出発点 と生成物 ( 到着点 のエネルギー差反応物 > 生成物 発熱反応反応物 < 生成物 吸熱反応

3 化学反応の定量的な予測 : 安定構造 遷移状態 活性化エネルギー 反応経路などポテンシャルエネルギー面 ( 構造変化とエネルギー変化の対応関係 の情報が不可欠 構造変化は 6個の自由度 ( 直線分子は 5個の自由度 をもつ多原子分子では全ての可能な構造に対してエネルギー計算を実行することはできない 化学反応において重要なのは安定構造と遷移状態およびそれらを結ぶ反応経路 量子化学計算により高精度で定常点 ( 安定構造 遷移状態 および反応経路の探索ができるようになってきている=エネルギー勾配法 ( エネルギー微分法 および基準振動解析 7 エネルギー勾配法 分子の運動の自由度 原子当たり xyz,, 方向の3つの自由度をもつ N 原子分子では の自由度がある xyz,, 方向の並進運動の自由度が3 ( 分子の構造は変化しない xyz,, 軸まわりの回転の自由度が3 ( 分子の構造は変化しない 非直線分子では振動の自由度は 6( 直線分子では振動の自由度は 5 メッシュ法によるポテンシャルエネルギー面上の定常点の探索直線 原子分子 : 振動の自由度 =3 5= 原子間距離を変えてエネルギー計算 E 非直線 3 原子分子 : 振動の自由度 =3 3 6= 3 つの原子間距離とつの結合角を独立に変化させてエネルギー計算 自由度当り0 点だけ考慮しても計 回のエネルギー計算が必要非直線 N 原子分子 : 振動の自由度 = 6 メッシュ法により可能な全ての構造をエネルギー計算することは非現実的 核座標に関するエネルギー微分を計算し定常点を自動的に探索ポテンシャルエネルギー面の構造を全て計算する必要はない ポテンシャルエネルギー面上の谷底や峠は定常点 ( 極小点 極大点 = { } 原子核の空間座標 = { X Y Z X Y Z XN YN ZN} 定常点 e では E( e = 0 for =,,, 核座標 に対するエネルギー勾配 (7- E( e F ( e = = 0 for =,,, 核座標 に働く力 (7-3

4 定常点では 全ての核座標に対するエネルギー勾配がゼロ = 全ての原子核に働く力がゼロ全ての変位方向に対して極小 安定構造 ( 谷底 つの変位方向に対して極大で他の変位方向に対しては極小 遷移状態 ( 峠 エネルギー勾配法 ある構造 の近傍の構造 のエネルギーは のまわりで 次まで aylor 展開すると E E E( = E( + ( + ( ( = j j = j= j 3 3 = ( ( N N E ( ( + G + Hj j j = = j= (7-3 E G = グラジエント ( 勾配 (7-4 E Hj = ヘシアン ( 力の定数 (7-5 j (7-3 の両辺を k で偏微分 E( N ( ( N N = E ( ( + G + Hj j j k k = = j= E( j = + G + Hj ( j j + Hj ( k = k = j= k = j= k = G δ k + H j ( δ ( k j j + Hj δjk = = j= = j= = Gk + Hkj ( j j + Hk ( (7-6 j= = = Gk + Hk ( + Hk ( ( = = k k ( = = G + H E ( 3 = G + N Hj j j j= ある構造 から だけ離れた近傍に定常点の構造 e があるとする = e (7-7 (7- より (7-6 で を e とすると E( e e = G + Hj ( j j = G + Hj j = 0 j= j= for =,,, (7-8 4

5 行列表記を用いると 核座標の列ベクトル = (7-9 3 N Δ 変位の列ベクトル e e Δ = = (7-0 e 3 N Δ 変位の列ベクトルを転置した行ベクトル Δ = ( (7- G グラジエント ( 勾配 の列ベクトル E E G = (7- E 3 N H ヘシアン ( 力の定数 の行列 E E E 3 3 N E E E H = (7-3 E E E N N N (7-3 は E( e = E ( + Δ G + Δ HΔ (7-4 (7-8 は G + H Δ = 0 Δ = H G (7-5 よって = + Δ = H G (7-6 e ある構造 でのエネルギー勾配 G ( 次微分 と力の定数行列 H ( 次微分 がわかれば 定常点の構造 e への方向ベクトル Δ が計算できる 多原子分子に対して定常点の構造 を予測できる e 5

6 実際には適当な初期構造から出発して繰り返し計算を行う 構造最適化初期構造 (0 を設定 H G を決定 ( 構造最適化の 回目の構造 = ( ( 0 ( 0 ( 0 = H G を決定 ( 構造最適化の 回目の構造 ( ( ( ( = H G ( n ( n ( n ( n を決定 ( 構造最適化の n 回目の構造 定常点 ( 全ての原子核に働く力がゼロ に到達するまで探索を繰り返す初期構造の設定が悪いと繰り返し計算の回数が多くなる 初期構造 エネルギー 極小 困難な場合もあり得る 極小 最安定 通常 初期構造に最も近い安定構造が得られる構造異性体など複数の安定構造が存在する場合はそれぞれの異性体に対応した初期構造を設定する必要がある 6

7 構造最適化の計算例 : 分子理論の展開永瀬茂 平尾公彦岩波書店 (00 HO H + O 反応 7

8 エネルギー微分の実用的な応用ポテンシャルエネルギー面の定常点の探索 E E 双極子モーメント F ( F : E 電場, 分極率 F 3 E 赤外スペクトルの強度 F, ラマンスペクトルの強度 E F E NM のケミカルシフト ( B : 磁場, I : 磁気モーメント BI 探索した定常点がエネルギーの極小点か極大点か定常点 e のエネルギー E ( を = e e E ( を e に関して 次まで aylor 展開すると Δ だけ離れた近傍 のエネルギー E ( と比較する E( = E ( e + Δ Ge + Δ HeΔ (7-7 定常点ではG e = 0なので第 3 項が正 E( e < E ( : e は極小点第 3 項が負 E( e > E ( : e は極大点第 3 項について対称行列は直交行列を用いて対角化できる ( 力の定数行列を対角化 = k δ (7-8 ( e UHU j j He を対角化する直交行列 U を用いて座標変換する ( U = U Q= U Δ = U Δ Δ = UQ (7-9 Δ HeΔ = ( UQ He( UQ = Q U HeUQ = Q ( U HeU Q k 0 0 Q (7-0 = ( Q Q 0 0 = kq = 0 0 k Q Q は正であるので力の定数行列 He の固有値 k の正負だけで判別できる k が全て正 全ての方向への変位 Q に対してエネルギー上昇 ( 極小点 安定構造 k のうち k だけ負 Q 方向への変位に対してのみエネルギー低下 ( 極大点 他方向への変位に対してはエネルギー上昇 ( 極小点 遷移状態 Q : 遷移ベクトル ( 遷移状態での反応座標 = 反応の進行方向 k : 力の定数, Q : 基準座標 ( 振動様式を与える 3 6 ν = k : 基準振動数 ( 赤外スペクトル, N h ν : ゼロ点振動エネルギー π = 8

9 8 反応経路の決定 ポテンシャルエネルギー面の谷底 ( 反応物 中間体 生成物 と峠 ( 遷移状態 はエネルギー勾配法により見つけることができる 反応物 遷移状態 中間体 遷移状態 生成物を滑らかに結ぶ道筋 ( 反応経路 を知りたい 遷移状態では力の定数行列の固有値がただ つだけ負になる 虚数の振動数をもつ基準座標は峠での反応座標 ( エネルギーの低下する方向 を示している 固有反応座標 (I = Intrnsc eacton oordnate 反応座標の方向にゆっくりと峠を下っていくと谷底に到達する道筋を一義的に見つけることができる 峠からボールを道筋からはずれないようにきわめて静かに転がしたときの軌跡 峠 谷底 谷底 峠を出発点として古典的運動方程式 ( 質量 加速度 = 力 を全ての核座標 に対して解く = = { } { X Y Z X Y Z X Y Z } d E M = d t N N N 短い時間範囲 δt で積分 δ d d d d 左辺 : t M d = 0 t M dt dt dt δ t dt 0 E δ 右辺 : d ( δ 0 = = δ t E E t t t 0 d d E M = δt dt δ t dt 0 最初は核が静止していると仮定 d = 0 dt 0 (8- よりδt 時間後には d E M = δt dt Md δ tdt = E (8- (8- (8-3 (8-4 (8-5 9

10 全ての原子核に対して上式が成立 Md Md Md = = = E E E 3 N この関係を満足する d,d,,d の比で分子が変形する 反応座標の方向質量荷重座標を用いると η = M dη dη dη = = = E E E η η η (8-6 (8-7 (8-8 d 無限小時間ごとに核速度 をゼロにして計算される古典的軌跡 = I dt 反応系の核運動を記述する平均の軌跡 = 質量荷重座標系では峠からの最急降下経路 ポテンシャルエネルギー面の等高線に常に垂直な曲線 谷底 I 古典的軌跡 峠 谷底 I における構造変化ホルムアルデヒド ヒドロキシカルベンにおける,- 水素原子シフトの固有反応座標計算 (I から得られる一連の構造 ホルムアルデヒド 遷移状態 ヒドロキシカルベン 0